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The field of laser-matter interaction traditionally deals with the response of atoms, 
molecules and plasmas to an external light wave. However, the recent sustained tech- 
nological progress is opening up the possibility of employing intense laser radiation to 
trigger or substantially influence physical processes beyond atomic-physics energy scales. 
Available optical laser intensities exceeding 10^^ W/cm^ can push the fundamental light- 
electron interaction to the extreme limit where radiation-reaction effects dominate the 
electron dynamics, can shed light on the structure of the quantum vacuum, and can 
trigger the creation of particles like electrons, muons and pions and their corresponding 
antiparticles. Also, novel sources of intense coherent high-energy photons and laser- 
based particle colliders can pave the way to nuclear quantum optics and may even allow 
for potential discovery of new particles beyond the Standard Model. These are the 
main topics of the present article, which is devoted to a review of recent investigations 
on high-energy processes within the realm of relativistic quantum dynamics, quantum 
electrodynamics, nuclear and particle physics, occurring in extremely intense laser fields. 
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I. INTRODUCTION 

The first realization of the laser in 1960 (Maiman, 
1960) is one of the most important technological break- 
throughs. Nowadays lasers are indispensable tools for 
investigating physical processes in different areas rang- 
ing from atomic and plasma physics to nuclear and high- 
energy physics. This has been possible mainly due to the 
continuous progress made along two specific directions: 
decrease of the laser pulse duration and increase of the 
laser peak intensity (Mourou and Tajima, 2011). On the 
one hand, multiterawatt laser systems with a pulse du- 
ration in the femtosecond time scale are readily available 
nowadays and different laboratories have succeeded in the 
generation of single attosecond pulses. Physics at the at- 
tosecond time scale has been the subject of the recent 
review Krausz and Ivanov, 2009. In this review it has 
been pointed out how pulses in the attosecond domain 
allow for the detailed investigation of the electron mo- 
tion in atoms and during molecular reactions. The pro- 
duction of ultrashort pulses is strongly connected with 
the increase of the laser peak intensity. This is not only 
because temporal compression evidently implies an in- 
crease in intensity at a given laser energy, but also be- 
cause higher intensities allow, in general, for controlling 
faster physical processes, which in turn can be exploited 
for generating correspondingly shorter light pulses. 

Not long after the invention of the laser, available in- 
tensities were already sufficiently high to trigger non- 
linear optical effects like second harmonic generation. 
It is, however, only after the experimental implementa- 
tion of the Chirped Pulse Amplification (CPA) technique 
(Strickland and Mourou, 1985) that it has been possi- 
ble to reach the intensity threshold of 10^^-10^^ W/cm^ 
corresponding to electric field amplitudes of the same 
order as the Coulomb field in atoms. At such inten- 
sities the interplay between the laser and the atomic 
field significantly alters the electron's dynamics in atoms 
and molecules and this can be exploited, for example, 
for generating high-frequency radiation in the extreme- 
ultraviolet (XUV) and soft-x-ray regions (high-order 
harmonic generation or HHG) (Agostini and DiMauro, 
2004; Protopapas et ai, 1997). HHG as well as atomic 
processes in intense laser fields have been recently re- 
viewed in Fennel et ai, 2010; Teubner and Gibbon, 2009; 
and Winterfeldt et al, 2008, with specific emphasis on 



the control of high-harmonic spectra by spatio-temporal 
shaping of the driving pulse (Winterfeldt et al., 2008), on 
harmonic generation in laser-plasma interaction (Teub- 
ner and Gibbon, 2009) and on the dynamics of clusters 
in strong laser fields (Fennel et ai, 2010). 

By increasing the optical laser intensity to the or- 
der of lO^'^-lO^^ W/cm^, another physically important 
regime in laser-matter interaction is entered: the rel- 
ativistic regime. In such intense electromagnetic fields 
an electron reaches relativistic velocities already within 
one laser period, the magnetic component of the Lorentz 
force becomes of the same order of magnitude of the elec- 
tric one, and the electron's motion becomes highly non- 
linear as a function of the laser's electromagnetic field. 
Although the increasing influence of the magnetic force 
causes a suppression of atomic HHG in the relativistic do- 
main, the highly nonlinear motion of the electrons in such 
strong laser fields is at the origin of numerous new effects 
as relativistic self-focusing in plasma and laser wakefield 
acceleration (Mulser and Bauer, 2010). In the recent re- 
views Ehlotzky et at, 2009; Mourou et ai, 2006; and 
Salamin et al., 2006, different processes occurring at rel- 
ativistic laser intensities are discussed. In particular, in 
Ehlotzky et ai, 2009, QED processes like Compton, Mott 
and M0ller scattering in a strong laser field are covered, in 
Mourou et al, 2006, technical aspects and new possibili- 
ties of the CPA techniques are reviewed together with rel- 
ativistic effects in laser-plasma interaction as, for exam- 
ple, self-induced transparency and wakefield generation, 
while in Salamin et ai, 2006, spin-effects as well as rel- 
ativistic multiphoton and tunneling recollision dynamics 
in laser-atom interactions are reviewed. Also in the same 
year another review was published on nonlinear collective 
photon interactions, including vacuum-polarization ef- 
fects in a plasma (Marklund and Shukla, 2006). Whereas, 
the physics of plasma-based laser electron accelerators 
is the main subject covered in Esarey et al., 2009 and 
Malka, 2011, with a special focus on the different phases 
involved (electron injection and trapping, and pulse prop- 
agation) and on the role of plasma instabilities in the 
acceleration process. Finally, in Ruffini et ai, 2010 dif- 
ferent processes related to electron-positron (e+-e^) pair 
production are reviewed with special emphasis on those 
occurring in the presence of highly-charged ions and in 
astrophysical environments. 

In the present article we address physical processes 
that mainly occur at optical laser intensities mostly larger 
than 10^^ W/cm^, i.e., well exceeding the relativistic 
threshold. After reporting on the latest technological 
progress in optical and x-ray laser technology (Sec. II), 
we review some basic results on the classical and quan- 
tum dynamics of an electron in a laser field (Sec. III). 
Then, we bridge to lower-intensity physics by reviewing 
more recent advances in relativistic ionization and HHG 
in atomic gases (Sec. IV). The main subject of the re- 
view, i.e., the response of fundamental systems like elec- 



3 



trons, photons and even the vacuum to uhra-intense radi- 
ation fields is covered in Sees. V-X. As wiU be seen, such 
high laser intensities represent a unique tool to investi- 
gate fundamental processes like multiphoton Compton 
scattering (Sec. V), to clarify conceptual issues like radi- 
ation reaction in classical and quantum electrodynamics 
(Sec. VI) and to investigate the structure of the quantum 
vacuum (Sec. VII). Also, other fundamental quantum- 
relativistic phenomena like the transformation of pure 
light into massive particles as electrons, muons and pi- 
ous (and their corresponding antiparticles) can become 
feasible and can even limit the attainability of arbitrar- 
ily high laser intensities (Sees. VIII-X). Finally, we will 
also review recent suggestions on employing novel high- 
frequency lasers and laser-accelerated particle beams to 
directly trigger nuclear and high-energy processes (Sees. 
XI- XIII). The main conclusions of the article will be pre- 
sented in Sec. XIV. 

Units with h = c = 1 and the space-time metric 
TyA"' = diag(-f 1, — 1, — 1, — 1) are employed throughout 
this review. 



II. NOVEL RADIATION SOURCES 

In this section we review the latest technical and ex- 
perimental progress in laser technology. We will discuss 
optical and x-ray laser systems separately. The latter 
are especially useful for e^-e" pair production, for di- 
rect laser-nucleus interaction, as well as as probes for 
vacuum-polarization effects (see in particular Sees. VII, 
VIII and XI). For overviews of feasible accelerators also 
of relevance for the present review see, e.g., Esarey et al, 
2009; Malka, 2011; Nakamura et al, 2010; and Wilson, 
2001 and the relevant original literature as quoted in the 
respective sections. 



A. Strong optical laser sources 

As has been mentioned in the Introduction, since the 
invention of the CPA technique (Strickland and Mourou, 
1985) laser peak intensities have been boosted by sev- 
eral orders of magnitude. Another amplification tech- 
nique called Optical Parametric Chirped Pulse Ampli- 
fication (OPCPA), based on the nonlinear interaction 
among laser beams in crystals, was suggested almost at 
the same time as the CPA and proved to be promising as 
well (Piskarskas et al., 1986). As a result of the increase 
in available laser intensities, exciting perspectives have 
been envisaged in different fields spanning from atomic to 
plasma and even nuclear and high-energy physics (Feder, 
2010; Gerstner, 2007; Mourou, 2010; Tajima et al, 2010). 

The group of G. Mourou at the University of Michi- 
gan (Michigan, USA) holds the record so far for the 
highest laser intensity ever achieved of 2 x 10^^ W/cm^ 



(Yanovsky et al, 2008), while no experiments have been 
performed at this intensity yet. This record intensity has 
been reached when the HERCULES laser was upgraded 
to become a 300 TW Ti:Sa system, amplified via CPA 
and capable of a repetition rate of 0.1 Hz. The 4-grating 
compressor allowed for a pulse duration of about 30 fs 
and adaptive optics together with a //I parabola en- 
abled to focus the beam down to a diameter of about 
1.3 /im. This experimental achievement on the laser in- 
tensity pushed the capabilities of a multiterawatt laser 
almost to the limit. 

The 1-PW threshold has been already reached and 
even exceeded in various laboratories. For example, the 
Texas Petawatt Laser (TPL) at the University of Texas at 
Austin (Texas, USA) has exceeded the Petawatt thresh- 
old thanks to the OPCPA technique, by compressing an 
energy of 186 J in a pulse lasting only 167 fs (TPL, 2011). 
The TPL has been employed for investigating laser- 
plasma interactions at extreme conditions, particularly 
relevant for astrophysics. Also, the two laser systems 
Vulcan (Vulcan, 2011) and Astra Gemini (Astra Gemini, 
2011) at the Central Laser Facility (CLE) in the United 
Kingdom provide powers of the order of 1 PW. The 
Vulcan facility can deliver an energy of 500 J in a pulse 
lasting 500 fs. It is a Nd:YAG laser system amplified 
via CPA and can provide intensities up to 10^"'^ W/cm^. 
Whereas, the Astra Gemini laser consists of two inde- 
pendent Ti:Sa laser beams of 0.5 PW each (energy of 
15 J and a pulse duration of 30 fs), with a maximum 
focused intensity of 10^^ W/cm^. The particular layout 
of the Astra Gemini laser renders this system especially 
versatile for unique applications in strong-field physics, 
where two ultrastrong beams are required. Two laser 
systems are likely to be updated to the Petawatt level 
in Germany. The first one is the Petawatt High-Energy 
Laser for heavy Ion experiments (PHELIX) Nd:YAG 
laser at the Gesellschaft fiir Schwerioncnforschung (GSI) 
in Darmstadt, capable now of delivering an energy of 
120 J in about 500 fs (PHELIX, 2011). The 1-PW 
threshold should be reached by increasing the pulse en- 
ergy to 500 J. Combined with the highly-charged ion 
beams at GSI, the PHELIX facility can be attractive for 
experimental investigations in strong-field quantum elec- 
trodynamics (QED). The second system to be updated 
to 1 PW is the Petawatt Optical Laser Amplifier for Ra- 
diation Intensive Experiments (POLARIS) laser in Jena 
(Hein et al, 2010). At the moment, a power of about 
100 TW (energy of 10 J for a pulse duration of 100 fs) 
has been reached and the goal of 1 PW power should 
be achieved by compressing 120 J in about 120 fs. The 
Scottish Centre for the Application of Plasma-based Ac- 
celerators (SCAPA) research center is one of the main 
initiatives within the Scottish Universities Physics Al- 
liance (SUPA) project dedicated to the high-power laser 
interaction with plasmas. A laser system will be devel- 
oped, which will generate pulses of 5-7 J energy and of 



4 



25-30 fs duration at a repetition rate of 5 Hz, correspond- 
ing to a peak power of 200-250 TW, with potential for 
future upgrades to the petawatt level (SCAPA, 2012). 

The 1-PW threshold has been also exceeded in Ti:Sa 
laser systems like those described in Sung et al., 2010 
(energy of 34 J for a pulse duration of 30 fs) and in 
Wang et al., 2011 (energy of 32.3 J for a pulse duration 
of 27.9 fs) and constructed at the Advanced Photonics 
Research Institute (APRI) at Gwangju (Republic of Ko- 
rea) and at the Beijing National Laboratory for Con- 
densed Matter Physics in Beijing (China), respectively. 
The BELLA (Berkeley Lab Laser Accelerator) is a Ti:Sa 
laser system under construction at the Lawrence Berke- 
ley National Laboratory (LBNL) at Berkeley (California, 
USA) which will also reach the 1-PW threshold by com- 
pressing 40 J in 40 fs at a repetition rate of 1 Hz (BELLA, 
2012). 

All the above systems require laser energies larger 
than 10 J and this limits the repetition rate of exist- 
ing petawatt lasers in the best situation to about 0.1 Hz 
(Sung et al, 2010). The Ti:Sa Petawatt Field Synthe- 
sizer (PES) system under development in Garching (Ger- 
many) aims to be the first high-repetition rate petawatt 
laser system with an envisaged repetition rate of 10 Hz 
(PES, 2011). By adopting the OPCPA technique the 
PES should reach the petawatt level by compressing an 
energy of about 5 J in 5 fs (Major et al., 2010). Eor 
a recent review on petawatt-class laser systems, see Ko- 
rzhimanov et al., 2011. 

Einally, we also want to mention other high-power 
lasers, mainly devoted to fast ignition and characterized 
by relatively long pulses of the order of 1 ps-1 ns. Among 
others we mention the OMEGA EP system at Rochester 
(New York, USA) (energy of 1 kJ for a pulse duration 
of 1 ps) (OMEGA EP, 2011) and the National Igni- 
tion Eacility (NIE) at the Lawrence Livermore National 
Laboratory (LLNL) at Livermore (California, USA) (en- 
ergy of 2 MJ distributed in 192 beams with a pulse du- 
ration of about 3-10 ns) (NIE, 2011). Another high- 
power laser facility is the PETawatt Aquitaine Laser 
(PETAL) in Le Barp close to Bordeaux (Erance), which 
is a multi-petawatt laser, generating pulses with energy 
up to 3.5 kJ and with a duration of 0.5 to 5 ps (PETAL, 
2011). The PETAL facility is planned to be coupled to 
the Laser MegaJoule (LMJ) under construction in Bor- 
deaux (France). In the LMJ a total energy of 1.8 MJ is 
distributed in a series of 240 laser beamlines, collected 
into eight groups of 30 beams with a pulse duration of 
0.2-25 ns (LMJ, 2011). 



1. Next-generation 10 PW optical laser systems 

The possibility of building a 10 PW laser system is 
under consideration in various laboratories. At the CLE 
in the United Kingdom the 10 PW upgrade of the Vul- 



can laser has already started (Vulcan lOPW, 2011). The 
new laser will provide beams with an energy of 300 J 
in only 30 fs via the OPCPA. A 10 PW laser system is 
in principle capable of unprecedented intensities larger 
than 10^"^ W/cm^ if the beam is focused to about 1 /im. 
The front-end stage of the Vulcan 10 PW is already com- 
pleted and it delivers pulses with about 1 J of energy at a 
central wavelength of 0.9 ^m, with sufficient bandwidth 
to support a pulse with duration less than 30 fs. 

Another 10 PW laser project is the ILE APOLLON 
to be realized at the Institut de Lumiere Extreme (ILE) 
in France (Chambaret et al., 2009). The laser pulses 
are expected to deliver an energy of 150 J in 15 fs at 
the last stage of amplification after the front end (en- 
ergy of 100 mJ in less than 10 fs), with a repetition rate 
of one shot per minute. Laser intensities of the order 
of lO^-* W/cm^ are envisaged at the ILE APOLLON sys- 
tem, entering the so-called ultrarelativistic regime, where 
also ions (rest energy of the order of 1 GeV) become rel- 
ativistic within one laser period of such an intense laser 
field. 

We mention here also the PEtawatt pARametric Laser 
(PEARL-10) project at the Institute of Applied Physics 
of the Russian Academy of Sciences in Nizhny Novgorod 
(Russia), which is an upgrade of the present 0.56-PW 
laser employing the OPCPA technique, to 10 PW (200 J 
of energy compressed in 20 fs) (Korzhimanov et al., 
2011). 



2. Multi- Petawatt and Exawatt optical laser systems 

The Extreme Light Infrastructure (ELI) (ELI, 2011) 
(see Fig. 1), the Exawatt Center for Extreme Light 
Studies (XCELS) (XCELS, 2012) and the High Power 
laser Energy Research (HiPER) facility at the CLE in 
the United Kingdom (HiPER, 2011) are envisaged laser 
systems with a power exceeding the 100 PW level. 

ELI is a large-scale laser facility consisting of four "pil- 
lars" (see Fig. 1): one devoted to nuclear physics, one 
to attosecond physics, one to secondary beams (photon 
beams, ultrarelativistic electron and ion beams) and one 
to high-intensity physics. This last one is of relevance 
here and it is supposed to comprise ten beams each with 
a power of 10-20 PW that, when combined in phase, 
should deliver a single beam of about 100-200 PW at a 
repetition rate of one shot per minute. The relatively 
high repetition rate is obtained by compressing in each 
beam alone 0.3-0.4 kJ of energy in a pulse of 15 fs. We 
mention that one of the aims of the ILE APOLLON sys- 
tem is to provide a prototype of the 10-20 PW beams, 
that will be then employed at ELI. In the high-field pil- 
lar of ELI ultrahigh intensities exceeding 10^^ W/cm^ 
are envisaged, which are well above the ultrarelativistic 
regime. At such intensities, it will be possible to test dif- 
ferent aspects of fundamental physics for the first time. 
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Lasers planned for the Extreme Light Infrastructure* 
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FIG. 1 (Color online) Summary of the four pillars of ELI. A 
power value of 10(x2) PW indicates the availability of two 
laser systems each with 10 PW power. Reprinted with per- 
mission from Feder, 2010. Copyright 2010, American Institute 
of Physics. 



The XCELS infrastructure is planned to be built in 
Nizhny Novgorod (Russia) and it will consist of 12 beams 
each with energy of 300-400 J and with duration of 
20-30 fs. The pulse resulting from the superposition of 
these beams is expected to have a power of 200 PW, a 
pulse duration of about 25 fs, and divergence less than 
3 diffraction limits (at a central wavelength of 0.91 /im). 
Apart from aiming to overcome the 100 PW threshold, 
the main priorities of XCELS are the creation of sources 
of attosecond and subattosecond, X-ray and 7-ray pulses, 
the development of laser based electron and ion acceler- 
ators with electron and ion energies exceeding 100 GeV 
and up to 10 GeV, respectively, the realization in labo- 
ratory of astrophysical and early-cosmological conditions 
and the investigation of the structure of the quantum 
vacuum. 

The main goal of the other large-scale facility HiPER 
is the first demonstration of laser-driven fusion, or fast 
ignition. To this end HiPER will deliver: 1) an energy of 
about 200 kJ distributed in 40 beams with a pulse dura- 
tion of several nanoseconds and a photon energy of 3 eV 
in the compression side; 2) an energy of about 70 kJ dis- 
tributed in 24 beams with a pulse duration of 15 ps and 
a photon energy of 2 eV in the ignition side. Employing 
HiPER for high-intensity physics would imply a feasible 
reconfiguration of the ignition side to deliver 10 kJ in 
only 10 fs via the OPCPA technique. This would render 
HiPER a laser facility with Exawatt (10^* W) power and 
with a potential intensity of 10^^ W/cm^. 

Finally, we briefly mention the GEKKO EXA facility 
conceptually under design in Osaka (Japan) (GEKKO 
EXA, 2011). This facility is expected to deliver pulses 
of 2 kJ energy and of 10 fs duration corresponding to 
200 PW and with an intensity up to 10^^ W/cm^. 




1 10 10^ 10'^ 10' 

Photon energy [eV] 



FIG. 2 (Color online) Comparison among the peak brilliances 
of the three facilities FLASH, LCLS and European XFEL as 
a function of the laser photon energy. An envisaged peak 
brilliance of 5 x 10^'^ photons/ (s mrad^ mm^ 0.1% bandwidth) 
at a photon energy of 12.4 keV for the SACLA facility is 
reported in European XFEL, 2011. See also European XFEL, 
2011. 



B. Brilliant x-ray laser sources 

Strong optical laser systems are sources of coherent 
radiation at wavelengths of the order of 1 /im, corre- 
sponding to photon energies of the order of 1 cV. Con- 
siderable efforts have been devoted in the past few years 
to develop coherent radiation sources at photon energies 
larger than 100 eV. The discovery of the Self-Amplified 
Spontaneous Emission (SASE) regime (Bonifacio et ai, 
1984) has opened the possibility of employing Free Elec- 
tron Lasers (FELs) to generate coherent light at such 
short wavelengths. In a FEL relativistic bunches of elec- 
trons pass through a spatially-periodic magnetic field 
(undulator) and emit high-energy photons. In the SASE 
regime the interaction of the electron bunch with its own 
electromagnetic field "structures" the bunch itself into 
slices (microbunches) each one emitting coherently even 
at wavelengths below 1 nm (FELs at such small wave- 
lengths are dubbed X-ray Free Electron Lasers (XFELs)). 

The Free-Electron Laser in Hamburg (FLASH) facil- 
ity at the Deutsches Elektronen-SYnchrotron (DESY) 
in Hamburg (Germany) (FLASH, 2011) is one of the 
most brilliant operating FEL's. It delivers short pulses 
(duration of about 10-100 fs) of coherent radiation in 
the extreme ultraviolet-soft x-ray regime (fundamental 
wavelength from 60 nm down to 6.5 nm corresponding 
to photon energies from 21 to 190 eV) at a repetition 
rate of 100 kHz. The intense electron beams available at 
FLASH (total charge of 0.5-1 nC at an energy of 1 GeV) 
allow for peak brilliance of the photon beam of about 
^q29_-^q30 photons/(smrad^mm2 0.1% bandwidth) (see 
Fig. 2), exceeding the peak brilliance of conventional 
synchrotron light sources by several orders of magnitude. 
The Linac Coherent Light Source (LCLS) at Stanford 
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(California, USA) uses the electron beams generated by 
the Stanford Linear Accelerator (SLAC) at the National 
Accelerator Laboratory to generate flashes of coherent 
x-ray radiation of unprecedented brilliance (LCLS, 2011) 
(see also Emma et ai, 2010). Since the electron beam en- 
ergy can be varied from 4.5 GeV to 14.4 GeV, accordingly 
the wavelength of LCLS can be tuned from 1.5 nm to 
0.15 nm (corresponding to photon energies from 0.8 keV 
to 8 keV). The peak brilliance of the LCLS is of the or- 
der of 10^^-10^^ photons/ (smrad^ mm^ 0.1% bandwidth) 
(see Fig. 2), the pulse duration is typically of about 
40-80 fs up to 500 fs, which can be decreased to 10 fs in 
the low-charge electron beam mode, and the repetition 
rate is 120 Hz. 

Another XFEL in operation is the SPring-8 Angstrom 
Compact free electron LAser (SACLA) at the RIKEN 
Harima Institute in Japan (SACLA, 2011). The electron 
accelerator, based on a conducting C-band high-gradient 
radiofrequency acceleration system, and the short-period 
undulator allow for a relatively compact facility of around 
700 m in length (compared, for example, with the about 
2 km of the LCLS). SACLA employs the 8 GeV electron 
beam of the Super Photon Ring - 8 GeV (SPring-8) ac- 
celerator (SPring8, 2011) and has generated x-ray beams 
with 0.08 nm wavelength (corresponding to a photon en- 
ergy of 15.5 keV) at a repetition rate of 60 Hz (see also 
the caption of Fig. 2). 

The European XFEL is under development at DESY in 
Hamburg (Germany) (European XFEL, 2011). It is ex- 
pected to deliver x-ray pulses with a peak brilliance up to 
about 5 X lO'^'^ photons/ (s mrad^ mm^ 0.1% bandwidth) 
at the unprecedented repetition rate of 27 kHz. The elec- 
tron accelerator provides an electron beam with maximal 
energy of 17.5 GeV able to generate laser pulses with a 
central wavelength of 0.05 nm, which corresponds to a 
photon energy of 24.8 keV and with a pulse duration of 
100 fs. Moreover, the European XFEL will be a versatile 
machine consisting of three photon beamlines: the SASE- 
1 and SASE-2, with linearly polarized photons with en- 
ergy in the range 3.1-24.8 keV and the SASE-3, with 
linearly or circularly polarized photons of energy in the 
range 0.26-3.1 keV. 

We also mention that coherent attosccond pulses of 
XUV radiation (photon energy of the order of 100 
eV) have been generated employing HHG in a gaseous 
medium (Agostini and DiMauro, 2004). This technique 
allows for the production of beams with central photon 
energy up to several keVs (Popmintchev et ai, 2009; San- 
sone et ai, 2006), though with intensities several orders 
smaller than XFELs. Less stable sources of coherent soft 
x-rays are the so-called x-ray lasers, which are based on 
the amplification of spontaneous emission by multiply 
ionized atoms in dense plasmas created by intense laser 
pulses (Suckewer and Jaegle, 2009; Wang et ai, 2008; 
Zeitoun et al, 2004). 



III. FREE ELECTRON DYNAMICS IN A LASER FIELD 



In this Section we review, for the benefit of the reader, 
some important basic results on the dynamics of a free 
electron in a laser field (see also the review Eberly, 1969) 
and link them to recent investigations on the subject. 
Results in the realm of classical and quantum electro- 
dynamics are considered separately. Radiation-reaction 
and electron self-interaction effects are not included here 
and their discussion is developed in Sec. VI. 



A. Classical dynamics 

The motion of a charged particle in a laser field is usu- 
ally associated to an oscillation along the laser polariza- 
tion direction. This is pertinent to the non-relativistic 
regime, while the charge dynamics in the relativistic do- 
main is enriched by new features like the drift along the 
laser propagation direction and other non-dipole effects 
(like the well-known figure-8 trajectory), as well as by 
the sharpening of the trajectory at those instants where 
the velocity along the polarization direction reverses. As 
a consequence, laser-driven relativistic free electrons also 
emit high harmonics of the laser frequency (see Sec. V). 

The classical motion of an electron (electric charge 
e < and mass m) in an arbitrary external electromag- 
netic field F^"(x) is determined by the Lorentz equation 
mdu^^/ds = eF'^'ui/, where u'^ = dx^/ds is the elec- 
tron four-velocity and s its proper time (Landau and 
Lifshitz, 1975). If the external field is a plane wave, 
the field tensor F^'^{x) depends only on the dimensional 
phase (p = (nQx), where Uq = (1, no), with no being the 
unit vector along the propagation direction of the wave. 
In this case, for an arbitrary four- vector v^^ = {v'^,v) 
it is convenient to introduce the notation U|| — tiq ■ v, 
vj_ = V — W||no a-nd w_ = {tiqv) — v*^ — The four- 
vector potential of the wave can be chosen in the Lorentz 
gauge as A''(0) = (0, A(<?!))), with A_((/)) = -A||(0) = 0. 
We indicate as = (e,p) = mw^ the (kinetic) four- 
momentum of the electron. Since a plane-wave field de- 
pends only on 0, the canonical momenta p±{(f>) + eA{(j}) 
and p- {4>) are conserved as they are the conjugated mo- 
menta to the cyclic coordinates x±_ and t + x\\, respec- 
tively. For p^{4>o) = = (£ro,po) = 77170(1,^0) being 
the initial condition for the electron's four-momentum at 
a given phase 4>o^ the above-mentioned conservation laws 
already allow to write the electron's four-momentum at 
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an arbitrary phase </> as (Landau and Lifshitz, 1975) 



==£0 - e- 



- A(0o)]^ 

2 po,- 

PJ-(0) =Po.± - e[A((/)) - A{(t)a)], 

/ ,N Po,± • [-4(0) - A(0o)] 

P||(0) =Po.|| - e 

Po,- 

[A(0) - A(0o)]^ 
2 po,- 



(1) 
(2) 
(3) 



where the on-shell condition e((/)) +p||(0) = [p5_(0) + 
TO^]/po,- was employed. For the paradigmatic case of 
a linearly polarized monochromatic plane wave, it is 
A^'{(j)) = A'^cos{ujo(f)), with A'^ = (0, E'qm/wo), where Eq 
is the electric field amplitude, wq the angular frequency 
and u the polarization direction (perpendicular to no). 

The above analytical solution indicates that even if an 
electron is initially at rest, it becomes relativistic within 
one laser period Tq = Iujq if the parameter 



Co 



lel^o 



771^0 



(4) 



is of the order of or larger than unity. In the rela- 
tivistic regime the magnetic component of the Lorentz 
force, which depends on the electron's velocity, becomes 
comparable to the electric one and the electron's dy- 
namics becomes highly nonlinear in the laser-field am- 
plitude. Thus, the parameter is known as classical 
nonlinearity parameter. An heuristic interpretation of 
the parameter is as the work performed by the laser 
field on the electron in one laser wavelength Aq = To in 
units of the electron mass, which clearly explains why 
relativistic effects become important at ^o ^ 1- Alter- 
natively, Eqs. (l)-(3) indicate that the figure-8 trajec- 
tory has a longitudinal (transverse) extension of the or- 
der of Aqi^o (-^oCo): implying that the electron trajectory 
deviates from the unidirectional oscillating one and be- 
comes nonlinear in the field amplitude at ^ 1- Note 
that numerically $o = e.G^ToilO^" W/cm2]Ao[Aini] = 
7.5 v^/o [1020 W/cm^]/wo[eV], where /o = £^0/47^ is the 
wave's peak intensity, and that is gauge- and Lorentz- 
invariant: the gauge invariance has to be intended with 
respect to gauge transformations which do not alter the 
dependence of the four- vector potential on ^ (see Heinzl 
and Ilderton, 2009 for a thorough analysis of this issue) . 
The solution in Eqs. (l)-(3) also indicates that in the 
ultrarelativistic regime at ^o ^ Ij the electron acquires 
a drift momentum along the propagation direction of the 
laser field which is proportional to in contrast to the 
transverse momentum which is proportional to ^q. In 
the case of an electron initially at rest, for example, the 
momentum p(oo) of the electron after the laser pulse 



has been switched off (A(oo) = 0) has the components 
Pi^(cx)) = eA(0o) andp||(oo) ^ <^ A^{f\)Q)l1m,. 

Realistic laser pulses, as those produced in laborato- 
ries, have a more complicated structure than a plane 
wave, essentially because they are spatially focused on 
the transverse planes and the area of the focusing spot 
changes along the laser's propagation axis. Generally 
speaking, if the radius of the minimal focusing area (spot 
radius) is much larger than the central wavelength of the 
laser pulse, then the pulse can be reasonably approxi- 
mated by a plane wave. A Gaussian beam in the paraxial 
approximation offers a more accurate analytical descrip- 
tion of a realistic laser pulse, which shows a Gaussian pro- 
file in the transverse planes (Salamin and Keitel, 2002). 
The dynamics of an electron in such a field cannot be de- 
rived analytically and a numerical solution of the Lorentz 
equation is required (Salamin tt al, 2002). 



B. Quantum dynamics 

In the realm of relativistic quantum mechanics, i.e., 
when e"'"-e~ pair production is negligible (see also Sec. 
VIII) and the single-particle quantum theory is appli- 
cable, the dynamics of an electron in an external elec- 
tromagnetic field with four- vector potential A''(a;) is de- 
scribed by the Dirac equation 



[tdf, - eA^{x)] - m}^ = 0, 



(5) 



where 7'' are the Dirac matrices and where 4'(x) is the 
four-component electron bi-spinor (Berestetskii et al., 
1982). Analogously to the classical case, if the exter- 
nal field is a plane wave, the Dirac equation can be 
solved exactly. If = (eo,Po) and cro/2 = ±1/2 are 
the electron's four- momentum and spin at — > — 00 and 
if A'' (—00) = 0, the positive-energy {eq > 0) solution 

Volkov, 1935) 



po,aoix) of Eq. (5) reads (Berestetskii et al., 1982; 



1 



2po,- 



-noi(0) 



(6) 



where in general v = ^^v^ for a generic four- vector , 
where Wpo^uo is a- positive-energy free bi-spinor (Berestet- 
skii et al., 1982), V is the quantization volume, and where 



Spo{x) = -{pox) 



i: 



eipoAic^')) e^A^' 



Po, 



2po, 



(7) 

is the classical action of an electron in the plane 
wave (Landau and Lifshitz, 1975). The above electron 
states are known as positive-energy Volkov states. The 
negative-energy states -p„^-a-„{x) can be formally ob- 
tained by the replacements pg ^ ^Pa and ctq —cro in 
Eq. (6) except for the energy in the square root (the re- 
sulting bi-spinor u_pq.__cto is the corresponding negative- 
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FIG. 3 (Color) Free wave packet evolution in a plane wave 
field. The solid gray line indicates the center of mass tra- 
jectory, coinciding essentially with the classical trajectory, 
and the laser pulse travels from left to right. The blue re- 
gions indicate the copropagating self-adaptive numerical grid. 
Time and space coordinates are given in "atomic units" , with 
1 a.u. = 24 as and 1 a.u. — 0.05 nm, respectively. From 
Bauke and Keitel, 2011. 



energy free bi-spinor ^ (Berestetskii et ai, 1982)). Al- 
though it has been shown long ago that positive- and 
negative-energy Volkov states form a complete set of or- 
thogonal states on the hypersurfaces 4> = const (Ritus, 
1985), the corresponding property on the hypersurfaces 
t = const is not straightforward and it has been proved 
only recently (see Ritus, 1985 and Zakowicz, 2005, and 
Boca and Florescu, 2010 for a proof of the orthogonal- 
ity and of the completeness of the Volkov states, respec- 
tively) . 

Since the Volkov states form a basis of the space of 
the solutions of Dirac equation in a plane wave, they can 
be employed to build electron wave packets and study 
their evolution. A pedagogical example of laser-induced 
Dirac dynamics is displayed in Fig. 3 for a plane wave 
with peak intensity of 6.3 x 10^"^ W/cm^ and central 
wavelength of 2 nm. The figure shows the drift of the 
wave packet in the propagation direction of the wave, its 
spreading and its shearing due to non-dipole effects. In 
Fillion-Gourdeau et al., 2012 an alternative method of 
solving the time-dependent Dirac equation in coordinate 
space is presented, which explicitly avoids the fermion 
doubling, i.e., the appearance of unphysical modes when 
the Dirac equation is discretized. 

As in the classical case, we shortly mention here the 
paradigmatic case of a monochromatic, linearly polar- 
ized plane-wave field A^(0) — Aq cos{u!o(p) . In this case 
the action Sp^, (x) can be written in the form Sp^ (x) = 



^ We point out that the discussed Volkov states "Pipo.io-o (^) ^re 
the so-called Volkov in-states, as they transform into free-states 
in the limit t — ^ — oo (Fradkin et al, 1991). Volkov out-states, 
which transform into free-states in the limit i — ^ oo, can be 
derived analogously and differ from the Volkov in-states only by 
an inconsequential constant phase factor (recall that A(oo} = 0). 



-{qox) + "oscillating terms", with (Ritus. 1985) 



(8) 



The four- vector Qq plays the role of an "effective" four- 
momentum of the electron in the laser field and it is indi- 
cated as electron "quasimomentum" . The corresponding 
electron 



"mass" = TO* = m^yl + Q/2 is known as 
electron's dressed mass. The results for the quasimo- 
mentum gp and the dressed mass to* in the case of a 
circularly polarized laser field with the same amplitude 
and frequency is obtained from the above ones with the 
replacement ~^ 2^g. The quasimomentum coincides 
classically with the average momentum of the electron 
in the plane wave. Correspondingly, the mass dressing 
depends only on the classical nonlinearity parameter 
and it is an effect of the quivering motion of the electron 
in the monochromatic wave (see also the recent review 
Ehlotzky et al, 2009). As we will see in Sec. V.A, it 
is important that conservation laws in QED processes 
in the presence of a monochromatic plane-wave field in- 
volve the quasimomentum for the incoming electrons 
rather than the four- momentum pg . The question of the 
electron dressed mass in pulsed laser fields has been in- 
vestigated in Heinzl et al., 2010a and Mackenroth and 
Di Piazza, 2011. 

In the realm of QED the parameter ^g can also be 
heuristically interpreted as the work performed by the 
laser field on the electron in the typical QED length 
Ac = 1/to sa 3.9 x 10"^"'^ cm (Compton wavelength) in 
units of the laser photon energy wg (see Eq. (4)). This 
qualitatively explains why multiphoton effects in a laser 
field become important at ^g > 1, such that the laser field 
has to be taken into account exactly in the calculations 
(Ritus, 1985). In the framework of QED this is achieved 
by working in the so-called Furry picture (Furry, 1951), 
where the e+-e~ field ^(x) is quantized in the presence 
of the plane-wave field. This amounts essentially in em- 
ploying the Volkov (dressed) states and the correspond- 
ing Volkov (dressed) propagators (Ritus, 1985) instead 
of free particle states and free propagators to compute 
the amplitudes of QED processes. In the Furry picture 
the effects of the plane wave are accounted for exactly 
and only the interaction between the e^-e~ field ^(z) 
and the radiation field J^'''{x) = d^'A^ix) - ^''A^'{x) 
is accounted for by means of perturbation theory. The 
complete evolution of the system "e"^-e~ field-|-radiation 
field" is obtained by means of the S'-matrix 



S = T 



cxp 



-le 



(^x^-^^'^A^ 



(9) 



where T is the time-ordering operator and ^{x) = 
(a;)7° . For an initial state containing only a single elec- 
tron with four-momentum pg, the quantitative descrip- 
tion of the interaction between the electron, the laser field 
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and the radiation field involves in particular the gauge- 
and Lorcntz-invariant quantum parameter 



Xo 



|e|v/-(^o,M.Po)' _ Po,- E, 







m F„ 



(10) 



where F„ = /\e\ = 1.3 x 10^^ V/cm = 4.4 x lO^-^ G 
being the critical electromagnetic field of QED (Ritus, 
1985). The definition of F^r indicates that a constant 
and uniform electric field with strength of the order of 
Fct, provides an e^-e~ pair with an energy of the or- 
der of its rest energy 2to in a distance of the order of 
the Compton wavelength Ac, implying the instability 



of the vacuum under e 



pair creation in the pres- 



ence of such a strong field (Heisenberg and Euler, 1936; 
Sauter, 1931; Schwinger, 1951). In Eq. (10) we con- 
sidered the case of a linearly polarized plane wave of 
the form A'^(0) = A'^ip{(f)), with 4'{4>) being an arbi- 
trary function with max |(i';/'(0)/(i(/)| = 1 and we intro- 
duced the tensor amplitude F^" = fcp — fep^Qj with 
fcp = cJo?^o • -f'^'" ultrarelativistic electron initially 
counterpropagating with respect to the plane wave it is 
Xo = 5.9 X 10-2eo[GeV]v//o[1020 W/cm^]. The parame- 
ter xo can be interpreted as the amplitude of the electric 
field of the plane wave in the initial rest-frame of the elec- 
tron in units of the critical field of QED and it controls 
the magnitude of pure quantum effects like the photon 
recoil in multiphoton Compton scattering and spin ef- 
fects. This is why it is known as "nonlinear quantum 
parameter" . 

Since the probability dP^/dVdt per unit volume and 
unit time of a quantum process is a gauge- and Lorentz 
invariant quantity, for those processes in a plane-wave 
field involving an incoming electron, as, e.g., multiphoton 
Compton scattering, it can depend only on the two pa- 
rameters ^0 and Xo (Ritus, 1985). For an electromagnetic 
field F^'^{x) ~ {E{x),B{x)) cither constant or slowly- 
varying, the quantity dPe/dVdt, calculated in the lat- 
ter case in the leading order with respect to the fields' 
derivatives, can in principle also depend on the two field 
invariants 

,^ix) = ^F^'^ix)F,A^) = -^lE^ix) - B^ix)], (11) 
^(x) = ^F^^''{x)F^,{x) = -E{x) ■ B{x) (12) 

which identically vanish for a plane wave. In the sec- 
ond equation Ff^„{x) = e^j^yapF"'^ {x)/2 is the dual 
field of F^^''{x) and e'^'^"^ is the four-dimensional com- 
pletely anti-symmetric tensor with e^i^a _ (since 
'^{x) is actually a pseudo-scalar function, the probabil- 
ity dPe/dVdt can only depend on ^^^(a;)). Note, how- 
ever, that ff |,^(x)|, |^(x )| < min(l,x^(a;))^;^,, with 
x{x) = |e|A/|(F^i/(a;)pg)2|/m^, then the dependence of 
dPe/dVdt on ^{x) and $^(x) can be neglected. In this 
case the probability dPe/dVdt essentially coincides with 
the analogous quantity calculated for a constant crossed 



field F^", with the replacement F^" Pi'^ix) (Ritus, 
1985). For a monochromatic plane wave with angular fre- 
quency Wo this occurs if ^ 1- As will be seen in Sec. 
V.A, this condition corresponds, e.g., to the formation 
time of multiphoton Compton scattering (~ m/|e|i?o) 
being much shorter than the laser period Tq. 

As has been mentioned, the iS-matrix in Eq. (9) 
describes all possible electrodynamical processes among 
electrons, positrons and photons. The above consid- 
erations can be easily adapted for discussing processes 
involving an initial positron. Whereas, the probability 
dP^/dVdt of a quantum process in a plane- wave field in- 
volving an incoming photon, as, e.g., multiphoton e^-e~ 
pair production, depends on the parameters ^o and 



^0 



|e|V-(fo,^.fc^)^ ^ k_E^ 
m? m Fcr 



(13) 



where — (w, fc) is the four-momentum of the incoming 
photon (see Ritus, 1985 and Sees. VII-IX). For a pho- 
ton counterpropagating with respect to the plane wave 
it is = 5.9 X 10-2c^[GeV] V^o[102" W/cm^]. In the 
case of multiphoton e'^-e~ pair production, the param- 
eter can be interpreted as the amplitude of the elec- 
tric field of the plane wave in units of the critical field 
Fcr in the ccnter-of-mass system of the created electron 
and positron (Ritus, 1985). The above remarks on pro- 
cesses occurring in a constant or slowly-varying back- 
ground field F^^[x) and involving an incoming electron, 
also apply to the case of an incoming photon once one 
replaces x{x) with >c{x) = \e\\J\{Ffj^y{x)k''Y\/w?. 



IV. RELATIVISTIC ATOMIC DYNAMICS IN STRONG 
LASER FIELDS 

When super-intense infrared laser pulses, as those de- 
scribed in Sec. II. A, impinge on an atom, the latter is 
immediately partly or fully ionized (Becker et at, 2002; 
Keitel, 2001; Protopapas et al., 1997). The ejected elec- 
trons experience the typical "zig-zag" motion of a free 
electron in both laser polarization and propagation direc- 
tions (see Eqs. (2)-(3) and Fig. 3) and will not, in gen- 
eral, return to the ionic core. With an enhanced binding 
force on the remaining electrons, the ionization dynamics 
becomes increasingly complex and may experience sub- 
tle relativistic and correlation effects. When the bind- 
ing force of the ionic core and that of the applied laser 
field eventually become comparable, the electrons may 
in special cases return to and interact with the parent 
ion (rcscattering (Corkum, 1993; Kuchiev, 1987; Schafer 
et al., 1993)). This interaction leads, for example, to the 
ejection of other electrons, to the absorption of energy in 
a scattering process or to the emission of high-harmonic 
photons in case of recombination. 
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A. Ionization 

Previously, atomic or molecular ionization was stud- 
ied with laser pulses of intensity below 10^^ W/cm^, and 
the relativistic laser-matter interaction was dominated 
by the plasma community. The pioneering experiment 
reported in Moore et ai, 1999 on the ionization behav- 
ior of atoms and ions in interaction with a laser with 
intensity of 3 x lO^'^ W/cm^ has thus attracted consid- 
erable interest. The laser magnetic field component was 
shown to alter the direction of ionization characteristi- 
cally (see Sec. III. A). This is because in the relativis- 
tic regime the ionized electron in a laser field acquires 
a large momentum along the laser propagation direction 
(see Eq. (3)) and photoelectrons are emitted mostly in 
that direction within a characteristic opening angle 6: 
ta.n9 ^ p±{oo)/p\\{oo) ~ 2/^o (see the discussion below 
Eq. (4)). With highly-charged ions becoming more easily 
available in a wide range of charges, e.g., via super-strong 
laser fields or by passing the ion beams through metallic 
foils, relativistic laser-induced ionization has been fur- 
ther studied (Chowdhury et ai, 2001; Dammasch et at, 
2001; DiChiara et al, 2008, 2010; Gubbini et ai, 2005; 
Palaniyappan et at, 2008; Yamakawa et at, 2003, 2004). 

In this situation rescattering is generally suppressed 
and multiple ionization of atoms and ions takes place 
mostly via direct ionization, especially including tunnel- 
ing. On the theoretical side attention was then focused 
on the relativistic generalization (Milosevic et ai, 2002; 
Popov et ai, 1997; Popov, 2004; Popov et al, 2006) of 
the so-called Perelomov-Popov-Terent'ev (PPT) theory 
or Ammosov-Delone-Krainov (ADK) model (Ammosov 
et ai, 1986; Perelomov, 1967), which describes atomic 
ionization in the quasistatic tunneling regime. While the 
common intuitive interpretation of the laser induced tun- 
neling fails in the relativistic regime (Reiss, 2008), a re- 
vised picture has been proposed in Klaiber et al., 2012. 
The Strong Field Approximation (SEA) (Faisal, 1973; 
Keldysh. 1965; Reiss, 1980), which treats in a universal 
way both the multiphoton and the tunneling regimes of 
strong-field ionization, has also been extended to the rel- 
ativistic regime (Reiss, 1990a,b). Both the PPT theory 
and the SEA assume that the direct ionization process oc- 
curs as a single-electron phenomenon and thus neglects 
atomic structure effects. 

When the tunneling process proceeds very fast, multi- 
electron correlation effects can occur due to the so-called 
shake-up processes. Thus, the detachment of one elec- 
tron from the atom or ion via tunneling modifies the 
self-consistent potential sensed by the remaining elec- 
trons and may result, consequently, in the excitation of 
the atomic core (inelastic tunneling). A strong excitation 
may also trigger the simultaneous escape of several elec- 
trons from the bound state through the potential barrier 
(collective tunneling). These effects were known to oc- 
cur also in the nonrelativistic regime (Zon, 1999, 2000) 
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FIG. 4 (Color) The total Rb"'+ (black lines) and Rb"+ (red 
lines) ion populations in a gaseous target as a function of 
the peak intensity of a linearly-polarized laser field with a 
wavelength of 0.8 fim and with a pulse duration of 5 fs. The 
solid lines display two-electron inelastic tunneling, the dashed 
lines one-electron inelastic tunneling and the dashed-dotted 
lines the results via the PPT theory. Adapted from Kornev 
et ai, 2009. 



but were concealed by competing rescattering effects. In 
Kornev et ai, 2009, it was shown that, in the relativistic 
regime, the role of inelastic and collective tunneling can 
significantly increase and the relativistic PPT rate has 
been generalized in this respect. In a linearly polarized 
field, the rate of inelastic collective tunneling of N 
equivalent electrons from the outer shell of an ion is de- 
scribed by the following universal formula (Kornev et al., 
2009): 

^(^) _ mVe M!(/ + 1/2)^ {1 + m,)\ 



-1 11 



(Ar_i) k; 2^^7VJ1^+i0F y (mj!)2(Z-TOj-)! 



/r,p ^ 2{v-l)N-M+l/2 

2„.3iV-l / \ ^-2NEa/3Eo 

\ J 

(14) 

where a = « 1/137 is the fine-structure constant, 
nil, ■ . ■ , ttln are the magnetic quantum numbers of the 



bound electrons, M 



E,=i ^ 



J' 



tum number of the electrons, 



r(0) 



I is the orbital quan- 

SjLi(^p°- ~ ^p°- is the jth ionization potential of 
a parent ion, is the energy of the core excitation, 
Ea — K^Fci is the atomic field, Z\e\ is the charge of 
the residual ion, X is the adimensional overlap integral 
(see Kornev et ai, 2009 for its precise definition) and 
Ck.1 w (2 1 vY I \/^2jnj , with v = Za/n. According to the 
calculations in Kornev et al., 2009, inelastic and collec- 
tive tunneling effects contribute significantly to the rel- 
ativistic ionization dynamics at intensities larger than 
lO^s W/cm^, thus changing the ionization probability by 
more than one order of magnitude (see Fig. 4). 

Spin effects of bound systems in strong laser fields were 
shown to moderately alter the quantum dynamics and 
its associated radiation via spin-orbit coupling in highly- 
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FIG. 5 (Color) (a) Experimental photoelectron spectra for 
argon at 7o = 1.2 x 10^'' W/cm^ and at an angle of 62° from 
the laser propagation direction. Analytical results are shown 
for all photoelectrons (continuous line) and for the L-shell 
(dashed line). The angular distributions are at an electron 
energy of (b) 60 keV, (c) 400 keV, and (d) 770 keV. From 
DiChiara et al, 2008. 



charged ions already at an intensity of ~' 10^^ W/cm^ 
(Hu and Keitel, 1999; Walser et al, 2002). More recently 
a nonperturbative rclativistic SFA theory has been de- 
veloped, describing circular dichroism and spin effects in 
the ionization of helium in an intense circularly polarized 
laser field (Bhattacharyya et al., 2011). Here, two-photon 
ionization has been studied in the nonrelativistic inten- 
sity range lO^'^-lO^^ W/cm^ with a photon energy of 45 
eV, yielding small relative spin-induced corrections of the 
order of 10"'^. 

A series of experiments has been devoted to the mea- 
surement of atomic multi-electron effects in relativisti- 
cally strong laser fields. In DiChiara et al., 2008, the en- 
ergy distribution of the ejected electrons and the angle- 
resolved photoelectron spectra for atomic photoioniza- 
tion of argon at Jo ~ 10^^ W/cm^ have been inves- 
tigated experimentally. Here, isolation of the single- 
atom response in the multicharged environment has been 
achieved by measuring photoelectron yields, energies, 
and angular distributions as functions of the sample den- 
sity. Ionization of the entire valence shell along with sev- 
eral inner-shell electrons was shown at Iq ~ 10^^-10^^ 
W/cm^. A typical spectrum in the case of linear polar- 
ization is displayed in Fig. 5. An extended plateau-like 
structure appears in the spectrum due to the electrons 
originating from the L-shell and the longitudinal com- 
ponent of the focused laser field. A surprising feature 
is observed in the energy-resolved angular distribution. 
In contrast to the nonrelativistic case with increasing 
rescatterings and, thus, angular-distribution widths at 
high energies, here azimuthally isotropic angular distri- 
butions are observed at low energies (~ 60 keV in Fig. 5), 
which become narrower for high-energy photoelectrons. 
The authors attribute the anomalous broad angular dis- 



tribution for low-energy electrons to electron-correlation 
effects. A similar experiment on the energy- and angle- 
resolved photoionization was later reported for xenon at a 
laser intensity of lO^^ W/cm^ (DiChiara et al, 2010). For 
energies below 0.5 MeV, the yield and the angular distri- 
bution were shown not to be described by a one-electron 
strong-field model, but rather involve most likely mul- 
tielectron and high-energy atomic excitation processes. 
A further experiment on relativistic ionization of the 
methane molecule at Iq 10^^-10^^ W/cm^ (Palaniyap- 
pan et al., 2008) indicated that molecular mechanisms of 
ionization play no role, and that C^"*" ions are produced 
at these intensities mostly via the cross-shell rescattering 
atomic ionization mechanism. All these experimental re- 
sults still await an accurate theoretical description. 

On the computational side, various numerical meth- 
ods have been developed to describe the laser-driven rel- 
ativistic quantum dynamics in highly-charged ions. A 
Fast-Fourier- Transform split-operator code was imple- 
mented in Mocken and Keitel, 2008 for solving the Dirac 
equation in 2-1-1 dimensions by employing adaptive grid 
and parallel computing algorithms. Another method has 
been developed in Selst0 et al., 2009 to solve the 3D 
Dirac equation by expanding the angular part of the 
wave-function in spherical harmonics. The latter was 
applied to hydrogenlike ions in intense high-frequency 
laser pulses with emphasis on investigating the role of 
negative-energy states. In Bauke et al., 2011, the clas- 
sical relativistic phase-space averaging method, general- 
ized to arbitrary central potentials, and the enhanced 
time-dependent Dirac and Klein-Gordon numerical treat- 
ments are employed to investigate the relativistic ioniza- 
tion of highly- charged hydrogenlike ions in short intense 
laser pulses. For ionization dynamics beyond the tunnel- 
ing regime, quantum mechanical and classical methods 
give similar results, for laser wavelengths from the near- 
infrared region to the soft x-ray regime. Furthermore 
a useful procedure has been developed, which employs 
the over-the-barrier ionization yields for highly-charged 
ions, to determine the peak laser field strength of short 
ultrastrong pulses in the range lo ^ 10^®-10^® W/cm^ 
(Hetzheim and Keitel, 2009). In addition, in this arti- 
cle the ionization angle of the ejected electrons is inves- 
tigated by the full quantum mechanical solution of the 
Dirac equation and the laser field strength is shown to be 
also linked to the electron emission angle. The magnetic 
field-induced tilt in the lobes of the angular distributions 
of photoelectrons in laser-induced relativistic ionization 
has also been discussed in Klaiber et al, 2007. 

There are also several new theoretical results for the 
ionic quantum dynamics in strong high-frequency laser 
fields, in the so-called stabilization regime, where the 
ionization rate decreases or remains constant also with 
increasing laser intensity. An unexpected nondipole ef- 
fect has been reported in F0rre et al., 2006 via numeri- 
cally solving the Schrodinger equation for a hydrogenic 
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FIG. 6 (Color online) Dipole (left) and nondipole (right) 
probability densities of the Kramers-Henneberger wave- 
function in the x-z plane for a i-polarized, 10-cycle sin-like 
pulse propagating in the positive z direction (upward), with 
Eo — 1.5 X 10^^ V/cm and ujq = 54 eV. The snapshots are 
taken at t = 0, t = To/2, and t = 1.8 To from top to bot- 
tom. The length of the horizontal line corresponds to about 
50 as ~ 2.7 nm, with as ~ Ac/a ~ 5.3 x 10~® cm being 
the Bohr radius . Note that the scale is logarithmic with four 
contours per decade. From F0rre et ai, 2006. 



atom beyond the dipole approximation. For this purpose 
the Kramers-Henneberger transformation (Henneberger, 
1968; Kramers, 1956) has been employed, i.e., the trans- 
formation to the instantaneous rest-frame of a classical 
free electron in the laser field, and the terms have 
been neglected in the Hamiltonian (the value of con- 
sidered was approximately 0.14). In Fig. 6, the result- 
ing angular distribution of the ejected electrons in the 
nondipole regime of stabilization displays a third unex- 
pected lobe anti-parallel to the laser propagation direc- 
tion, together with the two expected lobes along the laser 
polarization direction. As a classical explanation, a drift 
along the laser propagation direction was identified for 
the bound electron wave packet in the nondipole case (see 
the middle panel of Fig. 6). Inside the laser field the elec- 
tron has a velocity component along the positive z axis 
but this velocity tends to zero at the end of the pulse. 
Thus, the electromagnetic forces alone do not change the 
electron momentum along the propagation direction at 
the end of the pulse. The net effect of the Coulomb forces 
on the electron wave packet is consequently a momentum 
component along the negative z axis: the electron, which 
is most probably situated in the upper hemisphere over 
the pulse, undergoes a momentum kick in the negative z 
direction each time it passes close to the nucleus. A sim- 
ilar effect has been reported for molecules (F0rre et al, 
2007). 

Radiative recombination, being the time-reversed pro- 
cess of photoionization, of a relativistic electron with a 
highly-charged ion in the presence of a very intense laser 
field has been considered in Miiller et ai, 2009. It was 
shown that the strong coupling of the electron to the laser 



field may lead to a very broad energy spectrum of emit- 
ted recombination photons, with pronounced side wings, 
and to characteristic modifications of the photon angular 
distribution. 

Specific features of nondipole quantum dynamics in 
strong and ultrashort laser pulses have also been inves- 
tigated employing the so-called Magnus approximation 
(Dimitrovski et ai, 2009). The dominant nondipole ef- 
fect is found to be a shift of the entire wave-function to- 
wards the propagation direction, inducing a substantial 
population transfer into states with similar geometry. 

The recent experiment reported in Smeenk et al., 2011 
addresses the question of how the photon momenta are 
shared between the electron and ion during laser-induced 
multiphoton ionization. Theoretically, this problem re- 
quires a nondipole treatment, even in the nonrelativis- 
tic case, to take into account explicitly the laser photon 
momentum. Energy-conservation of £-photon ionization 
here means that lujo ~ Ip + Up + K , where Ip is the ion- 
ization energy of the atom. Up = / AmijjQ is the pon- 
deromotive energy and K is the electron's kinetic energy. 
The experimental results in Smeenk et al., 2011, obtained 
using laser fields with wavelength of 0.8 //m and 1.4 ^m in 
the intensity range of lO-'^^-lO^^ W/cm^ has show that the 
fraction of the momentum, corresponding to the number 
of observed photons needed to overcome the ionization 
energy Ip, is transferred to the created ion rather than 
to the photoelectron. The electron carries only the mo- 
mentum corresponding to the kinetic energy K, while 
the ponderomotive energy and the corresponding por- 
tion of the momentum are transferred back to the laser 
field. This experiment shows that the tunneling concept 
for the ionization dynamics is only an approximation. In 
fact, the quasistatic tunneling provides no mechanism to 
transfer linear momentum to the ion, a conclusion that 
agrees with recent concerns in Reiss, 2008. 



B. Recollisions and high-order harmonic generation 

Tunneling in the nonrelativistic regime is generally fol- 
lowed by recollisions with the parent ion along with var- 
ious subsequent effects (Corkum, 1993; Kuchiev, 1987; 
Schafer et al, 1993). A characteristic feature of strong- 
field processes in the relativistic regime is the suppression 
of recollisions due to the magnetically induced relativistic 
drift of the ionized electron in the laser propagation direc- 
tion (see Sec. III. A). Although relativistic effects become 
significant when the parameter exceeds unity, signa- 
tures of the drift in the laser propagation direction can be 
observed already in the weakly relativistic regime ^ 1- 
The drift will have a significant impact on the electron's 
rescattering probability if, at the instant of recollision, 
the drift distance dy in the laser propagation direction is 
larger than the electron's wave packet size a„pj| in that 
direction (Palaniyappan et al., 2006). The drift distance 
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is approximately given by dy ~ Aq^q/^ (see, e.g., Eq. 
(3)). Instead, the wave packet size a^^pji can be estimated 
from Oujpji ~ v\\At, where v\\ is a typical electron velocity 
along the laser propagation direction and At is the excur- 
sion time of the electron in the continuum. The velocity 
I'll can be related to the tunneling time Ttun via the time- 
energy uncertainty: ~ 1/Ttun- In turn, one can 
estimate the tunneling time Ttun as Ttun ~ hun/vb, where 
^tun ~ Ip/\e\Eo is the tunneling length and vt ~ y^2Ip/m 
the velocity of the bound electron. In the above esti- 
mate, it was assumed that the work carried out by the 
laser field along the tunneling length equals Ip. Thus, at 
the rcscattcring moment At Tq, the wave packet size 

a^jp II is of the order of Aq 



e\Er 



rrr'Ip and the role of 



the drift can be characterized by means of the parameter 



{d\\/ a^jip^w) as estimated by 



so 
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IGwo 
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The condition r > 1 determines the parameter region 
over which the signature of the drift becomes conspicu- 
ous. 

As an alternative view on the relativistic drift, the ion- 
ized electron here misses the ionic core when it is ionized 
with zero momentum. Nevertheless, the recollision will 
occur if the electron is ionized with an appropriate initial 
momentum mQ^j^^ see Eq. (3)), opposite to the 

laser propagation direction. The probability Pi{pd) of 
this process is exponentially damped, though, due to the 
nonzero momentum (see, e.g., Salamin et al, 2006): 



Pi{Pd) ^ cxp 



2 (2m/p)3/2 

3 mlelfi 







Amir, 



(16) 



The drift term in the exponent proportional to will be 
important if ^y2mIpp'^/m\e\EQ > 1, which is equivalent 
to the condition r > 3. At near- infrared wavelengths 
(wo ~ 1 eV) and for the ionization energy Ip = 13.6 eV of 
atomic hydrogen, it becomes relevant at laser intensities 
Jo approximately above 3 x 10^^ W/cm^. Then, HHG 
and other recollision phenomena are suppressed. 

The attainability of relativistic recollisions would, how- 
ever, be very attractive for ultrahigh HHG (Kohler et al., 
2012) as well as for the realization of laser controlled 
high-energy (Hatsagortsyan et al., 2006) and nuclear pro- 
cesses (Chelkowski et al., 2004; Milosevic et al., 2004). 
Various methods for counteracting the relativistic drift 
have been proposed, such as by utilizing highly-charged 
ions (Hu and Keitel, 2001; Keitel and Hu, 2002) which 
move relativistically against the laser propagation di- 
rection (Chirila et al., 2004; Mocken and Keitel, 2004), 
by employing Positronium (Ps) atoms (Henrich et al, 
2004), or through preparing antisymmetric atomic (Fis- 
cher et al., 2007) and molecular (Fischer et al., 2006) or- 
bitals. Here the impact of the drift of the ionized electron 
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FIG. 7 (Color) The HHG setup with two counterpropagating 
APTs. After ionization by the laser pulse 1, the ejected elec- 
tron is driven in the same pulse (light blue) , propagates freely 
after the pulse 1 has left (gray dashed) and is driven back to 
the ion by the laser pulse 2 (dark blue). From Kohler et al., 
2011. 



is reduced by the increase of the laser frequency in the 
system's center of mass, an equally strong drift via two 
constituents with equal mass or via appropriate initial 
momenta from antisymmetric orbitals, respectively. 

On the other hand, the laser field can also be modi- 
fied to suppress the relativistic drift by employing tightly 
focused laser beams (Lin et al., 2006), two counter- 
propagating laser beams with linear polarization (Kei- 
tel et al., 1993; Kylstra et al., 2000; Taranukhin, 2000; 
Taranukhin and Shubin, 2001, 2002) or equal- handed cir- 
cular polarization (Milosevic et al, 2004). In the first 
two cases, the longitudinal component in the tightly fo- 
cused laser beam may counteract the drift, or the Lorentz 
force may be eliminated in a small area near the antin- 
odes of the resulting standing wave, respectively. In the 
third case involving circularly polarized light, the rela- 
tivistic drift is eliminated because the electron velocity 
is oriented in the same direction as the magnetic field. 
This setup is well suited for imaging attosecond dynam- 
ics of nuclear processes but not for HHG because of the 
phase- matching problem (Liu et al., 2009). In the weakly 
relativistic regime the Lorentz force may also be com- 
pensated by a second weak laser beam polarized along 
the direction of propagation of the strong beam (Chirila 
et al., 2002). Furthermore, the relativistic drift can be 
significantly reduced by means of special tailoring of the 
driving laser pulse, which strongly reduces the time when 
the electron's motion is relativistic with respect to a si- 
nusoidal laser pulse (Klaiber et al., 2006, 2007). Two 
consecutive laser pulses (Verschl and Keitel, 2007a) or 
a single laser field assisted by a strong magnetic field 
can also be used to reverse the drift (Verschl and Kei- 
tel, 2007b). In addition two strong Attosecond Pulse 
Trains (APTs) (Hatsagortsyan et at, 2008) or an infrared 
laser pulse assisted by an APT (Klaiber et al., 2008) 
have been employed to enhance relativistic recollisions. 
In fact, due to the presence of the APT the ionization 
can be accomplished by one XUV photon absorption and 
the relatively large energy ujx of the XUV-photon with 
uox = Ip + P^il^m can compensate the subsequent mo- 
mentum drift pd ^ Ti^g /4 in the infrared laser field. 

The main motivation for the realization of relativis- 
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tic recollisions is the extension of HHG towards the hard 
x-ray regime with obvious benefits for time-resolved high- 
resolution imaging. In the past couple of decades, nonrel- 
ativistic atomic HHG (Corkum, 1993; Lewenstein et al., 
1994) has been developed as a reliable source of coher- 
ent XUV radiation and attosecond pulses (Agostini and 
DiMauro, 2004) opening the door for attosecond time- 
resolved spectroscopy (Krausz and Ivanov, 2009). Non- 
relativistic HHG in an atomic gas medium allows already 
to generate coherent x-ray photons up to keV energies 
(Sansone et al., 2006) and to produce XUV pulses shorter 
than 100 as (Goulielmakis et al, 2008). The most favor- 
able conversion efficiency for nonrelativistic keV harmon- 
ics is anticipated with mid-infrared driving laser fields 
(Chen et ai, 2010; Popmintchev et al., 2009). However, 
progress in this field has slowed down, especially because 
of the inhibition, alluded to above, of recollisions due 
to optical driving-field intensity above 3 x 10^® W/cm^. 
This indicates the limit on the cut-off frequency ujc of 
nonrelativistic HHG to uJc ~ 3.17C/p ~ 10 keV. 

Another factor hindering HHG at high intensities is 
the less favorable phase-matching. In strong laser fields, 
outer-shell electrons are rapidly ionized and produce a 
large free electron background causing a phase mismatch 
between the driving laser wave and the emitted x-rays. 
The feasibility of phase-matched relativistic HHG in a 
macroscopic ensemble was first investigated in Kohler 
et al., 2011. Here, the driving fields are two counterprop- 
agating APTs consisting of 100 as pulses with a peak 
intensity of the order of 10^^ W/cm^ (see Fig. 7). The 
electron is driven to the continuum by the laser pulse 1 
in Fig. 7, followed by the usual relativistic drift. There- 
after, the laser pulse 2 overtakes the electron, reverses 
the drift and imposes the rescattering, yielding a much 
higher HHG signal than for a conventional laser field at 
the same cut-off energy. Here phase-matching can be ful- 
filled due to an additional intrinsic phase specific to this 
setup, depending on the time delay between the pulses 
and on the pulse intensity. The latter, being unique for 
this laser setup, mainly affects the electron excursion 
time and varies along the propagation direction. The 
phase- matching is achieved by modifying the laser in- 
tensity along the propagation direction and by balanc- 
ing the phase slip due to dispersion with the indicated 
intrinsic phase. Note, however, that HHG in the rela- 
tivistic regime has been observed experimentally rather 
efficiently in laser plasma interactions (Dromey et al, 
2006). 



V. MULTIPHOTON THOMSON AND COMPTON 
SCATTERING 

In this section we discuss one of the most fundamental 
processes in QED in a strong laser field: the emission of 
radiation by an accelerated electron. After reporting on 



recent theoretical investigations on this process, we dis- 
cuss its possible applications for producing high-energy 
photon beams. 



A. Fundamental considerations 

When an electron is wiggled by an intense laser wave, 
it emits electromagnetic radiation. This process occurs 
with absorption of energy and momentum by the electron 
from the laser field and it is named as multiphoton Thom- 
son scattering or multiphoton Compton scattering, de- 
pending on whether quantum effects, like photon recoil, 
are negligible or not. Multiphoton Thomson and Comp- 
ton scattering in a strong laser field have been studied 
theoretically since a long time (see Salamin and Faisal, 
1998; Sarachik and Schappert, 1970; and Sengupta, 1949 
for multiphoton Thomson scattering and Brown and Kib- 
ble, 1964; Goldman, 1964; and Nikishov and Ritus, 1964a 
for multiphoton Compton scattering). The classical cal- 
culation of the emitted spectrum is based on the analyt- 
ical solution in Eqs. (l)-(3) of the Lorentz equation in 
a plane wave and the substitution of the corresponding 
electron trajectory in the Lienard-Wiechert fields (Jack- 
son, 1975; Landau and Lifshitz, 1975). Whereas, as we 
have discussed in Sec. III.B, the quantum calculation of 
the amplitude of the process is performed in the Furry 
picture of QED. As a result, the total emission probabil- 
ity depends only on the two Lorentz- and gauge-invariant 
parameters (see Eq. (4)) and xo (see Eq. (10)). 

The parameter has already been discussed in Sec. 
III. A. In the contest of multiphoton Compton scattering 
this parameter controls in particular the effective order 
£oS of the emitted harmonics, which, for an ultrarela- 
tivistic electron, can be estimated in the following way. 
In order to effectively emit a frequency uj' , the formation 
length If of the process must not exceed the coherence 
length Icoh, because, otherwise, interference effects would 
hinder the emission. Since an electron with instanta- 
neous velocity (3 and energy e = m'^ = mj ^\ — /S^ ^ m 
mainly emits along the direction of /3, within a cone with 
apex angle -d ^ 1/7^ 1, the formation length If can be 
estimated from If ^ q/"1, with g being the instantaneous 
radius of curvature of the electron trajectory (Jackson, 
1975). On the other hand, ^coh = ■k/uj'{1 - (icosd) ^ 
(Baier et al, 1998; Jackson, 1975). By requir- 
ing that If < /coh, we obtain the following estimate 
for the largest-emitted frequency (cut-off frequency) w^: 
uj'^ ^ j^/g. Now, in the average rest-frame of the elec- 
tron, i.e., in the reference frame where the average elec- 
tron velocity along the propagation direction of the laser 
vanishes (see Eq. (3)), it is 7* ^ ^0 (corresponding to the 
energy e* ^ m^o) and l*f ~ AJ/^o, where the upper in- 
dex * indicates the variable in this frame. Consequently, 
^c* ^ ^0^0 ^^'^ the effective order of the emitted har- 
monics is £cS ^ £,0 (note that £cS is a Lorentz scalar). As 



15 



the order of the emitted harmonics corresponds quantum 
mechanically to the number of laser photons absorbed by 
the electron during the emission process, the parameter 
^0 is also said to determine the "multiphoton" character 
of the process. 

On the other hand, the nonlinear quantum parameter 
Xo (see Eq. (10)) in the contest of multiphoton Compton 
scattering controls the importance of quantum effects as 
the recoil of the emitted photon. In fact, we can esti- 
mate classically the importance of the emitted photon 
recoil from the ratio i^'Je and our considerations above 
exactly indicate that w^/e ~ (^ixf^jra ~ xo- Thus, multi- 
photon Thomson scattering is characterized by the con- 
dition Xo ^ 1, while multiphoton Compton scattering 
by Xo ^ 1- This result can also be obtained in the case 
of a monochromatic laser wave starting from the energy- 
momentum conservation relation 



(17) 



in the case in which i laser photons are absorbed in the 
process (Ritus, 1985). Here and g''^ are the quasimo- 
menta of the initial and final electron (see Eq. (8)) and 
k'^ ~ ijj'ri!^ is the four-momentum of the produced pho- 
ton (n'^ = 0). From this expression it is easy to obtain 
the energy w' of the emitted photon as 



^WoPo, 



(18) 



By reminding that l^s ~ ^'^'^ by estimating the typical 
emission angle of the photon (Mackenroth and Di Piazza, 
2011), it is possible to show that w' ^ xo£o at ^ 1- As 
it has also been throughroughly investigated analytically 
and numerically in Boca and Florescu, 2011 and Seipt 
and Kampfer, 2011a, multiphoton Compton and Thom- 
son spectra coincide in the limit xo ~^ 0, although in 
Seipt and Kampfer, 2011a differences have been observed 
numerically in the detailed structure of the classical and 
quantum spectra also for xo ^ 1- The most important 
difference between classical and quantum spectra is cer- 
tainly the presence of a sharp cut-off in the latter as 
an effect of the photon recoil: the energy of the photon 
emitted in a plane wave is limited by the initial energy 
of the electron^. This does not occur classically, as there 
the frequency of the emitted radiation does not have the 
physical meaning of photon energy. The dependence of 
the energy cut-off on the laser intensity has been recently 



In the case of a plane-wave background field, this limitation 
rather concerns the quantity fc'_ of the emitted photon, as 
k'_ =po,-~p!_ <P0,-- However, for an ultrarelativistic electron 
with po,— ^ ^"€0 and initially counterpropagating with respect 
to the laser field, it is po.- ~ 2eo, k'_ fa 2a;' and p'_ ?s 2e' (Baier 
ei al., 1998). 



recognized as a possible experimental signature of multi- 
photon Compton scattering (Harvey et ai, 2009). 

First calculations on multiphoton Thomson and Comp- 
ton scattering have mainly focused on the easiest case of a 
monochromatic background plane wave, either with cir- 
cular or linear polarization. The main results of these 
investigations, like the dependence of the emitted fre- 
quencies on the laser intensities have been recently re- 
viewed in Ehlotzky et ai, 2009. The complete description 
of the multiphoton Compton scattering process with re- 
spect to the polarization properties of the incoming and 
the outgoing electrons, and of the emitted photon in a 
monochromatic laser wave has been presented in Ivanov 
et ai, 2004. Recently significant attention has been de- 
voted to the investigation of multiphoton Thomson and 
Compton scattering in the presence of short and even 
ultrashort plane-wave pulses (we recall that such pulses 
have still an infinite extension in the directions perpen- 
dicular to the propagation direction). In Boca and Flo- 
rescu, 2009 multiphoton Compton scattering has been 
considered in the presence of a pulsed plane wave. The 
angular-resolved spectra are practically insensitive to the 
precise form of the laser pulse for woTq ^ 20, with tq 
being the pulse duration. The main differences with re- 
spect to the monochromatic case are: 1) a broadening of 
the lines corresponding to the emitted frequencies; 2) the 
appearance of sub-peaks, which are due to interference 
effects in the emission at the beginning and at the end of 
the laser pulse. On the one hand, the continuous nature 
of the emission spectrum in a finite pulse in contrast to 
the discrete one in the monochromatic case has a clear 
mathematical counterpart. In both cases, in fact, the 
total transverse momenta P± with respect to the laser 
propagation direction and the total quantity are con- 
served in the emission process (see Sec. HI). However, in 
the monochromatic case the following additional conser- 
vation law holds (for a linearly polarized plane wave, see 
Eq. (17)) 



Jo+Po.ii + 



-2£ujn = 



-k'. 



2pL 



, (19) 



so that the resulting four-dimensional energy-momentum 
conservation law allows only for the emission of the dis- 
crete frequencies in Eq. (18). On the other hand, the 
appearance of sub-peaks has been in particular investi- 
gated in Heinzl et ai, 2010c, where it has been found 
that the number Ns-p of sub-peaks within the first har- 
monic scales linearly with the pulse duration tq and with 
^g: A's-p — 0.24^gTo[fs]. In this paper the effects of spa- 
tial focusing of the driving laser pulse are also discussed. 
The authors investigate in particular the dependence of 
the deflection angle aout undergone by the electron af- 
ter colliding head-on with a Gaussian focused beam as a 
function of the impact parameter b (see Fig. 8). 

By further decreasing the laser pulse duration, it has 
been argued that effects of the relative phase between the 
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FIG. 8 (Color online) Deflection Oout of an electron initially 
counterpropagating with respect to a laser field with an en- 
ergy of 3 J and a pulse duration of 20 fs, as a function of the 
impact parameter b for different laser waist radii wo- From 
Heinzl et ai, 2010c. 



pulse profile and the carrier wave (the so-called Carrier 
Envelope Phase (CEP)) should become visible in multi- 
photon Thomson and Compton scattering. In Boca and 
Florescu, 2009 the case of ultrashort pulses with woTq > 4 
has been discussed and effects of the CEP on the har- 
monic yield in specific frequency ranges have been ob- 
served. In Mackenroth et ai, 2010 the dependence of the 
angular distribution of the emitted radiation in multi- 
photon Thomson and Compton scattering on the CEP of 
few-cycles pulses has been exploited to propose a scheme 
to measure the CEP of ultrarelativistic laser pulses (in- 
tensities larger than 10^° W/cm^). The method is es- 
sentially based on the high directionality of the photon 
emission by an ultrarelativistic electron, because the tra- 
jectory of the electron, in turn, also depends on the laser's 
CEP. Accuracies in the measurement of the CEP of the 
order of a few degrees are theoretically envisaged. Multi- 
photon Compton scattering in one-cycle laser pulses has 
been considered in Mackenroth and Di Piazza, 2011 and 
a substantial broadening of the emission lines with re- 
spect to the monochromatic case has been observed. The 
high-directionality of radiation emitted via multiphoton 
Thomson scattering has also been employed as a diag- 
nostic tool in Har-Shemesh and Di Piazza, 2012, where 
a new rather precise method has been proposed to mea- 
sure the peak intensity of strong laser fields (intensities 
between 10'^° W/cm^ and 10^'^ W/cm^) from the angular 
aperture of the photon spectrum. 

The study of multiphoton Thomson and Compton 
scattering in short laser pulses has also stimulated the in- 
vestigation of scaling laws for the photon spectral density 
(Boca and Oprea, 2011; Heinzl et ai, 2010c; Seipt and 
Kampfer, 2011a, b). For example, in Heinzl et ai, 2010c a 
scaling law has been found for backscattered radiation in 
the case of head-on laser-electron collisions, which sim- 
plifies the averaging over the electron-beam phase space. 



A more general scaling law has been determined in Seipt 
and Kampfer, 2011b, which relaxes the previous assump- 
tions on head-on collision and on backscattered radiation 
employed in Heinzl et al, 2010c. Moreover, in Seipt and 
Kampfer, 2011a a simple relation is determined between 
the classical and the quantum spectral densities. Finally, 
in Boca and Oprea, 2011 it is found that in the ultra- 
relativistic case 70 ^ 1 the angular distribution of the 
emitted radiation, integrated with respect to the photon 
energy, only depends on the ratio ^0/70 and not on the 
independent values of ^0 and 70 (see also Mackenroth 
et ai, 2010). 

In the above-mentioned publications the spectral prop- 
erties of the emitted radiation in the classical and quan- 
tum regimes have been considered. In Kim et al., 2009 
and Zhang et al., 2008, instead, the temporal properties 
of the emitted radiation in multiphoton Thomson scatter- 
ing have been investigated. In both papers the feasibility 
of generating single attosecond pulses is discussed. 

Photoemission by a single-electron wave packet via 
Thomson scattering in a strong laser field has been dis- 
cussed in Peatross et ai, 2008. It was shown that the 
partial emissions from the individual electron momentum 
components do not interfere when the driving field is a 
plane wave. In other words, the size of the electron wave 
packet, even when it spreads to the scale of the wave- 
length of the driving field, does not affect the Thomson 
emission. 

Finally, we shortly mention that multiphoton effects 
in Thomson and Compton scattering have been mea- 
sured in various laboratories. The second-harmonic ra- 
diation was first observed in the collision of a 1 keV elec- 
tron beam with a Q-switched Nd:YAG laser, although 
the laser intensity was such that ^0 ~ 0-01 (Englert 
and Rinehart, 1983), and then in the interaction of a 
mode-locked Nd:YAG laser (<^o = 2) with plasma elec- 
trons (Chen et ai, 1998). Multiphoton Thomson scat- 
tering of laser radiation in the x-ray domain has been 
reported in Babzien et ai, 2006 (see Pogorelsky et ai, 
2000 for a similar proof-of-principle experiment). Single- 
shot measurements of the angular distribution of the sec- 
ond harmonic (photon energy 6.5 keV) at various laser 
polarizations have been carried out by employing a 60 
MeV electron beam and a subterawatt CO2 laser beam 
with ^0 ~ 0.35. In the prominent SLAC experiment 
(Bula et at, 1996) multiphoton Compton emission was 
detected for the first time. In this experiment an ultra- 
relativistic electron beam with energy of about 46.6 GeV 
collided with a terawatt Nd:glass laser with an intensity 
of 10^^ W/cm^ (^0 « 0.8 and xo ~ 0.3) and four-photon 
Compton scattering has been observed indirectly via a 
nonlinear energy shift in the spectrum of the outcoming 
electrons. 
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B. Thomson- and Compton-based sources of high-energy 
photon beams 

The single-particle theoretical analysis presented above 
indicates that high-energy photons can be emitted via 
multiphoton Thomson and Compton scattering of an ul- 
trarclativistic electron. For example, an electron with 
initial energy Eq ^ m colliding head-on with an optical 
laser field (wq ~ 1 cV) of moderate intensity (^o ^ 1) is 
barely deflected by the laser field {g ^ A070/C0 ^ -^o) 
and potentially emits photons with energies aj[keV] < 
3.8 X 10~^£g[MeV]. This feature has boosted the idea 
of so-called Thomson- and Compton-based sources of 
high-energy photons as a valid alternative to conven- 
tional synchrotron sources, the main advantages of the 
former being the compactness, the wide tunability, the 
shortness of the photon beams in the femtosecond scale 
and the potential for high brightness. Unlike the exper- 
iments on multiphoton Thomson and Compton scatter- 
ing where laser systems with ^0 1 a-re generally em- 
ployed, Thomson- and Compton-based photon sources 
preferably require lasers with ^0 ^ 1, such that multipho- 
ton effects are suppressed and shorter bandwidths of the 
photon beam are achieved. On the other hand, the elec- 
tron beam quality is crucial for Thomson- and Compton- 
based radiation sources. In particular, the brightness of 
the photon beam scales inversely quadratically with the 
electron beam emittancc, and linearly with the electron 
bunch current density. 

Proof-of-principle experiments have demonstrated 
Thomson- and Compton-based photon sources by cross- 
ing a high-energy laser pulse with a picosecond relativis- 
tic electron beam from a conventional linear electron ac- 
celerator (ChoufFani et ai, 2002; Leemans et ai, 1996; 
Pogorelsky et at, 2000; Sakai et ai, 2003; Schoenlein 
et al, 1996; Ting et al, 1995, 1996). We also mention 
the benchmark experiment carried out at LLNL, where 
photons with an energy of 78 keV have been produced 
with a total flux of 1.3 x 10^ photons/shot, by collid- 
ing an electron beam with an energy of 57 MeV with a 
Ti:Sa laser beam with an intensity of about 10^® W/cm^ 
(Co ~ 0.5) (Gibson et al, 2004). 

Another achievement in the development of Thomson- 
and Compton-based photon sources has been the exper- 
imental realization of a compact all-optical setup, where 
the electrons are accelerated by an intense laser. In 
the first experiment with an all-optical setup (Schwoerer 
et al., 2006), x-ray photons in the range of 0.4 keV to 2 
keV have been generated. In this experiment the elec- 
tron beam was produced by a high-intensity Ti:Sa laser 
beam (/q « 2 x 10^^ W/cm^) focused into a pulsed he- 
lium gas jet. The characteristic feature of the all-optical 
setup is that the electron bunches and, consequently, the 
generated x-ray photon beams have an ultrashort dura- 
tion (^ 100 fs) and a linear size of the order of 10 /im. 
Another advantage is that the electrons can be precisely 



synchronized with the driving laser field. In order to 
further improve the all-optical setup, design parameters 
for a proof-of-concept experiment have been analyzed in 
Hartemann et al., 2007. For the calculation of the Comp- 
ton scattering parameters, a 3D Compton scattering code 
has been used, which was extensively tested for Comp- 
ton scattering experiments performed at LLNL (Brown 
et al., 2004; Brown and Hartemann, 2004; Hartemann 
et al, 2004, 2005) (see Sun and Wu, 2011 for an alter- 
native numerical simulation scheme). It is shown that 
x-ray fluxes exceeding 10^^ s"-'^ and a peak brightness 
larger than 10^^ photons/ (smrad^ mm^ 0.1% bandwidth) 
can be achieved at photon energies of about 0.5 MeV. A 
few years later the Compton-based photon source Mo- 
noEnergetic Gamma-ray (MEGa-ray) has been designed 
at LLNL (Gibson et al., 2010). Production of gamma- 
rays ranging from 75 keV to 0.9 MeV has been demon- 
strated with a peak spectral brightness of 1.5 x 10^^ pho- 
tons/(s mrad^ mm^ 0.1% bandwidth) and with a flux of 
1.6 X 10^ photons/shot. An experimental setup for high- 
flux gamma-ray generation has been constructed in the 
Saga Light-Source facility in Tosu (Japan), by collid- 
ing a 1.4 GcV electron beam with a CO2 laser (wave- 
length 10.6 /im) (Kaneyasu et al., 2011). A flux of about 
3.2 X 10'' photons/s gamma photons with energy larger 
than 0.5 MeV has been obtained. 

In the basic setups of Thomson- and Compton-based 
photon sources the electrons experience the intense laser 
held for a time-interval much shorter than that needed 
to cross the whole laser beam, the former being of the or- 
der of the laser's Rayleigh length divided by the speed of 
light. Thus, the quest for a more intense laser pulse at a 
given power in order to increase the photon yield implies 
a tighter focusing and therefore a shorter effective inter- 
action time, which in turn causes a broadening of the 
photon spectrum. In Debus et al, 2010 the Traveling- 
Wave Thomson Scattering (TWTS) setup is proposed, 
which allows the electrons to stay in the focal region 
of the laser beam during the whole crossing time (see 
Fig. 9). This is achieved by employing cylindrical op- 
tics to focus the laser held only along one direction (red 
lines in Fig. 9) and, depending on the angle between 
the initial electron velocity and the laser wave-vector, by 
tilting the laser pulse front. As a result, an interaction 
length ^ 1 cm-1 m can be achieved and, correspondingly 
very large photon fluxes (up to 5 x 10^° photons/shot 
at 20 keV). An alternative way of reaching longer effec- 
tive laser-electron interaction times has been proposed in 
Karagodsky et al, 2010, where a planar Bragg structure 
is employed to guide the laser pulse and realize Thom- 
son/Compton scattering in a waveguide. In this way, 
the yield of x-rays can be enhanced by about two orders 
of magnitude with respect to the conventional free-space 
Gaussian-beam configuration at given electron beam and 
injected laser power in both configurations. However, 
there are two constraints specific to this setup. On the 
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X-ray pulse 




FIG. 9 (Color) Schematic setup of TWTS with the red hnes 
indicating the laser focal lines. In the notation of Debus et al, 
2010 (j) is the angle between the initial electrons' velocity and 
the laser's wave vector and a is the angle between the laser 
pulse front and the laser propagation direction. Adapted from 
Debus et al, 2010. 



one hand, the electron beam has to have a small angular 
spread in order to be injected into the planar Bragg struc- 
ture without causing wall damage. On the other hand, 
the laser field strength has to be such that ^ 8 x 10"''' 
to avoid surface damage. 

Finally, in Hartemann et al, 2008 a setup has been 
proposed to obtain bright GeV gamma-rays via Comp- 
ton scattering of electrons by a thermonuclear plasma. 
In fact, a thermonuclear deuterium-tritium plasma pro- 
duces intense blackbody radiation with a temperature 

20 keV and a photon density ^ 10^^ cm~^ (Tabak 
et ai, 1994). When a thermal photon with energy 
a; ~ 1 keV counterpropagates with respect to a GeV 
electron (70 ^ 10'^), a Doppler-shiftcd high-energy pho- 
ton uj' ^ JqU! ^ 1 GeV can be emitted on axis, i.e., in 
the same direction of the incoming electron. Since cj' 
has to be smaller than the initial electron energy Sq, a 
kincmatical photon pileup is induced in the emitted pho- 
ton spectrum at Eq (Zeldovich and Sunyaev, 1969) (see 
Fig. 10). This results in a quasimonochromatic GeV 
gamma-ray beam with a peak brightness > lO'^" pho- 
tons/(s mrad^ mm^ 0.1% bandwidth), comparable with 
that of the FLASH (see Sec. II.B). 

The unique features of Thomson- and Compton-based 
photon sources render them a powerful experimental de- 
vice. For example, they can be employed for medi- 
cal radio-isotope production and photo-fission, and for 
studying nuclear resonance fiuorescence for in situ iso- 
tope detection (Albert et al., 2010). In this respect, such 
photon sources will represent the main experimental tool 
for nuclear-physics investigations at the Romanian pil- 
lar of ELI (see Fig. 1). Other possible applications, 
at photon energies beyond the MeV threshold, include 
the production of positron beams (Hugenschmidt et al, 
2012; Omori et at, 2003) as well as the investigations of 
high-energy processes occurring in gamma-gamma and 
gamma- lepton collisions (Telnov, 1990). 
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FIG. 10 (Color online) Normalized on-axis brightness for dif- 
ferent values of the rapidity p = cosh~^ 70 at a plasma tem- 
perature of 20 keV. See Hartemann et al., 2008 for the mean- 
ing of the blue circles at p = 6. Adapted from Hartemann 
et ai, 2008. 



VI. RADIATION REACTION 

The issue of "radiation reaction" (RR) is one of the 
oldest and most fundamental problems in electrodynam- 
ics. Classically it corresponds to the determination of 
the equation of motion of a charged particle, an electron 
for dcfinitencss, in a given electromagnetic field F'^'^{x). 
In fact, the Lorcntz equation mdu'^/ds = eF^^'^u^ (see 
Sec. III. A) does not take into account that the electron, 
while being accelerated, emits electromagnetic radiation 
and loses energy and momentum in this way. The first 
attempt of taking into account the reaction of the radia- 
tion emitted by the electron on the motion of the electron 
itself (from here comes the expression "radiation reac- 
tion"), was accomplished by H. A. Lorentz in the nonrel- 
ativistic regime (Lorentz, 1909). Starting from the known 
Larmor formula Vl — (2/3)e^a'^ for the power emitted 
by an electron with instantaneous acceleration a, Lorentz 
argued that this energy-loss corresponds to a "damping" 
force Fr = {2/i)e'^da/dt acting on the electron. The 
expression of the damping force was generalized to the 
relativistic case by M. Abraham in the form (Abraham, 
1905) 

In order to solve the problem of a radiating electron sclf- 
consistently, P. A. M. Dirac suggested in Dirac, 1938 to 
start from the coupled system of Maxwell and Lorcntz 
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equations 



d\FT,f_Lv + dfjFx^iyx + d^Fr^Xfi = 



(21) 



mo- 



ds 



where Ff{x) = F^'^(x) + F^''{x), with i^^''(x) be- 
ing the "self" electromagnetic field generated by the 
electron four-current j^{x) = e/ ds5{x — x{s))u^ and 
where the meaning of the symbol mo for the electron 
mass will be clarified below. In order to write an "ef- 
fective" equation of motion for the electron which in- 
cludes RR, one "removes" the degrees of freedom of 
the electromagnetic field (Teitelboim, 1971). This is 
achieved in Landau and Lifshitz, 1975 at the level of 
the Lagrangian of the system electron-|-electromagnetic 
field and an interesting connection of the RR problem 
with the derivation of the so-called Darwin Lagrangian 
is indicated. By working at the level of the equations 
of motion (21), one first employs the Green's function 
method and formally determines the retarded solution 
-^Trctl^) of (inhomogeneous) Maxwell's equations: 
F^'^^^^lx) = F'^'ix) + F^'^^^^{x) (Teitelboim, 1971). Sub- 
stitution of rct('^) Lorentz equation eliminates 
the electromagnetic field's degrees of freedom, but it is 
not straightforward because F^^^^^^x) has to be calcu- 
lated at the electron's position, where the electron cur- 
rent diverges. This difficulty is circumvented by model- 
ing the electron as a uniformly charged sphere of radius a 
tending to zero. After performing the substitution, and 
by neglecting terms which vanish in the limit a — > 0, one 
obtains the equation (mo + 6m)du'^ /ds = eF'^'^Ui, + F^ 
with dm = (4/3)e^/a being formally diverging. How- 
ever, it is important to note that the only diverging term 
in the limit a — > is proportional to the electron four- 
acceleration. At this point a sort of "classical renor- 
malization principle" is employed, saying that what one 
measures experimentally as the physical electron mass m 
is the overall coefficient of the four-acceleration du'^/ds. 
Therefore, one sets m = mo -I- dm and obtains the so- 
called Lorentz- Abraham-Dirac (LAD) equation: 

du^' 2 /d^u^' du" du^ ,\ , , 

We point out that renormalization in quantum field the- 
ory is based on the fact that the bare quantities, like 
charge and mass, appear in the Lagrangian density of the 
theory, which is not an observable physical quantity. On 
the other hand, the bare electron mass mo, which is for- 
mally negatively diverging for m to be finite, appears here 
in the system of equations (21), which should "directly" 
provide classical physical observables, like the electron 
trajectory. On the other hand, it is also known that the 
LAD equation is plagued with physical inconsistencies 
like, for example, the existence of the so-called "runaway" 



solutions with an exponentially-diverging electron accel- 
eration, even in the absence of an external field (see the 
books Hartemann, 2001 and Rohrlich, 2007 for reviews 
on these issues). 

It was first shown in Landau and Lifshitz, 1975 that in 
the nonrelativistic limit the RR force given by Eq. (20) 
is much smaller than the Lorentz force, if the typical 
wavelength A and the typical field-amplitude F of the 
external electromagnetic field fulfill the two conditions 



A > aAc, 



(23) 



where F„ is the critical electromagnetic field of QED (see 
Sec. III.B). This allows for the reduction of order in the 
LAD equation, i.e., for the substitution of the electron 
acceleration in the RR force via the Lorentz force divided 
by the electron mass. In order to perform the analogous 
reduction of order in the relativistic case, the conditions 
(23) have to be fulfilled in the instantaneous rest-frame 
of the electron (Landau and Lifshitz, 1975). The result 
is the so-called Landau-Lifshitz (LL) equation 



m— — =eF'"'u^ + -e^ 
ds 3 



+ {F'''^u,){F,,xu^)u^ 

m^ 771^ 

(24) 

The LL equation is not affected by the shortcomings of 
the LAD equation: for example, it is evident that if the 
external field vanishes, so does the electron acceleration. 
Most importantly, the conditions (23) in the instanta- 
neous rest-frame of the electron have always to be fulfilled 
in the realm of classical electrodynamics, i.e., if quantum 
effects are neglected. In order for this to be true, in fact, 
the two weaker conditions A » Ac and F <C F^- have 
to be fulfilled in the instantaneous rest-frame of the elec- 
tron: the first guarantees that the electron's wave func- 
tion is well localized and the second ensures that pure 
quantum effects, like photon recoil or spin effects are 
negligible (Baier et ai, 1998; Berestetskii et ai, 1982; 
Ritus, 1985) (see also Sec. III.B). This observation led 
F. Rohrlich to state recently that the LL equation is the 
"physically correct" classical relativistic equation of mo- 
tion of a charged particle (Rohrlich, 2008). Rohrlich's 
statement is also supported by the findings in Spohn, 
2000, where it is shown that the physical solutions of 
the LAD equation, i.e., those which are not runaway- 
like, are on the critical manifold of the LAD equation 
itself and are governed there exactly by the LL equa- 
tion. On the other hand, since the LL equation is derived 
from the LAD equation, one may still doubt on its rigor- 
ous validity, due to the application in the latter equation 
of the "suspicious" classical mass-renormalization proce- 
dure. However, this procedure is avoided in Gralla et al., 
2009 by employing a more sophisticated zero-size limit- 
ing procedure, where also the charge and the mass of the 
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particle are sent to zero but in such a way that their ra- 
tio remains constant. The authors conclude that at the 
leading-order level the LL equation represents the self- 
consistent perturbative equation of motion for a charge 
without electric and magnetic moment. The motion of 
a continuous charge distribution interacting with an ex- 
ternal electromagnetic field is also investigated by a self- 
consistent model and at a more formal level in Burton 
et al, 2007. 

From the original derivation of the LL equation from 
the LAD equation in Landau and Lifshitz, 1975, it is 
expected that the two equations should predict the same 
electron trajectory, possibly with differences smaller than 
the quantum effects. This conclusion has been recently 
confirmed by analytical and numerical investigations in 
Hadad et ai, 2010 for an external plane-wave field with 
linear and circular polarization and in Bulanov et al, 
2011 for different time-dependent external electromag- 
netic field configurations. An effective numerical method 
to calculate the trajectory of an electron via the LL equa- 
tion, which explicitly maintains the relativistic covari- 
ance and the mass-shell condition = 1, has been ad- 
vanced in Harvey et al., 2011a. An alternative numerical 
method for determining the dynamics of an electron in- 
cluding RR effects has been proposed in Mao et al., 2010. 

We should emphasize that the LL equation is not the 
only equation which has been suggested to overcome the 
inconsistencies of the LAD equation. A list of alternative 
equations can be found in the recent review (Hammond, 
2010) (see also Seto et at, 2011). A phenomenological 
equation of motion, including RR and quantum effects 
related to photon recoil, has been suggested in Sokolov 
et al, 2010a, 2009 (see also Sec. VLB). The authors write 
the differential variation of the electron momentum as 
due to two contributions: one arising from the external 
field and one corresponding to the recoil of the emitted 
photon. The resulting equation can be written as the 
system 



m—— = 
dr 
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(25) 



where r is the time in the "momentarily comoving 
Lorentz frame" of the electron where the spatial compo- 
nents of instantaneously vanish, Iqed is the quantum 
radiation intensity (Ritus, 1985) and Xl = (2/3)aa;g^g. 
The expression of Iqed in the case of a plane wave is 
employed, which is valid only for an ultrarclativistic elec- 
tron in the presence of a slowly- varying and undercritical 
otherwise arbitrary external field (see Sec. III.B). 

It has also to be stressed that the original LAD equa- 
tion is still the subject of extensive investigation (the 
first study of the LAD equation in a plane- wave field was 
performed in Hartemann and Kerman, 1996). In Ferris 
and Gratus, 2011, for example, the origin of the Schott 



term in the RR force, i.e., the term proportional to the 
derivative of the electron acceleration (see Eq. (20)), 
is thoroughly investigated and in Kazinski and Shipulya, 
2011 the asymptotics of the physical solutions of the LAD 
equation at large proper times are obtained. Whereas, in 
Noble et al., 2011 a kinetic theory of RR is proposed, 
based on the LAD equation and applicable to study sys- 
tems of many particles including RR (this last aspect 
is also considered in Rohrlich, 2007). Enhancement of 
RR effects due to the coherent emission of radiation by a 
large number of charges is discussed in Smorenburg et al. , 
2010. In this respect, we mention here that, in order to 
investigate strong laser-plasma interactions at intensities 
exceeding 10^"^ W/cm^, RR effects have been also imple- 
mented in Particle- In-Cell (PIC) codes (Tamburini et al., 
2012, 2010; Zhidkov et al, 2002) by modifying the Vlasov 
equation for the electron distribution function according 
to the LL equation. Specifically, in Zhidkov et al, 2002 
it is shown that in the collision of a laser beam with 
intensity Iq ~ 10^^ W/cm^ with an overdense plasma 
slab, about 35% of the absorbed laser energy is con- 
verted into radiation and that the effect of RR amounts 
to about 20%. One-dimensional (Tamburini et al, 2010) 
and three-dimensional (Tamburini et al, 2012) PIC sim- 
ulations have shown that RR effects strongly depend on 
the polarization of the driving field: while for circular po- 
larization they are negligible even at Iq ^ 10^^ W/cm^, 
at those intensities they are important for linear polariza- 
tion. The simulations also show the beneficial effects of 
RR in reducing the energy spread of ion beams generated 
via laser-plasma interactions (see also Sec. XII. A). A dif- 
ferent beneficial effect of RR on ion acceleration has been 
found in Chen et al, 2011 for the case of a transparent 
plasma: RR strongly suppresses the backward motion of 
the electrons, cools them down and increases the num- 
ber of ions to be bunched and accelerated. Finally, the 
system in Eq. (25) has been implemented in a 3D PIC 
code in Sokolov et al, 2009 showing that a laser pulse 
with intensity 10^^ W/cm^ loses about 27% of its energy 
in the collision with a plasma slab. The same system of 
equations has been employed to study the penetration 
of ultra-intense laser beams into a plasma in the hole- 
boring regime (Naumova et al, 2009) and to investigate 
the process of ponderoniotive ion acceleration at ultra- 
high laser intensities in overcritical bulk targets (Schlegel 
et al, 2009). 



A. The classical radiation dominated regime 

As it was already observed in Landau and Lifshitz, 
1975, the fact that the RR force in the LL equation has 
to be much smaller than the Lorentz force in the instan- 
taneous rest-frame of the electron does not exclude that 
some components of the two forces can be of the same 
order of magnitude in the laboratory system. This occurs 
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if the condition a7^F/Fcr '--^ 1 is fulfilled at any instant, 
with 7 being the relativistic Lorentz factor of the elec- 
tron and F the amplitude of the external electromagnetic 
field. For an ultrarelativistic electron this condition can 
be fulfilled also in the realm of classical electrodynam- 
ics (quantum recoil effects are negligible if jF/Fa <C 1) 
and it characterizes the so-called Classical Radiation- 
Dominated Regime (CRDR). The CRDR has been in- 
vestigated in Shen, 1970 for a background constant and 
uniform magnetic field. In Koga et al., 2005 an equivalent 
definition of the CRDR in the presence of a background 
laser field has been formulated, as the regime where the 
average energy radiated by the electron in one laser pe- 
riod is comparable with the initial electron energy. By 
estimating the radiated power Vl from the relativistic 
Larmor formula -Pi = -{2/3)e'^{dUf,/ds){du'' /ds) (Jack- 
son, 1975) with du'^/ds {e/m)F^^^Uy, one obtains that 
Vl ^ olxq£,oSq- Therefore, the conditions of being in the 
CRDR are 

Rc = "Xo^o ~ 1, Xo < 1, (26) 

where, as has been seen in Sec. III.B, the second condi- 
tion ensures in particular that the quantum effects like 
photon recoil are negligible. The same condition Rc ~ 1 
has been obtained in Di Piazza, 2008 by exactly solving 
the LL equation (24) for a general plane- wave background 
field. The analytical solution shows in fact that for an 
ultrarelativistic electron the main effect of RR is due to 
the last term in Eq. (24). As a consequence, while the 
quantity w_ (</)) is constant if the equation of motion is 
that due to Lorentz, it decreases here with respect to (f> 
as U-{(j)) = ?io,-//i((/)), where Uq is the four-velocity at 
an initial 00 and 

„«.l + |££/%.(^)\ (27) 

3 Wo \ dtp J 

where the four-potential of the wave has been assumed 
to have the form A^((/)) = ^oV'(0) (see Sec. III.B, be- 
low Eq. 10). This effect has been recently suggested in 
Harvey et al., 2011b as a possible signature to measure 
RR (see also Lehmann and Spatschek, 2011). The two 
conditions in Eq. (26) are, in principle, compatible for 
sufficiently large values of ^o- For example, for an opti- 
cal (wo = 1 eV) laser field with an average intensity of 
10^4 W/cm^ and for an electron initially counterpropa- 
gating with respect to the laser field with an energy of 
20 MeV, it is xo 0.16 and Rc = 1.3. This example 
shows that, in general, it is not experimentally easy to 
enter the CRDR at least with presently-available laser 
systems. In Di Piazza et al., 2009a a different regime has 
been investigated, which is parametrically less demand- 
ing than the CRDR but in which the effects of RR are 
still large. In this regime the change in the longitudi- 
nal (with respect to the laser field propagation direction) 
momentum of the electron due to RR in one laser period 



is of the order of the electron's longitudinal momentum 
itself in the laser field. As a result, it is found that in 
the ultrarelativistic case and for a few-cycle pulse, if the 
conditions 

Rc > > (28) 

^so 

are fulfilled, then the electron is reflected in the laser 
field only if RR is taken into account (see also Harvey and 
Marklund, 2012 for a recent investigation of the electron's 
dynamics in the two complementary regimes 270 ^ ^0 
including RR effects). This can have measurable effects 
if one exploits the high directionality of the radiation 
emitted by an ultrarelativistic electron (sec Sec. V.A). 
The results in Di Piazza et al., 2009a show in fact that 
the apex angle of the angular distribution of the emitted 
radiation, with and without RR effects included, may 
differ by more than 10° already at an average optical 
laser intensity of 5 x 10^^ W/cm^ (^0 ~ 150) and at an 
initial electron energies of 40 MeV (270 « 156) for which 
Rc ~ 0.08. Small RR effects on photon spectra emitted 
by initially bound electrons had already been predicted 
via numerical integration of the LL equation in Keitel 
et al., 1998 well below the CRDR. 

B. Quantum radiation reaction 

The shortcomings of the classical approaches to the 
problem of RR suggest that it can be fully understood 
only at the quantum level. In the seminal paper Moniz 
and Sharp, 1977 the origin of the classical inconsisten- 
cies, like the existence of runaway solutions of the LAD 
equation, were clarified in the nonrclativistic case. The 
authors first show that such inconsistencies are also ab- 
sent in classical electrodynamics if one considers charge 
distributions with a typical radius larger than the classi- 
cal electron radius ro = aXc ~ 2.8 x 10~^^ cm. Going 
to the nonrclativistic quantum theory and by analyzing 
the Heisenberg equations of motion of the electron in an 
external time-dependent field, the authors conclude that 
the quantum theory of a pointlike particle does not ad- 
mit any runaway solutions, provided that the external 
field varies slowly along a length of the order of Ac (this 
is an obvious assumption in the realm of nonrelativis- 
tic theory, as time-dependent fields with typical wave- 
lengths of the order of Ac would in principle allow for 
e+-e~ pair production, see Sec. VIII). From this point 
of view a classical theory of RR has only physical mean- 
ing as the classical limit {h — )- 0) of the corresponding 
quantum theory and the authors indicate that the result- 
ing equation of motion is the nonrclativistic LL equation 
with the bare mass mo (it is shown that the electrostatic 
self-energy of a point charge vanishes in nonrclativistic 
quantum electrodynamics). On the other hand, if one 
considers the quantum equations of motion of a charge 
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distribution and perforins the classical limit before the 
pointlikc limit, then the classical equations of motion of 
the charge distribution are, of course, recovered and, once 
the point-like limit is then performed, runaway solutions 
appear again. The nonrelativistic form of the LL equa- 
tion has been also recovered from quantum mechanics 
in Krivitskii and Tsytovich, 1991 by including radiative 
corrections to the time-dependent electron momentum 
operator in the Heisenberg representation, and by calcu- 
lating the time-derivative of the average momentum in a 
semiclassical state. 

The situation in the relativistic theory is less straight- 
forward because relativistic quantum electrodynamics, 
i.e., QED, is a field theory fundamentally different from 
classical electrodynamics. The first theory of relativis- 
tic quantum RR goes back to W. Hcitler and his group 
(Heitler, 1984; Jauch and Rohrlich, 1976). However, the 
evaluations of the QED amplitudes in Heitler's theory in- 
volve the solution of complicated integral equations and 
it has given a practical result only in the calculation of 
the total energy emitted by a nonrelativistic quantum 
oscillator, with and without RR. 

At first sight one would say that RR effects are auto- 
matically taken into account in QED, because the elec- 
tromagnetic field is treated BjS 0j collection of photons 
that take away energy and momentum, when they are 
emitted by charged particles. However, photon recoil is 
always proportional to h, making it a purely quantum 
quantity with no classical counterpart. Moreover, if one 
calculates the spectrum of multiphoton Compton scat- 
tering in an external plane- wave field, for example, and 
then performs the classical limit Xo ~^ Oi one obtains 
the corresponding multiphoton Thomson spectrum cal- 
culated via the Lorentz equation and not via the LAD 
or the LL equation (see also Sec. V.A). Finally, it has 
also been seen that in classical electrodynamics the RR 
effects may not be a small perturbation on the Lorentz 
dynamics and they cannot be obtained as the result of a 
single limiting procedure. Otherwise they would always 
appear as a small correction. 

In order to understand what RR is in QED, it is more 
convenient to go back to Eq. (21) and to notice that 
the LAD, namely the LL, equation is equivalent to the 
coupled system of Maxwell and Lorentz equations. If one 
determines the trajectory of the electron via the LL equa- 
tion and then calculates the total electromagnetic field 
FJ^'^{x) via the Lienard-Wiechert four-potential (Landau 
and Lifshitz. 1975), one has solved completely the classi- 
cal problem of the radiating electron in the given electro- 
magnetic field. As has been discussed in Sec. III.B, the 
solution of the analogous problem in strong-field QED 
would correspond to completely determine the iS-matrix 
in Eq. (9), as well as the asymptotic state \t — >■ -|-oo) 
for the given initial state |t — > — oo) = |e~), which rep- 
resents a single electron. The first-order term in the 
perturbativc expansion of the S'-matrix corresponds to 



the process of multiphoton Compton scattering and then, 
classically, to the Lorentz dynamics. Whereas, all high- 
order terms give rise to radiative corrections and to high- 
order coherent and incoherent (cascade) processes, and 
determine what we call "quantum RR" . Here, by high- 
order coherent processes is meant those involving more 
than one basic QED process (photon emission by an 
electron/positron or e~^-e~ photoproduction) but all oc- 
curring in the same formation region. Analogously, in 
higher-order incoherent or cascade processes each basic 
QED process occurs in a different formation region. Now, 
in the case of a background plane wave at ^ 1 and 
Xo ^ Ij the quantum effects are certainly important but 
the radiative corrections and higher-order coherent pro- 
cesses scale with a and can be neglected (Ritus, 1972). 
Also, if xo does not exceed unity then the photons emit- 
ted by the electron are mainly unable, by interacting 
again with the laser field, to create e'^-e~ pairs, as the 
pair production probability is exponentially suppressed 
(see also Sees. VIII and IX). Therefore, it can be con- 
cluded that at ^0 ^ 1 and Xo ^ RR- in QED corre- 
sponds to the overall photon recoil experienced by the 
electron when it emits many photons consecutively and 
incoherently (Di Piazza et at, 2010a). 

A qualitative understanding of the above conclusion 
can be attained by assuming that Xo ^ 1 a-nd by esti- 
mating the average number N-y of photons emitted by 
an electron in one laser period at 3> 1. Since the 
probability of emitting one photon in a formation length 
is of the order of a and since one laser period contains 
about ^0 formation lengths (Ritus, 1985) then ^ a^o- 
Also, the typical energy uj' of a photon emitted by an 
electron is of the order w' ^ Xo^Oi then the average en- 
ergy £ emitted by the electron is £" ~ a^oXo^o ~ Rc^o- 
This estimate is in agreement with the classical result ob- 
tained from the LL equation. In other words, the classical 
limit of RR in this regime corresponds to the emission of 
a higher and higher number of photons all with an en- 
ergy much smaller than the electron energy, in such a 
way that even though the recoil at each emission is al- 
most negligible, the cumulative effect of all photon emis- 
sions may have a finite nonnegligible effect. Note that ut' 
and N-y arc both pure quantum quantities and only their 
product £ has a classical analogue in the limit xo ~^ 0- 
These considerations have allowed for the introduction 
in Di Piazza et al, 2010a of the Quantum Radiation- 
Dominated Regime (QRDR), which is characterized by 
multiple emission of photons already in one laser period. 
This regime is then characterized by the conditions 

Rq = a^o « 1, Xo > 1- (29) 

Quantum photon spectra have been calculated numeri- 
cally in Di Piazza et al., 2010a without RR, i.e., by in- 
cluding only the emission of one photon (four-momentum 
/c''^), and with RR, i.e., by including multiple-photon 
emissions (and by integrating with respect to all the four- 
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FIG. 11 (Color online) Quantum photon spectra as a function 
of Tjj' = k'_/po^- calculated with (solid, black hne) and with- 
out (long dashed, red line) RR and the corresponding classical 
ones with (short dashed, blue line) and without (dotted, ma- 
genta line) RR. The error bars in the quantum spectrum with 
RR stem from numerical uncertainties in multidimensional 
integrations. The numerical parameters in our notation are: 
eo = 1 GeV, = 1.55 eV and h = 10^^ W/cm^ {Rq = 1.1 
and xo = 1-8). Adapted from Di Piazza et ai, 2010a. 

momenta of the emitted photons except one indicated as 
k'^). The results show that in the QRDR the effects of 
RR are essentially three (see Fig. 11): 1) increase of the 
photon yield at low photon energies; 2) decrease of the 
photon yield at high photon energies; 3) shift of the max- 
imum of the photon spectrum towards low photon ener- 
gies. Figure 11 also shows that the classical treatment of 
RR (via the LL equation) artificially overestimates the 
above effects, the reason being that quantum corrections 
decrease the average energy emitted by the electron with 
respect to the classical value (Ritus, 1985). However, at 
Xo ^ 1, i-e., when the recoil of each emitted photon is 
much smaller than the electron energy, then the quantum 
spectra converge into the corresponding classical ones. 
As mentioned in Sec. VI, a semiclassical phenomenolog- 
ical approach to RR in the quantum regime has been 
proposed in Sokolov et al, 2009, 2010b. 

Finally, the quantum modifications induced by the 
electron's self-field onto the Volkov states (see Eq. (6)) 
have been recently investigated in Meuren and Di Pi- 
azza, 2011. It is found that the classical expression of 
the electron quasimomcntum in a linearly polarized 
plane wave (see Eq. (8)) admits a correction depend- 
ing on the quantum parameter xo ^-^d also that self-field 
effects induce a peculiar dynamics of the electron spin. 

VII. VACUUM-POLARIZATION EFFECTS 

QED predicts that photons interact with each other 
also in vacuum (Berestetskii et al., 1982). Effects aris- 
ing from this purely quantum interaction are referred to 
as vacuum polarization effects. This is in contrast to 
classical electrodynamics where the linearity of Maxwell's 
equations in vacuum forbids self-interaction of the elec- 




FIG. 12 Vacuum polarization diagram in an external back- 
ground electromagnetic field. The thick electron lines indicate 
electron propagators calculated in the Furry picture, account- 
ing exactly for the presence of the background field. 

tromagnetic field in the vacuum itself. The possibility of 
photon-photon interaction in vacuum, within the frame- 
work of QED, can be understood qualitatively by observ- 
ing that a photon may locally "materialize" into an e^-e~ 
pair which, in turn, interacts with other photons. For the 
same reason a background electromagnetic field can in- 
fluence photon propagation (sec Fig. 12 and Berestetskii 
et ai, 1982). In the latter case the extension If of the 
region where this transformation occurs, i.e., its forma- 
tion length, depends, in principle, on the structure of 
the background field (Baier and Katkov, 2005). How- 
ever, in some cases it can be estimated qualitatively via 
the Heisenberg uncertainty principle from the typical mo- 
mentum p flowing in the e~^-e^ loop in Fig. 12. We con- 
sider, for example, a constant background electromag- 
netic field (or a slowly-varying one, at leading order in 
the space-time derivatives of the field itself) . In this case, 
if the energy to of the incoming photon (see Fig. 12) is at 
most of the order of m, then the momentum p fiowing in 
the e+-e~ loop is of the order of m and If ^ 1/p ^ Ac- 
If a; 3> m the analysis is more complicated and the for- 
mation length strongly depends on the structure of the 
background field. 

From the theoretical point of view it is convenient to 
distinguish between low-energy vacuum-polarization ef- 
fects if a; ^ TO and high-energy ones if to. 

A. Low-energy vacuum-polarization effects 

The scattering in vacuum of a real photon by another 
real photon is possibly the most fundamental vacuum- 
polarization process (Berestetskii et al., 1982) and it 
has not yet been observed experimentally. The total 
cross section of the process depends only on the Lorentz- 
invariant parameter rj = {kik2)/rn?, with fcf and ki^ be- 
ing the four-momenta of the colliding photons or, equiv- 
alently, on the energy w* of the two colliding photons in 
their center-of-momentum system (77 = 2u;*^/to^). This 
process has been investigated in Euler, 1936 in the low- 
energy limit 77 ^ 1 and then in Akhiezer, 1937 in the 
high-energy limit jy 3> 1 . The complete expression of the 
cross section a^^^-^^ was calculated in Karplus and Neu- 
man, 1950 and can also be found in Berestetskii et al., 
1982 (sec also Fig. 13). In the low-energy limit ?/ ^ 1 
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FIG. 13 Cross section of real photon-photon scattering as a 
function of the energy of the colhding photons in their center- 
of-momentum system in units of the electron mass. Adapted 
from Berestetskii et at, 1982. Copyright Elsevier (1982). 
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and in the high-energy one 77 ^ 1 by (Baier et ai, 1974; 
Berestetskii et ai, 1982) 
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where C(a;) is the Riemann zeta function (Olver et ai, 
2010). In terms of the center-of-momentum energy 



uj* , the cross section becomes cr. 
10-*^6(w*[eV])*^ at w* < m and a 
10-36/(w*[GeV])2 at uj* > m 



= 7.4 X 
= 5.4 X 



The steep dependence 



of (Tjj^jj on 77 for ^ 1 is the main reason why real 
photon-photon scattering has, so-far, eluded experimen- 
tal observation (sec the reviews Marklund and Shukla, 
2006 and Salamin et ai, 2006 for experiments and ex- 
perimental proposals until 2005 aiming to observe real 
photon-photon scattering in vacuum). 

However, various proposals have been put forward re- 
cently in order to observe this process by colliding strong 
laser beams which contain a large number of photons. 
A common theoretical starting point of all these pro- 
posals is the "effcctive-Lagrangian" approach (Dittrich 
and Gies, 2000; Dittrich and Renter, 1985). In this ap- 
proach the interaction among photons in vacuum is de- 
scribed via an effective Lagrangian density of the elec- 
tromagnetic field. By starting from the total Lagrangian 
density of the classical electromagnetic field and of the 
quantum e+-e~ Dirac field, one integrates out the de- 
grees of freedom of the latter field and is left with a 
Lagrangian density depending only on the electromag- 
netic field. As has been seen above, the formation re- 
gion of photon-photon interaction at low energies is of 



the order of Ac, therefore if the classical electromag- 
netic field F^'^{x) = {E{x),B{x)) comprises only wave- 
lengths much larger than Xq, the interaction is approxi- 
mately pointlike and the effective Lagrangian density is 
accordingly a local quantity. Also, since the effective La- 
grangian density is a Lorentz invariant, it can depend 
only on the electromagnetic field invariants ^{x) and 
'^'^{x) already introduced in Sec. III.B. The complete ex- 
pression of the effective Lagrangian density was reported 
for the first time in Heisenberg and Euler, 1936 and Weis- 
skopf, 1936 (see also Schwinger, 1951) and it is known as 
the Euler-Heisenberg Lagrangian density. Here we arc in- 
terested only in the experimentally-relevant low-intensity 
limit \^{x)\,\'i^{x)\ <^ Fj?j. and the leading-order Euler- 
Heisenberg Lagrangian density J^eh{x) reads (Dittrich 
and Gies, 2000; Dittrich and Renter, 1985) 

^i.^(-) = -^^(-) + ^ Ff^ ■ (32) 

Different experimental observables have been suggested 
to detect low-energy vacuum-polarization effects. Those 
will be reviewed next. 



1. Experimental suggestions for direct detection of 
photon-photon scattering 

The most direct way to search for photon-photon scat- 
tering events in vacuum by means of laser fields is to let 
two laser beams collide and to look for scattered photons. 
However, by employing a third "assisting" laser beam, if 
one of the final photons is kinematically allowed to be 
emitted along this beam with the same frequency and po- 
larization, then the number of photon-photon scattering 
events can be coherently enhanced (Varfolomeev. 1966). 
In this laser-assisted setup the "signal" of photon-photon 
scattering is, of course, the remaining outgoing photon. 
In Lundin et ai, 2007 and Lundstrom et ai, 2006 an 
experiment has been suggested to observe laser-assisted 
photon-photon scattering with the Astra-Gemini laser 
system (see Sec. II. A). The authors found a particular 
"three-dimensional" setup, which turns out to be espe- 
cially favorable for the observation of the process (see 
Fig. 14). The number A^^ of photons scattered in one 
shot for an optimal choice of the geometrical factors and 
of the polarization angles between the incoming and the 
assisting beams is found to be 



A^, « 0.25 



Pl[PW]P2[PW]P3[PW] 

{X,[^^m]r 



(33) 



Here, Pi and P2 are the powers of the incoming beams, 
P3 is the power of the assisting beam and A4 is the wave- 
length of the scattered wave to be measured. By plug- 
ging in the feasible values for Astra- Gemini Pi = P2 ~ 
0.1 PW and P3 = 0.5 PW, one obtains A^^ « 0.07, i.e.. 
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FIG. 14 (Color) Schematic "three-dimensional" setup for 
laser-assisted photon-photon scattering involving two incom- 
ing beams (in blue), an assisting one (in red) and a scattered 
one (in violet). From Lundstrom et al, 2006. 



roughly one photon scattered every 15 shots (for the pa- 
rameters of Astra Gemini the wavelength of both incom- 
ing beams is chosen as 0.4 /im, that of the assisting beam 
as 0.8 ^m, so that the wavelength A4 = 0.276 /im of the 
scattered photon is different from those of the incoming 
and assisting beams). 

The quantum interaction among photons in vacuum 
has been exploited in King et at, 2010a to propose, for 
the first time, a double-slit setup comprised only of light 
(see also Marklund, 2010). In this setup two strong par- 
allel beams collide head-on with a counterpropagating 
probe pulse. The photons of the probe have the choice 
to interact either with one or with the other strong beam, 
and, when scattered, they are predicted to build an in- 
terference pattern with alternating minima and maxima 
typical of double-slit experiments (see Fig. 15). Also, 
if one of the slits is closed, i.e., if the probe collides 
with only one strong beam, the interference fringes dis- 
appear. The key idea behind this setup is that the 
vacuum-scattered beam (intensity /^), although propa- 
gating along the same direction as the probe, has a much 
wider angular distribution than the latter, offering the 
possibility of detecting vacuum-scattered photons out- 
side the focus of the undifFracted probe beam. For a 
strong- field intensity of /o ~ 5 x lO^"' W/cm^ that may 
be in the near future be available at ELI or at HiPER 
and for a probe beam with wavelength Xp = 0.527 /im 
and intensity /p = 4 x 10^^ W/cm^, it is predicted that 
about four photons per shot would contribute to build up 
the interference pattern in the observable region (it is the 
region outside the circle in Fig. 15, where Id > 100 /p). 
The diffraction of a probe beam in vacuum by a single fo- 
cused strong laser pulse is also investigated in Tommasini 
and Michinel, 2010 in the case of almost counterpropa- 
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FIG. 15 (Color) Intensity Id of the vacuum-scattered wave 
for a probe beam propagating along the positive y direction 
and colliding with two strong beams aligned along the x axis. 
The crosses correspond to coordinates Xn according to the 
classical prediction x„ = {n+ l/2)\pd/D, where n in an inte- 
ger number, Xp is the wavelength of the probe field, d is the 
distance between the interaction region and the observation 
screen and D is the distance between the centers of the two 
strong beams. The numerical values of the parameters can be 
found in King et al., 2010a. Adapted from King et al., 2010a. 



gating beams. For optimal lasers parameters and at a 
strong-laser power of 100 PW the diffracted vacuum sig- 
nal is predicted to be measurable in a single shot. The 
effects on photon-photon scattering of the temporal pro- 
file of the laser pulses have been recently investigated 
in King and Keitel, 2012, showing a suppression of the 
number of vacuum scattered photons with respect to the 
infinite-pulse (monochromatic) case. 

The concept of Bragg scattering has been exploited in 
Kryuchkyan and Hatsagortsyan, 2011 to observe the scat- 
tering of photons by a modulated electromagnetic-field 
structure in vacuum. If a probe wave passes through 
a series of parallel strong laser pulses and if the Bragg 
condition on the impinging angle is fulfilled, the num- 
ber of diffracted photons can be strongly enhanced. At 
a fixed intensity for each strong beam the enhancement 
factor with respect to laser-assisted photon-photon scat- 
tering is equal to the number of beams in the periodic 
structure. However, in experiments usually the total en- 
ergy of the laser beams is fixed and an enhancement by 
a factor of two is predicted. By considering Nq equal 
Gaussian pulses propagating along the x direction and 
their centers separated by a distance D > 2wz from each 
other, the resulting photon-photon scattering probabil- 
ity will be proportional to the phase- matching factor 
with 



sm^{6k,D/2) ' 



(34) 



where the vector 6k = k2 — fci is the difference between 
the wave vectors of the reflected and incident waves. The 
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Bragg condition is satisfied for Skz = 2Trl/D, with I be- 
ing an integer. By employing ten optical laser beams 
with a wavelength of 1 //m and each with an intensity 
of 2.3 X 10^^ W/cm^, about five vacuum-scattered pho- 
tons are predicted per shot. Finally, an enhancement 
of vacuum-polarization effects in laser-laser collision has 
been predicted in Monden and Kodama, 2011 by em- 
ploying strong laser beams with large angular aperture. 
For example, it is predicted that the number of vacuum- 
radiated photons will be enhanced by two orders of mag- 
nitude, if the angular aperture of the colliding beams is 
increased from 53° to 103°. 

Other experimental suggestions to measure photon- 
photon scattering in vacuum can be found in Eriksson 
et ai, 2004 and Tommasini et ai, 2008. 



2. Polarimetry-based experimental suggestions 

The expression of the Euler-Heiseberg Lagrangian den- 
sity in Eq. (32), suggests to interpret a region where only 
an electromagnetic field is present as a material medium 
characterized by a polarization Peh ~ d^EH / dE — 
E/Att and a magnetization Meh = B^eh / dB + B/Att 
(Jackson, 1975) given by 



Peh = 



1807r2R2 



M 



EH 



1807r2p2, 



[2{E^ - B^)E + 7{E ■ B)B] , (35) 
[2{B^ - E^)B + 7{E ■ B)E] . (36) 



Note that an arbitrary single plane wave cannot "polar- 
ize" the vacuum, as in this case Peh and Meh identi- 
cally vanish. Equations (35) and (36) indicate that the 
presence of an electromagnetic field in the vacuum alters 
the vacuum's refractive index. The situation is even more 
complicated because of the vectorial nature of the back- 
ground electromagnetic field which polarizes the vacuum 
and introduces a privileged direction in it. As a result, 
the vacuum's refractive index is altered in a way that 
depends, in general, on the mutual polarizations of the 
probe electromagnetic field and of the background field: 
the polarized vacuum behaves as a bircfringent medium. 
For example, in the case of an arbitrary constant electro- 
magnetic field {E,B), the refractive indices nEH,i/2 of 
a wave propagating along the direction n and polarized 
along one of the two independent directions Ui = £/\£\ 
and U2 = B/\B\, with S = E - {n ■ E)E + n x B savd 
B = B {n B)B — n X E are given by (Dittrich and 
Gies, 2000) 



nEH.i = 1 



nEH,2 = 1 



4a (n X £;)2 + (n x B)^ - 2n ■ {E x B) 
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(37) 



7a [n X Ef + [n x Bf - 2n ■ {E x B) 
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respectively. 

The birefringence of the polarized vacuum is exploited 
in Heinzl et ai, 2006 to show that if a linearly-polarized 
probe x-ray beam (wavelength Ap) propagates along a 
strong optical standing wave, then it emerges from the 
interaction elliptically polarized with ellipticity e given 
by 



2a Jq,r Iq 

15 A^) Irr 



(39) 



(38) 



where k ~ 1 is a geometrical factor and Iq.r is the 
Rayleigh length of the intense laser beam. If this beam 
is generated by a laser like ELI {Iq ~ 10^^ W/cm^), val- 
ues of the ellipticities of the order of 10~^ are predicted 
at Ap = 0.1 nm. Recent advances on x-ray polarime- 
try allow for measurement of ellipticities of the order of 
10~^ at a wavelength of 0.2 nm (Marx et al, 2011). In 
Ferrando et al., 2007 a phase-shift has been found the- 
oretically to be induced by vacuum-polarization effects 
when two laser beams cross in the vacuum, which is pre- 
dicted to be measurable at laser intensities available at 
ELI or at HiPER. 

When an electromagnetic wave with wavelength A im- 
pinges upon a material body, the features of the scattered 
radiation depend on the so-called diffraction parameter 
& = l\/\d (Jackson, 1975). Here, l± is the spatial di- 
mension of the body perpendicular to the propagation 
direction of the incident wave and d is the distance of 
the screen, where the radiation is detected, from the in- 
teraction region. The near region, ^ 3> 1, is known as 
the "refractive-index limit", because the effects of the 
presence of the body can be described as if the wave 
propagates through a medium with a given refractive 
index. However, if < 1 then the diffraction effects 
become important and description of the wave-body in- 
teraction only in terms of a refractive index is in gen- 
eral not possible. This aspect has been pointed out in 
Di Piazza et al., 2006 within the context of light-light 
interaction in vacuum (see also Di Piazza et al., 2007a). 
Tight focusing required to reach high intensities usually 
renders the interaction region so small that diffraction 
effects may become substantial at typical experimental 
conditions. In some cases diffractive effects reduce by an 
order of magnitude the values of the ellipticity calculated 
via the refractive-index approach and also induce a ro- 
tation of the main axis of the polarization ellipse with 
respect to the initial polarization direction of the probe 
field (Di Piazza et al., 2006). Tight focusing of the strong 
polarizing beam requires quite a detailed mathematical 
description employing a realistic focused Gaussian beam, 
while a simpler description was employed for the usually 
weakly-focused probe beam. This prevented the appli- 
cability of the results in the so-called far region where 
^ <C 1 and where the spatial spreading of the probe field 
is also important. This assumption has been recently re- 
moved in King et al., 2010b, where it was also pointed 
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out that by considering the diffraction of a probe beam 
by two separated beams instead of that by a single stand- 
ing wave, an increase in the elhpticity and in the rotation 
of polarization angle by a factor 1.5 is expected. 

A different method based on the phase-contrast Fourier 
imaging technique has been suggested in Homma et al., 
2011 to detect vacuum birefringence. This technique pro- 
vides a very sensitive tool to measure the absolute phase 
shift of a probe beam when it crosses an intense laser 
field. Numerical simulations demonstrate the feasibility 
of measuring vacuum birefringence also by employing an 
optical probe field and a 100-PW strong laser beam. 

Photon "acceleration" in vacuum due to vacuum polar- 
ization has been studied in Mendonga et al, 2006. This 
effect corresponds to a shift of the photon frequency when 
it passes through a strong electromagnetic wave. If 
is the four-momentum of a probe photon with energy 
LUp when it enters a region where a strong laser beam 
is present, then, due to vacuum-polarization effects, Wp 
becomes sensitive to the gradient of the intensity of the 
strong beam. As a result, a frequency up-shift (down- 
shift) is predicted at the rear (front) of the strong beam. 

According to Eqs. (37) and (38) the phase velocity 
of light in vacuum is smaller than unity. This circum- 
stance has been exploited in Marklund et al., 2005, where 
Cherenkov radiation by ultrarclativistic particles mov- 
ing with constant velocity in a photon gas has been pre- 
dicted, if the speed of the particle exceeds the phase ve- 
locity of light. Finally, in Zimmer et al., 2012 an induced 
electric dipole moment of the neutron has been proposed 
as a signature of the polarization of the QED vacuum. 



3. Low-energy vacuum-polarization effects in a plasma 



where 



no = 



\ 



Z 



(41) 

is the refractive index of the plasma without vacuum- 
polarization effects (Mulser and Bauer, 2010). Equation 
(40) already indicates the possibility of enhancing the ef- 
fects of vacuum polarization by working at laser frequen- 
cies Wo such that no ^ 1, i.e., close to the effective plasma 
critical frequency. This region of parameters is in general 
complex to investigate, due to the arising of different in- 
stabilities. However, the idealized situation investigated 
in Di Piazza et al., 2007b shows, at least in principle, the 
possibility of enhancing the vacuum-polarization effects 
by an order of magnitude at a given intensity Iq with 
respect, for example, to the results in Di Piazza et al., 
2006. The effects of the presence of an additional strong 
constant magnetic field have been analyzed in Lundin 
et al., 2007. The general theory presented in this pa- 
per covers different waves propagating in a plasma as 
Alfven modes, whistler modes, and large-amplitude laser 
modes. We also mention the recent paper Bu and Ji, 
2010, in which the photon acceleration process has been 
investigated in a cold plasma and the reference Brodin 
et al., 2007, where vacuum-induced photon splitting in 
a plasma is studied. Finally, we mention the possibility 
of testing nonlinear vacuum QED effects in waveguides. 
In Brodin et al., 2001 the generation of new modes in 
waveguides due to vacuum-polarization effects has been 
predicted. More recently signatures of nonlinear QED 
effects in the transmitted power along a waveguide have 
been analyzed in Ferraro, 2010. 



Equations (37) and (38) indicate that vacuum- 
polarization effects elicited by a plane wave with intensity 
Jo would alter the vacuum refractive index by an amount 
of the order of (a/457r)(/o//ci )- It was first realized in 
Di Piazza et al., 2007b that this aspect can be in prin- 
ciple significantly improved in a plasma. For the sake 
of simplicity the case of a cold plasma was considered 
and the vacuum-polarization effects were implemented in 
the inhomogeneous Maxwell's equations as an additional 
"vacuum four-current" (Di Piazza et al., 2007b). Now, 
unlike in the vacuum, the field invariant ^(x) for a single 
monochromatic circularly polarized plane wave does not 
vanish in a plasma. Thus, vacuum-polarization effects in 
a plasma already arise in the presence of a single travel- 
ing plane wave. In Di Piazza et al., 2007b the vacuum- 
corrected refractive index n of a two-fluid electron-ion 
plasma (ion mass, density and charge number given by 
rrii. Hi and Z, respectively) in the presence of a circularly 
polarized plane wave has been found as 
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B. High-energy vacuum-polarization effects 

Generally speaking the treatment of vacuum- 
polarization effects for an incoming photon with energy 
uj in the presence of a background electromagnetic field 
with a typical angular frequency cob cannot be performed 
in an effective Lagrangian approach if ujujb/rn^ ^ 1: the 
incoming photon "sees" the nonlocality of its interaction 
with the background electromagnetic field through its 
local "transformation" into an e~^-e~ pair (see Fig. 12). 
The technical difficulty in treating vacuum-polarization 
effects at high energies arises from the fact that the 
interaction between the virtual e+-e^ pair and the 
background field has to be accounted for exactly. This 
has been accomplished for the background field of a 
nucleus with charge number Z such that Za ^ 1 and 
Delbriick scattering and photon splitting in such a field 
have also been observed experimentally (see the reviews 
Lee et al., 2003 and Milstein and Schumacher, 1994). 

As we have recalled in Sec. III.B, the Dirac equa- 
tion in a background plane wave described by the four- 
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vector potential A^{(j)) can be solved exactly and ana- 
lytically. Accordingly, the exact electron (Volkov) prop- 
agator G{x,y\A) in the same background field has also 
been determined (see, e.g., Ritus, 1985). The so-called 
"operator technique", developed in Baier et al, 1976a, b, 
turns out to be very convenient for investigating vacuum- 
polarization effects at high energies (the operator tech- 
nique for a constant background field was developed 
in Baier et al., 1975a, b; and Schwinger, 1951). In 
this technique a generic electron state \E'(a;) in a plane- 
wave field and the propagator G{x,y\A) are intended 
as the configuration representation of an abstract state 
I^E") and of an operator G{A) such that ^{x) = {xl"^) 
and G{x,y\A) = {x\G{A)\y). In particular, since the 
propagator G{x,y\A) is the solution of the equation 
{l''[id^ - eAf,{(j))] - m}G{x,y\A) = 6{x - y), then the 
abstract operator G{A) is simply 

G(A) = -— , (42) 

"/^[Pti - eAfj,{(j))\ - m 

with being the four-momentum operator. Evaluation 
via the operator technique of the matrix element corre- 
sponding to a generic vacuum-polarization process is then 
carried out by manipulating abstract operators, which is 
easier than by working with the corresponding quantities 
in configuration space. 

In Di Piazza et al., 2007c the operator technique has 
been employed to calculate the rate of photon splitting 
in a strong laser field for an incoming photon with four- 
momentum fc^. This was the first investigation of a QED 
process involving three Volkov propagators. The calcu- 
lated rate is valid for an arbitrary plane-wave field, pro- 
vided that radiative corrections can be neglected, i.e., at 
OL>?J^ <g; 1, with xtq = {k^/'m){EQ/ F„), in the most un- 
favorable regime ^o, ^ 1 (Ritus, 1972). In fact, as we 
have mentioned in Sec. III.B, QED processes involving 
the collision of a photon and an intense plane wave are 
controlled by the two Lorentz- and gauge-invariant pa- 
rameters ^0 and >fo- In Di Piazza et al., 2007c it turned 
out to be more convenient to perform a parametric study 
of the photon-splitting rate by varying the two parame- 
ters ^0 and rjo = ujok-fm^ (note that xq = jyoCo)- By 
employing the Furry theorem (Berestetskii et al., 1982), 
it is shown that photon splitting in a laser field only oc- 
curs with absorption of an odd number of laser photons. 
In particular, if the strong field is circularly polarized and 
if it counterpropagates with respect to the incoming pho- 
ton, conservation of the total angular momentum along 
the propagation direction of the beams implies that for 
770 ^ 1 photon splitting can occur only via absorption of 
one or of three laser photons. 

In Di Piazza et al., 2008a a physical scenario has been 
advanced in which nonperturbative vacuum-polarization 
effects can be in principle observed (here "nonperturba- 
tive" means "high-order" in the quantum nonlinearity 
parameter, see below). In this scenario a high-energy 



proton collides head-on with a strong laser field. The 
quantum interaction of the Coulomb field of the proton 
with the laser field allows for a merging of laser photons 
into a single high-energy photon. The use of a proton, 
instead of an electron, for example, is required in order to 
suppress the background process of multiphoton Thom- 
son/Compton scattering, where again many photons of 
the laser can be directly absorbed by the proton and con- 
verted into a single high-energy photon (see Sec. V.A). 
In fact, the kinematics of the two processes, vacuum- 
mediated laser photon merging and multiphoton Thom- 
son/Compton scattering is the same, except that in the 
treatment in Di Piazza et al., 2008a only an even num- 
ber of laser photons can merge in the vacuum-mediated 
process. The probability of £-photon Thomson/Compton 
scattering of a particle with charge Q and mass M de- 
pends on the parameter ^o,c = \Q\Eo/ Mloq and scales 
as ^ofc at ^o,c ^ 1- Thus, the use of a "heavy" particle 
like a proton (mass rup = 1.8 x 10^ to = 938 MeV) is 
essential to suppress this background process. Note that 
in order to have ^o,p = \s\Eo/mpUJo ~ 1 for a proton, 
a laser field with intensity of the order of 10^"* W/cm^ 
is required. It is found that for an optical background 
field such that ^0 ^ 1, the amplitude of the (2£)-photon 
merging process depends only on the nonlinear quantum 
parameter 

(2£) _ Eq 2£{1 + Vp)uJn l~cost? 

^°^P " F„ m 1 + VpCOS^' ^ ' 

where Vp is the proton velocity and ?? is the angle between 
the direction of the emitted photon and the propagation 
direction of the plane wave. By colliding a proton beam, 
of energy available at the Large Hadron Collider (LHC) 
of the order of 7 TeV (LHC, 2011) with an optical laser 
beam of intensity of 3 x 10^^ W/cm^, it is demonstrated 
that the vacuum-mediated merging of two laser photons 
and the analogous two-photon Thomson/Compton scat- 
tering have comparable rates, implying that the inclu- 
sive signal should be twice the one expected without 
vacuum-polarization effects. By accounting for the de- 
tails of the proton beams available at LHC and of the 
laser system PES (sec Sec. II. A), about 670 two-photon 
merging events and about 5 four-photon merging events 
are expected per hour. In the discussed setup most of the 

photons are emitted almost in the direction of the proton 

(2) 

velocity « tt) in which Xo p ^- It is also indicated in 
Di Piazza et al., 2Q08a that the use of the perturbative 
expression of the laser-photon merging rate at leading 
order in Xo.p ^ 1 would lead to an error of about 30%. 
Other setups for observing vacuum-polarization effects in 
laser-proton collisions have been discussed in Di Piazza 
et al., 2008b, involving, for example, XFEL or single in- 
tense XUV pulses. Einally, the process of Delbriick scat- 
tering in a combined Coulomb and laser field has been 
studied in Di Piazza and Milstein. 2008. Here an incom- 
ing photon is scattered by the Coulomb field of a nucleus 
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FIG. 16 Feynman diagrams corresponding to processes (i) 
(part a)), (ii) (part b)) and (iii) (part c)), respectively. The 
thick continuous lines in parts a) and b) indicate Volkov 
positive- and negative-energy states. The crossed vertex in 
part b) stands for the Coulomb electromagnetic field. The 
diagram in part c) is related to the vacuum current j'is.d^), 
that one has to determine in order to calculate the e^-e~ pair 
yield (Dittrich and Gies, 2000). The double line indicates the 
electron propagator calculated in the Furry picture including 
exactly the background standing wave. 

and by a strong laser field. While the presence of the 
laser field is taken into account exactly in the calcula- 
tions, only leading-order effects in the nuclear parameter 
Za are accounted for. Analogously to Di Piazza et al, 
2008a, it is found that high-order nonlinear corrections 
in the parameter >cq to the cross section of the process 
already become important at sa 0.2. For example, 
these corrections amount to about 50% at — 0.35. 



VIII. ELECTRON-POSITRON PAIR PRODUCTION 

One of the most important predictions of QED has 
been the possibility of transforming light into matter 
(Dirac, 1928). If two photons with four-momenta k'^ 
and k2 collide at an angle such that the parameter ij = 
{kik2)/m'^ exceeds two, the creation of an e+-e~ pair be- 
comes kincmatically allowed (Brcit- Wheeler e+-e~ pair 
production (Breit and Wheeler, 1934)). Shortly after the 
realization of the laser in 1960, theoreticians started to 
study possibilities for the creation of e~^-e~ pairs from 
vacuum by very strong laser fields (Nikishov and Ritus, 
1964a; Reiss, 1962; Yakovlev, 1965). Because of con- 
straints from energy-momentum conservation, a single 
plane-wave laser field cannot create pairs from vacuum, 
no matter how intense it is. For a single plane wave, in 
fact, all the photons propagate along the same direction 
and the parameter 77 vanishes identically for any pair of 
photons in the plane wave. Thus, an additional source of 
energy is therefore required to trigger the process of pair 
production in a plane wave. There arc essentially three 
different possibilities (sec also Fig. 16): 

(i) pair production by a high-energy photon propa- 
gating in a strong laser field (multiphoton Brcit- 
Wheeler pair production); 

(ii) pair production by a Coulomb field in the presence 
of a strong laser field; 

(iii) pair production by two counterpropagating strong 
laser beams forming a standing light wave. 



These processes share the common feature to possess 
different interaction regimes which are mainly character- 
ized by the value of the parameter ^q- When ^0 1 
the presence of the laser field can be taken into account 
perturbatively and this yields a pair production rate R 
of the form i?e+-e- ~ "^^o^'"; where £m is the minimum 
integer number which kincmatically allows the process. 
For the process (i) it is £m{kok) > 2m?, with k^ and k^ 
being the four-momentum of the laser photon and the 
incoming photon, respectively. For the process (ii) it is 
^toWq > 2to, where Wg = uiqUc,- is the laser angular fre- 
quency in the rest-frame of the charge which produces the 
Coulomb field and which has four- velocity u^. Finally, 
for the process (iii) it is im<^o > 2m. Due to the specific 
dependence of the pair production rate on ^q, this regime 
of pair production is called multiphoton regime. In con- 
trast, when ^0^1 the presence of the laser field has to 
be taken into account exactly by performing the calcu- 
lations in the Furry picture. As the condition ^0 3> 1 
is realized for vanishing laser frequencies at a fixed laser 
amplitude, this regime is called quasistatic regime and 
the pair-production rate here is governed by a different 
parameter which depends on the process at hand. For 
process (i), for example, the form of the rate depends 
on the physical parameter introduced in Sec. III.B 
and it scales as ^ mx^^^ exp(—8/3>CQ) if >fo ^ 1 a-nd as 
~ m>?J^ if >fo ^ 1 (Nikishov and Ritus, 1964a; Reiss, 
1962). For process (ii), we distinguish the case in which 
the incoming particle is an electron, from that in which 
it is a heavier particle like a nucleus with charge num- 
ber Z . In the first case the pair production rate depends 
on the parameter already introduced in Sec. III.B, 
and the recoil due to the pair creation on the electron 
has to be taken into account (see also Sec. VIII. A). In 
the second case the motion of the nucleus is usually as- 
sumed not to be altered by the pair creation process and 
the nucleus itself is described as a background Coulomb 
field (see also Sec. VIII. B). The pair-production rate 
depends on the parameter xo,n = Un,-{EQ/ Fa-), with 
being the four-velocity of the nucleus, and on the 
nuclear parameter Za. Specifically, the pair-production 
rate scales as to(.Zq;)^ exp(— 2-\/3/xo,n) if Xo,n ^ 1 and 
as TO(ZQ!)^xo,n Inxo.n if Xo,n ^ 1 (Milstein et al, 2006; 
Yakovlev, 1965). Finally, for process (iii), the rate Re+-e~ 
has been derived mainly by approximating the stand- 
ing wave as an oscillating electric field (see also Sec. 
VIII. C). It is found that R^+.e- depends on the ratio 
To = Eq/ Fc-[, with Eq being the amplitude of the stand- 
ing wave in the (fixed) laboratory frame, and that it 
scales as toTq exp(— tt/Tq) if Tq ^ 1 and as toTq if 
To > 1 (Brezin and Itzykson, 1970; Popov, 1971, 1972). 
As expected, these scalings coincide with the correspond- 
ing ones in a constant electric field Eq (Schwinger, 1951). 

The physical meaning of the three parcLineters ^s'o, Xo,n 
and To can be qualitatively understood in the follow- 
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ing way. For process (i) the dressing of the electron 
and positron mass (see Sec. III.B) modifies the thresh- 
old of e+-e~ pair production at ^ 1 according to 
£m{kok) > 2w?^Q. Now, analogously to multiphoton 
Thomson and Compton scattering (see Sec. V.A), the 
typical number of laser photons absorbed in pair pro- 
duction via photon-laser collision is of the order of <^q 
(Nikishov and Ritus, 1964a) and the threshold condition 
becomes xq > 1. Concerning the process (ii), the ap- 
pearance of the parameter xo.n in the quasistatic limit 
^0^1 can be understood by noting that the quantity 
Un.-Ef) is the amplitude of the laser field in the rest- frame 
of the nucleus and that a constant and uniform electric 
field with strength of the order of F^r = m^/|e| supplies 
a e^-e^ pair with its rest energy 2m along the typical 
length scale of QED Ac = 1/m (see also Sec. III.B). 
This last observation also demonstrates the presence of 
the parameter Tq for process (iii). The typical exponen- 
tial scaling of the pair production rate for ^ 1, and 
at ^ 1, Xo,n ^ 1 and Tq ^ 1 for processes (i), (ii) 
and (iii), respectively, is reminiscent of a quantum tun- 
neling process. Thus, one refers to this regime also as 
tunneling pair production (note, however, that the no- 
tion of "tunneling" in laser-induced processes should be 
regarded with special care beyond the dipole approxima- 
tion (Klaiber et aL, 2012; Reiss, 2008)). 

In all processes discussed above, the laser field is al- 
ways participating directly in the pair creation step and 
fundamental properties of the quantum vacuum under 
extreme high-field conditions are probed. However, as 
will be seen shortly, lying at the border of experimental 
feasibility, the expected pair yields are generally rather 
small. It is worth mentioning here that lasers can also 
be applied for abundant generation of e+-e~ pairs (Chen 
et al., 2009, 2010). When a solid target is irradiated by an 
intense laser pulse, a plasma is formed and electrons are 
accelerated to high energies. They may emit radiation 
by bremsstrahlung which efficiently converts into e~^-e~ 
pairs through the Bethe-Heitler process. The laser field 
plays an indirect role in the pair production here by serv- 
ing solely as a particle accelerator. The prolific amount 
of antimatter generated this way may lead to interesting 
applications in various fields of science (Miiller and Kei- 
tel, 2009). Abundant production of e^-e~ pairs and of 
high-energy photons in the collision of a multipetawatt 
laser beam and a solid target has been recently investi- 
gated in Nakamura et aL, 2011 and Ridgers et aL, 2012. 
In particular, in Nakamura et aL, 2011 it has been shown 
that almost all the laser pulse energy is converted after 
the collision into a well collimated high-power gamma- 
ray flash. Whereas, the numerical simulations in Ridgers 
et aL, 2012 indicate that about 35 % of the energy of 
a 10 PW laser pulse after the laser-target interaction is 
converted into a gamma-ray burst and that simultane- 
ously a pure e"^-e~ plasma is produced with a maximum 
positron density of lO'^^ m~'^. 



A. Pair production in photon-laser and electron-laser 
collisions 

Among the pair-production processes mentioned 
above, only laser-induced pair production for < 1 has 
been observed experimentally. Its feasibility has been 
shown, in fact, in the pioneering E-144 experiment at 
SLAC (Bamber et aL, 1999; Burke et aL, 1997) (see 
Reiss, 1971 for a corresponding theoretical proposal). 
The experiment relied on collisions of the 46.6 GeV elec- 
tron beam from SLAC's linear accelerator with a coun- 
terpropagating intense laser pulse of photon energy of 
Wo = 2.4 eV and intensity 1.3 x 10^*^ W/cm^ (^q ~ 0.4). 
In the rest-frame of the electrons, the laser intensity and 
frequency are largely Doppler up-shifted to the required 
level and the pair generation probability is effectively en- 
hanced. In principle both reactions (i) and (ii) contribute 
to pair production in this kind of collisions. Based on 
separate simulations of both production channels, reac- 
tion (i) was found to dominate. The high-energy photon 
originates from multiphoton Compton backscattering of 
a laser photon off the electron beam. 

Despite the significance of the SLAC experiment, a uni- 
fied description of pair creation in electron-laser collisions 
has been presented only recently, treating the competing 
mechanisms (i) and (ii) within the same formalism (Hu 
et aL, 2010). Good agreement with the experimental re- 
sults has been obtained. Moreover, it was shown that the 
SLAC study observed the onset of nonperturbative pair 
creation dynamics, which adds even further significance 
to this benchmarking experiment (sec also Reiss, 2009). 
A formal treatment of the process has also been given in 
Ilderton, 2011, where special emphasis is put on effects 
stemming from the finite duration of the laser pulse. 

Figure 17 shows a survey of various combinations of 
incoming electron energies and optical laser intensities 
which give rise to an observable pair yield. It covers the 
range from the perturbative few-photon regime (^o ~ 0.1 
at Iq « 10^'' W/cm^) to the highly nonperturbative do- 
main (^0 ~ 10 at Iq ~ 10^^ W/cm^), where the con- 
tributions from thousands of photon absorption chan- 
nels need to be included. We note that few-GeV elec- 
tron beams can be produced today using compact laser- 
plasma accelerators (Leemans et aL, 2006) (sec also Sec. 
XII). Future pair creation studies may therefore rely on 
all-optical setups, where a laser-generated electron beam 
collides with a counterpropagating laser pulse. Another 
all-optical setup for pair creation by a seed electron ex- 
posed to two counterpropagating laser pulses was put 
forward in Bell and Kirk, 2008, which will be discussed 
in Sec. IX. 

Pair creation studies could also be conducted as a non- 
standard application of the 17.5 GeV electron beam at 
the upcoming European XFEL beamline at DESY (Euro- 
pean XFEL, 2011), which will normally serve to generate 
coherent x-ray pulses. However, in combination with a 
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FIG. 17 (Color online) Transition from the perturbative to 
the fully nonperturbative regimes of e"^-e~ pair creation in 
electron-laser collisions. The laser photon energy is 2.4 eV. 
Adapted from Hu et al, 2010. 



table-top 10-TW optical laser system, it would also be 
very suitable to probe the various regimes of pair pro- 
duction. In particular, the production channel (ii) could 
be investigated by a suitable choice of beam parameters 
(Hu et al, 2010). 

Other aspects of pair creation by a high-energy photon 
and a strong laser field have been investigated in recent 
years. In Heinzl et al., 2010a, process (i) was considered 
in the case where the laser pulse has finite duration. It 
was found that the finite pulse duration is imprinted on 
the spectra of created particles. Pair production by a 
high-energy photon and an ultrashort laser pulse was also 
considered in Tuchin, 2010. Quantum interference effects 
can arise in photon-induced pair creation in a two-mode 
laser field of commensurate frequencies (Narozhny and 
Fofanov, 2000). 

In addition, the fundamental process (i) may allow for 
applications as a novel tool in ultrashort pulse spectrom- 
etry. A corresponding detection scheme for the char- 
acterization of short gamma-ray pulses of GeV photons 
down to the zeptosecond scale, called Streaking at High 
Energies with Electrons and Positrons (SHEEP), has 
been proposed in Ipp et al., 2011. The basic concept 
of SHEEP is based on e^-e" pair production in vacuum 
by a photon of the test pulse, assisted by an auxiliary 
counter-propagating intense laser pulse. In contrast to 
conventional streak imaging, two particles with opposite 
charges, electron and positron, are created in the same 
relative phase within the third streaking pulse that co- 
propagates with the test pulse. By measuring simulta- 
neously the energy and momentum of the electrons and 
the positrons originating from different positions within 
the test pulse, its length and, in principle, even its shape 
can be reconstructed. The time resolution of SHEEP for 
different classes of tests, streaking and strong pulses can 
range from femtosecond to zeptosecond duration. 



B. Pair production in nucleus-laser collisions 

While in electron-laser collisions the contribution of 
reaction (ii) to pair production is in general small, it be- 
comes accessible to experimental observation when the 
projectile electrons are replaced by heavier particles such 
as protons or other nuclei. The two-step production pro- 
cess via multiphoton Compton scattering will then be 
strongly suppressed by the large projectile mass. The re- 
cent commissioning of the LHC at CERN has stimulated 
substantial activities on pair production in combined 
laser and nuclear Coulomb fields, which may be viewed as 
a generalization of the well-known Bethe-Heitler process 
to strong fields (multiphoton Bethe-Heitler pair produc- 
tion) . The large Lorentz factors 7„ of the ultrarelativistic 
nuclear beams lead to efficient enhancement of the laser 
parameters in the projectile rest-frame. 

Indeed, when a proton beam with Lorentz factor 7^ sa 
3000, as presently available at LHC, collides head-on with 
a superintense laser beam of intensity /q ~ 10^^ W/cm^, 
the Lorentz-boosted laser field strength approaches the 
critical value i^cr- This circumstance motivated the first 
calculations of nonperturbative pair production in colli- 
sions of a relativistic nucleus with a superintense near- 
optical laser beam (Miiller et al., 2003a, b). Smaller pro- 
jectile Lorentz factors may be sufficient, when ultrastrong 
XFEL pulses are employed (Avetissian et al., 2003). The 
calculations were based on an S-matrix treatment and as- 
sumed laser fields of circular polarization. Later on, also 
the case of linear field polarization was studied (Kamihski 
et al, 2006; Krajewska et al, 2006; Miiller et al, 2004; 
Sieczka et al., 2006). This case is rendered more involved 
due to the appearance of generalized Bessel functions, 
which are of very high order when ^0 S> 1. The underly- 
ing S-matrix element is generally of the form 



Xx)Vn{r)^-p^.,^a+{x). 

(44) 

It describes the transition of an electron from the 
negative-energy Volkov state 5'_p^^_o-^ (cc) to a positive- 
energy Volkov state ^_ [x), which is mediated by the 
Coulomb potential Vn{r) = Z\e\/r of the projectile nu- 
cleus. An alternative approach to the problem based on 
the polarization operator in a plane electromagnetic wave 
has been developed in Milstein et al, 2006. It allows 
to obtain total production rates analytically. Both ap- 
proaches rely on the strong-field approximation and in- 
clude the laser field exactly to all orders, whereas the 
nuclear field is treated at leading order in Za. 

Since the high-intensity Bethe-Heitler process has not 
been observed in experiment yet, in recent years physi- 
cists have proposed scenarios which may allow to realize 
the various interaction regimes of the process by present- 
day technology. Few-photon Bethe-Heitler pair produc- 
tion in the perturbative domain could be realized in colli- 
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sions of the LHC proton beam with an XUV pulse of an- 
gular frequency uix ~ 100 eV and of moderate intensity 
Ix ^ 10^^ W/cm2 (Miiller, 2009). Corresponding radi- 
ation sources of table-top dimension are available nowa- 
days in many laboratories. They are based on HHG from 
atomic gas jets or solid surfaces (see Sec. IV. B). The rate 
i?g-i._g_ of pair creation by two-photon absorption close to 
the energetic threshold (i.e., = Un-UJx ^ rn, for the 
angular frequency of the XUV pulse in the rest-frame 
of the nucleus) is given by (Milstein et ai, 2006) 

= 45^^^")'^"^^ V ' ^^^^ 

with J = for linear polarization and j ~ 2 for circular 
polarization. 

In the quasistatic regime of the process sizable pair 
yields require superintense laser fields from a petawatt 
source in conjunction with an LHC proton beam (Miiller 
et al., 2003b; Sieczka et ai, 2006). Such experiments 
will become feasible when petawatt laser pulses are made 
available by high-power devices of table-top size, rather 
than by immobile large-scale facilities as they exist at 
present. A method to enable tunneling pair production 
with more compact multiterawatt laser systems has been 
proposed in Di Piazza et al., 2009b, 2010b. It relies on 
the application of an additional weak XUV field, which 
is superimposed on a powerful optical laser wave. In this 
two-color setup, the energy threshold for pair creation 
can be overcome by the absorption of one photon from 
the high-frequency field and several additional photons 
from the low-frequency field. As a result, by choosing the 
XUV frequency cox such that the parameter S ~ (2m — 
ujx)/m fulfills the conditions < 6 1, the tunneling 
barrier can be substantially lowered and even controlled. 
The pair production rate in the quasistatic regime for 
< 5 <C 1 depends essentially only on the parameters S 
and Xo,m and on the classical nonlinearity parameter 
of the XUV field. For a circularly polarized strong laser 
field it becomes (Di Piazza et ai, 2009b) 

Re^-e- = ^^miZafexxlnVC^exp (-^^) (46) 

for Co — Xo,n/2i5^^^ <C 1, to be compared with the usual 
scaling ~ (— 2-\/3/xo,n) in the absence of the XUV field. 

A related process is laser-assisted Bethe-Heitler pair 
creation, where the high-frequency photon energy sat- 
isfies U!x > 2m. A pronounced channeling of the e~^-e~ 
pair due to the forces exerted by the laser field after their 
creation was found (Lotstedt et al, 2008, 2009). Muhi- 
photon Bethe-Heitler pair creation in a two-color laser 
wave was investigated in Roshchupkin, 2001. 

Analytical formulas for positron energy spectra and 
angular distributions in the tunneling regime of the pro- 
cess were obtained in Kuchiev and Robinson, 2007. For 
pair production at ^ 1 no analytical expressions are 



known because of the intermediate nature of this param- 
eter regime. However, by performing a fitting procedure 
to numerically obtained results, a total pair production 
rate scaling as m(Za)^ exp(— 3.49/xo,n) was obtained in 
Miiller et ai, 2009b, which closely resembles the tunnel- 
ing exponential behaviour m(Za)^ exp(— 2-\/3/xo,n)- 

While in Eq. (44) the influence of the projectile is de- 
scribed by an external Coulomb field, the projectile can 
also be treated quantum particle which allows to 
study nuclear recoil effects (Krajewska and Kamiiiski, 
2010, 2011; Miiller and Miiller, 2009). Besides, in laser- 
nucleus collisions, bound-free pair creation can occur 
where the electron is created in a deeply bound atomic 
state of the nucleus. The process was studied first for 
circular laser polarization (Matveev et ai, 2005; Miiller 
et al., 2003c) and later on also for linear polarization 
(Deneke and Miiller, 2008), including contributions from 
the various atomic subshells. 



C. Pair production in a standing laser wave 

Purely light-induced pair production can occur when 
two noncopropagating laser waves are superimposed. 
The simplest field configuration consists of two coun- 
terpropagating laser pulses of equal frequency and in- 
tensity. The resulting field is a standing wave which is 
inhomogencous both in space and time and a theoreti- 
cal treatment of the process is very challenging. In or- 
der to render the problem tractable and since the pro- 
duction process mainly occurs where the electric field 
component of the background field is stronger than the 
magnetic one (Dittrich and Gies, 2000), in the standard 
approach the resulting standing light wave is approxi- 
mated by a purely electric field oscillating in time. This 
approximation is expected to be justified for a strong 
(Iq > 10^° W/cm^), optical laser field where the typ- 
ical spatial scale of the field variation Aq ~ 1 /^m is 
much larger than the pair formation length m/\e\Eo = 
XcF„/Eo ~ 2.6 X 10-2 /xm/\//o[1020 W/cm^] (Ritus, 
1985). Note the analogy between the formation length 
m/\e\EQ for pair production and the tunneling length 
^tun ^ Ip/\e\Eo in atomic ionization (see Sec. IV. B), 
with Ip being the ionization potential energy. 

Pair production in an oscillating electric field is a 
generalization of the Schwingcr mechanism (Schwinger, 
1951) to time-dependent fields and it has been considered 
by many theoreticians, starting from the seminal works 
Brezin and Itzykson, 1970 and Popov, 1971, 1972 (for a 
comprehensive list of references until 2005, see Salamin 
et ai, 2006). 

While in the laser-electron and laser-nucleus collisions 
of the previous subsections the Doppler boost of the laser 
parameters due to a highly relativistic Lorentz factor 
could be exploited, in laser-laser collisions this is not pos- 
sible so that high field strengths Eq or high frequencies 
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FIG. 18 (Color) Number of pairs produced by the x-ray as- 
sisted Scliwinger mechanism for two different values of the x- 
ray angular frequency bJx , indicated as uj in the figure, (blue 
and black solid lines) and the ratio of these catalyzed pairs to 
those produced by the standard Schwinger mechanism (blue 
and black dashed lines), both as functions of the optical field 
strength in units of Fct- Adapted from Dunne et al, 2009. 

Wo are required in the laboratory frame. Theoreticians 
are therefore aiming to find ways for enhancing the pair 
production probabihty in order to render the process ob- 
servable in the foreseeable future. 

A first possibility to facilitate the observability of pair 
creation in a standing optical laser wave is to superim- 
pose an x-ray photon (or any other high-frequency com- 
ponent) onto the high-field region (Dunne et ai, 2009; 
Monin and Voloshin, 2010; Schiitzhold et ai, 2008). In 
this way, the Schwinger mechanism is catalyzed so that 
the usual exponential suppression ~ exp(— tt/Tq) is sig- 
nificantly lowered. For example, in the limit when the 
x-ray energy approaches the threshold value 2m, the pair 
production rate i?e+-e- becomes 

Re+.e- ~ mexp (^-^^Y^^ (47) 

assuming that the x-ray propagation direction is perpen- 
dicular to the electric field vector of the strong optical 
field. An overview of the pair production enhancement 
effect due to the x-ray assistance is shown in Fig. 18. 

Another proposal to enhance the pair yield is the ap- 
plication of multiple colliding laser pulses instead of only 
two (Bulanov et ai, 2010a). It has been demonstrated 
that the threshold laser energy necessary to produce a 
single pair, decreases when the number of colliding pulses 
is increased. The results are summarized in Table I. Pair 
production exceeds the threshold when eight laser pulses, 
with a total energy of 10 kJ, are simultaneously focused 
on one spot. Doubling (tripling) the number of pulses 
leads to an enhancement by two (six) orders of magni- 
tude. The threshold energy drops from 40 kJ for two 
pulses to 5.1 kJ for 24 pulses, clearly indicating that the 
multiple-pulse geometry is strongly favorable. Besides, it 
was noticed that the pre-exponcntial volume factor in the 



TABLE I Number N^+_^- of e^-e~ pairs produced by differ- 
ent numbers n of laser pulses, with a total energy W of 10 kJ. 
The threshold value total energy Wth needed to produce one 
e^-e~ pair is shown in the third column. The precise collision 
geometry and the pulse parameters can be found in Bulanov 
et al., 2010a. Adapted from Bulanov et at, 2010a. 



_n 7V,+.,- at yy = 10 kJ Wth [kJ] 

2 < 10"^* 40 

4 < 10"* 20 

8 4.0 10 

16 1.8 X 10^ 8 

24 4.2 X 10® 5.1 



pair creation probability can be very large and partially 
compensate for the exponential suppression (Narozhny 
et ai, 2006). 

Fine details of pair production in a time-dependent os- 
cillating electric field are being studied nowadays because 
they might serve as characteristic signatures to discrim- 
inate the process of interest from potentially stronger 
background processes. For example, it was found that 
the momentum spectrum of the created particles is highly 
sensitive to a subcycle structure of the field (Hebenstreit 
et ai, 2009) and that in the presence of an alternating- 
sign time-dependent electric field, coherent interference 
effects are observed in the Schwinger mechanism (Akker- 
mans and Dunne, 2012). The observation in Hebenstreit 
et al, 2009 found an elegant mathematical explanation 
via the Stokes phenomenon (Dumlu and Dunne, 2010). 
Further effects stemming from the precise shape of the 
external field were analyzed in Dumlu, 2010 and Dumlu 
and Dunne, 2011. Also, the oscillating dynamics of the 
e'^-e" plasma created by a uniform electric field, in- 
cluding backreaction effects, was investigated (Apostol, 
2011; Benedetti et ai, 2011; Han et ai, 2010) (sec also 
Bialynicki-Birula and Rudnicki, 2011 and Kim and Schu- 
bert, 2011). 

In addition to pair creation in superstrong laser pulses 
of low frequency, the process is also extensively discussed 
in connection with the upcoming XFEL facilities (see, 
e.g., Alkofer et al, 2001 and Ringwald, 2001). Here the 
question arises as to what extent the spatial field depen- 
dence may influence the pair creation process, both in 
terms of total probabilities and particle momentum dis- 
tributions. According to Noether's theorem, pair produc- 
tion in a time-dependent oscillating electric field occurs 
with conservation of the total momentum, as well as of 
the total spin. The problem therefore reduces effectively 
to a two-level system since the field couples negative and 
positive-energy electron states of same momentum and 
spin only. The production process exhibits resonance 
when the energy gap is an integer multiple of the laser fre- 
quency, leading to a characteristic Rabi hopping between 
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FIG. 19 (Color online) Probability spectrum of e^-e~ pair 
production in two counterpropagating laser pulses, with the 
laser magnetic field included (black triangles) and neglected 
(red crosses). In the first case, the labeling {£r-£i) signifies the 
number of absorbed photons from the right-left propagating 
wave; in the second case, the peak labels denote the total pho- 
ton number {£). A vanishing initial momentum (i.e., positron 
momentum) and = 1 have been assumed. Adapted from 
Ruf et ai, 2009. 

the negative and positive-energy Dirac continua (Popov, 
1971). Due to the electron dressing by the oscillating 
field, the resonant laser frequencies are determined by 
the equation Iljq = 2(e), where (e) is the time- averaged 
electron energy in the time-dependent oscillating electric 
field. Accordingly, when the particle momentum is var- 
ied, several resonances occur corresponding to different 
photon numbers £. This gives rise to a characteristic ring 
structure in the momentum distribution (Mocken et at, 
2010). 

Modifications of these well-established properties of 
the pair creation process, when the spatial field depen- 
dence and, thus, the laser magnetic-field component are 
accounted for, have been revealed in Ruf et al., 2009. 
Utilizing an advanced computer code for solving the cor- 
responding Dirac equation numerically, it was shown that 
the positions of the resonances are shifted, several new 
resonances occur, and the resonance lines are split due to 
the influence of the spatial field dependence (see Fig. 19). 
The basic reason for these effects is that, in contrast to 
a uniform oscillating electric field, the photons in the 
counterpropagating laser pulses carry momentum along 
the beam axis. Therefore not only the total number £ 
of absorbed laser photons matters, but also how many 
of them have originated from the laser pulse travelling 
to the right and left, respectively. For example, for the 
multiphoton order £ = 5 two different resonance frequen- 
cies exist now, corresponding to £r = 3, £[ = 2 on the 
one hand, and to £r = 4, = 1, on the other. Due 
to the photon momentum, the former two-level scheme 
is also broken into a V^-type three-level scheme. This 
causes a splitting of the resonance lines, in analogy with 
the Autler-Townes effect known from atomic physics. 



For another numerical approach to space-time depen- 
dent problems in quantum field theory, we refer to the 
review Cheng et al, 2010. Moreover, Schwinger pair pro- 
duction in a space-time dependent electric field pulse has 
been treated very recently within the Wigner formalism 
(Hebenstreit et ai, 2011). Here, a self-bunching effect 
of the created particles in phase space, due to the spa- 
tiotemporal structure of the pulse, was found. 

Finally, we mention that, unlike a plane- wave field, a 
spatially focused laser beam is capable to produce e~^-e~ 
pairs from vacuum and this process has been investigated 
in Narozhny et ai, 2004 for different field polarizations. 
Spontaneous pair production may, in principle, also occur 
in a nuclear field for charge numbers Z exceeding a crit- 
ical value Zc, which depends on the nuclear model. For 
example, Zc = 173 for a uniformly charged sphere with 
radius 1.2 x 10"-^^ cm (Berestetskii et ai, 1982). See 
the reviews Baur et ai, 2007 and Zeldovich and Popov, 
1972 for more detailed information also on e^-e^ pair 
production in heavy ion collisions. 



D. Spin effects and other fundamental aspects of 
laser-induced pair creation 

A particularly interesting aspect of tunneling pair cre- 
ation is the electron and positron spin-polarization. In 
general, pronounced spin signatures in a field-induced 
process may be expected when the background field 
strength approaches the critical value i^cr (Kirsebom 
et al, 2001; Walser et ai, 2002). This indeed coincides 
with the condition for a sizable yield of tunneling pair 
production. Studies of spin effects in pair production by 
a high-energy photon and a strong laser field were per- 
formed in Tsai, 1993 and Ivanov et ai, 2005, based on 
considerations on helicity amplitudes and on the spin- 
polarization vector, respectively. Characteristic differ- 
ences between fermionic and bosonic particles have been 
revealed with respect to pair creation in an oscillating 
electric field (Popov, 1972) and in recent studies of the 
Klein paradox (Cheng et ai, 2009b; Krekora et al, 2004; 
Wagner et ai, 2010a). In the latter case it was shown 
that the existence of a fermionic (bosonic) particle in the 
initial state leads to suppression (enhancement) of the 
pair production probability due to the different quan- 
tum statistics. The enhancement in the bosonic case may 
be even exponential due to an avalanche process (Wag- 
ner et ai, 2010b). Concerning tunneling pair creation in 
combined laser and nuclear Coulomb fields, it has been 
shown that the internal spin-polarization vector is pro- 
portional ill nia,gnitudG to Xo,n 

and, to leading order, 
directed along the transverse momentum component of 
the electron (Di Piazza et al, 2010c). A helicity analy- 
sis of pair production in laser-proton collisions revealed 
that: 1) right-handed leptons are emitted in the labora- 
tory frame under slightly smaller angles with respect to 
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the proton beam than left-handed ones; 2) the rate of 
pair creation of spin- 1/2 particles exceeds by almost one 
order of magnitude the corresponding quantity for spin-0 
particles (Miiller and Miiller, 2011). 

Other fundamental aspects of laser-induced pair cre- 
ation have been recently investigated as well. They com- 
prise various kinds of e^-e~ correlations (Fedorov et al., 
2006; Krajewska and Kamihski, 2008; Krekora et al., 
2005), multiple pair creation (Cheng et al., 2009a), ques- 
tions of locality (Cheng et a/., 2008) and vacuum decay 
times (Labun and Rafelski, 2009), and consistency re- 
strictions on the maximum laser field strength to guar- 
antee the validity of the external-field approximation 
(Gavrilov and Gitman, 2008). 



IX. QED CASCADES 

As it was discussed in the previous section, the E-144 
experiment at SLAG is the only one, so far, where laser- 
driven multiphoton e^-e^ pair production has been ob- 
served. Considering that about 100 positrons have been 
detected in 22000 shots, each comprising the collision of 
about 10^ electrons with the laser beam, the process re- 
sults to be rather inefficient. One could attribute this 
to the relatively low intensity Iq of the laser system of 
1.3 X 10^® W/cm^ (^0 ~ 0.4 as the laser photon en- 
ergy was cjQ = 2.4 cV). However, the extremely high 
energy Eq of the electron beam (about 46.6 GeV) en- 
sured that that the nonlinear quantum parameter xo was 
about unity (xo ~ 0.3). A recent investigation (Sokolov 
et al., 2010a) has pointed out in general that in the 
mentioned setup, i.e., an electron beam colliding with 
a strong laser pulse, RR effects prevent the development 
of a cascade or avalanche process with an eflacient, pro- 
lific production of e+-e~ pairs even at much larger laser 
intensities such that 3> 1. By an avalanche or cas- 
cade process we mean here a process in which the in- 
coming electrons emit high-energy photons in the laser 
field, which can interact with the field itself generating 
e^-e" pairs, which, in turn, emit photons again and so 
on (of course a cascade process may also be initiated by 
a photon beam rather than by an electron beam). The 
above result has been obtained by numerically integrat- 
ing the kinetic equations, which describe the evolution of 
the electron, the positron and the photon distributions in 
a plane-wave background field from a given initial elec- 
tron distribution, and by accounting for the two basic 
processes that couple these distributions, i.e., multipho- 
ton Compton scattering and multiphoton Breit- Wheeler 
pair production. The physical reason why an avalanche 
process cannot develop in a single plane-wave field can 
be understood in the following way. In the ultrarelativis- 
tic case ^ li the above-mentioned basic processes in 
a plane-wave field are essentially controlled by the pa- 
rameters Xo and kq, respectively (see also Sees. III.B, 



V.A and VIII). Now, at the j*** step in which an elec- 
tron/positron emits a photon or a photon transforms into 
an e+-e~ pair, the initial quantity p'^^ or fcg''l is con- 
served and it is distributed over the two final particles (an 
electron/positron and a photon in multiphoton Comp- 
ton scattering and an e^-e^ pair in multiphoton Breit- 
Wheeler pair production). Thus, both resulting particles 
at each step will have a value of their own parameter 
X^''' = {p'o^^l/m){Eo/F„) or = / m){E^ / F„) 

smaller than that of the incoming particle. Moreover, 
due to the special symmetry of the plane-wave field, the 
quantities p\^\ or feg"*! are also rigorously conserved be- 
tween two steps (see also Sec. III. A). Then, the avalanche 
ends when the parameters x\^l a-nd >c^} at a certain step 
k are smaller than unity for i S [1, . . . , Nk], with Nk be- 
ing the number of particles at that step. 

The question arises as to whether other field configura- 
tions exist, where an avalanche process can be efficiently 
triggered (see the review Aharonian and Plyasheshnikov, 
2003 for the development of QED cascades in matter, 
photon gas and magnetic field). A positive answer to 
this question has first been given in Bell and Kirk, 2008: 
even the presence of a single electron initially at rest 
in a standing wave generated by two identical counter- 
propagating circularly polarized laser fields can prime an 
avalanche process already at field intensities of the order 
of 10^^ W/cm^. We note that in the presence of a single 
plane wave the same process would require an intensity 
of the order of !„, because for an electron initially at rest 
Xo = Eq/Fc^. From this point of view, the authors of Bell 
and Kirk, 2008 explain qualitatively the advantage of em- 
ploying two counterpropagating laser beams by means of 
an analogy taken from accelerator physics: a collision be- 
tween two particles in their center-of-momenta is much 
more efficient than if one of the particles is initially at 
rest, because much more of the initial energy can be 
transferred, for example, to create new particles. The 
authors approximated the standing wave by a rotating 
electric field (see Sec. VIII. C). In such a field and for 
an ultrarelativistic electron the controlling parameter is 
Xo = {p±/'m'){Eo/Fcr), where is the component of the 
electron momentum perpendicular to the electric field. 
By estimating p± ^ m^o (see also Eq. (2)), one ob- 
tains Xo ~ ^o[10^'^ W/cm^]/a;o[eV] (here ujq, Eq and Jq 
are the standing-wave's angular frequency, electric field 
amplitude and intensity). The investigation in Bell and 
Kirk, 2008 is based on the analysis of the trajectory of 
the electron in the rotating electric field including RR 
effects via the LL equation. Since the momentum of the 
electron oscillates around a value of the order of m^O: the 
electron emits high-energy photons efficiently that can in 
turn trigger the cascade (see Fig. 20). The possibility of 
describing the evolution of the electron via its classical 
trajectory, can be justified as follows. When an elec- 
tron interacts with a background electromagnetic field 
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FIG. 20 (Color online) The number of e^-e~ pairs (N±) and 
the number of photons (A'^) created by an initial single elec- 
tron in a rotating electric field as a function of the field inten- 
sity. The other plotted quantities are described in Bell and 
Kirk, 2008. From Bell and Kirk, 2008. 



like that of a laser, quantum effects are essentially of two 
kinds (Baler tt al., 1998): the first one is associated with 
the quantum nature of the electron motion and the sec- 
ond one with the recoil undergone by the electron when it 
emits a photon. For an ultrarelativistic electron it can be 
shown that, while the first kind of quantum effects is neg- 
ligible, the second kind is large and has to be taken into 
account (Baler et al, 1998). Thus, the basic assumption 
is that, since the background laser field is strong, the elec- 
tron is promptly accelerated to ultrarelativistic energies, 
the motion between two photon emissions is essentially 
classical and, if necessary, only the emissions have to be 
treated quantum mechanically by including the photon 
recoil. On the other hand, photons are assumed to prop- 
agate in the field along straight lines. 

The model employed in Bell and Kirk, 2008 was im- 
proved in Kirk et al. , 2009 by considering colliding pulsed 
fields with finite time duration and a realistic represen- 
tation for the synchrotron spectrum emitted by a rela- 
tivistic electron. The results in Bell and Kirk, 2008 were 
essentially confirmed and numerical simulations with lin- 
early polarized beams have shown a general insensitiv- 
ity of the cascade development to the polarization of 
the beams. Another interesting finding in Kirk et al., 
2009 is that the electrons in the standing wave tend to 
migrate to regions where the electric field vanishes and 
then they do not contribute to the pair-production pro- 
cess anymore. In both papers Bell and Kirk, 2008 and 
Kirk et al., 2009 the emission of radiation by the electron 
was treated classically, i.e., the electron was supposed 
to loose energy and momentum continuously although 
in Kirk et al., 2009 the damping force-term in the LL 



equation was evaluated by employing the total emitted 
power calculated quantum mechanically. The stochastic 
nature of the emission of a photon has been taken into 
account in Duclous et al., 2011. Analogously, the energy 
of the emitted photon is chosen randomly following the 
synchrotron spectral distribution and the momentum of 
the photon is always chosen to be parallel to that of the 
emitting electron. By contrast, in the pair production 
process by a photon, since the photon is not deflected 
by the laser field, it is assumed that after it has propa- 
gated one wavelength in the field, it decays into an e~^-e~ 
pair. The main result in Duclous et al, 2011 is that at 
relatively low intensities of the order of 10^'^ W/cm^ the 
pair production rate is increased if the quantum nature 
of the photon emission is taken into account. The reason 
is that, due to the stochastic nature of the emission pro- 
cess, the electron can propagate for an unusually large 
distance before emitting. In this way it may gain an un- 
usually large amount of energy and consequently emit a 
high-energy photon, that can be more easily converted 
into an e+-e^ pair. Moreover, the discontinuous nature 
of the (curvature of the) electron trajectory is shown to 
slow down the tendency of the electrons to migrate to 
regions where the electric field vanishes. 

The intensity threshold of the avalanche process in a 
rotating electric field has also been investigated in Fedo- 
tov et al, 2010. Denoting by tacc the time an electron 
needs to reach an energy corresponding to xo = 1 start- 
ing from rest in the given field, by ig (^7) the electron 
(photon) lifetime under photon emission {e~^-e~ pair pro- 
duction) and by tgsc the time after which the electron es- 
capes from the laser field, the authors give the following 
conditions for the occurrence of the avalanche process: 
^acc ^ te,tj <^ tcsc- Estimates based on the classical 
electron trajectory without including RR effects, lead to 
the simple condition Eq > aFcr for the avalanche to be 
primed in an optical field. The above estimate corre- 
sponds to an intensity of about 2.5 x 10^^ W/cm^, i.e., 
one order of magnitude larger than what was found in 
Bell and Kirk, 2008. However, the main result of Fedo- 
tov et al., 2010 concerns the limitation, brought about 
on the maximal laser intensity that can be produced 
in the discussed field configuration by the starting of 
the avalanche process. In fact, the energy to acceler- 
ate the electrons and the positrons participating in the 
cascade has to come from the background electromag- 
netic field. By assuming an exponential increase of the 
number of electrons and positrons, it is found that al- 
ready at laser intensities of the order of 10^^ W/cm^ the 
created electrons and positrons have an energy which ex- 
ceeds the initial total energy of the laser beams. This 
hints at the fact that at such intensities the colliding 
laser beams are completely depleted due to the avalanche 
process. The results obtained from qualitative estimates 
in Fedotov et al., 2010 have been scrutinized in Elkina 
et al., 2011 by means of more realistic numerical meth- 
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ods based on kinetic or cascade equations. In general, 
if f:f{r,p,t) {f^{r,k,t)) is the electron/positron (pho- 
ton) distribution function (upper and lower sign for elec- 
tron and positron, respectively) in the phase-space {r,p) 
((r,fc)) and e = m? + [uj — \k\), their evolution 
in the presence of a given classical electromagnetic field 
{E{r,t), B{r,t)) is described by the kinetic equations 



y dfc w^i,d{r,p + k ^ k,t)f^{r,p + k,t) 
-/t(^'P'*) J dkwy.s.dir,p^ k,t) 
+ J dk Wcrcir,k ^ p,t)f^{r,k,t), 

: J dpw,,d{r,p^k,t)[U{r,p,t) + f^{r,p,t)] (49) 
-fy{r,k,t) J dpwcrcir,k ^ p,t), 



(48) 



where F^^r^p, t) = e[E{r, t) + {p/e) x B{r, t)] and where 
Wra,d{r,p — k,t) {wcre{r,k — )• p,t)) is the local prob- 
ability per unit time and unit momentum that an elec- 
tron/positron (photon) with momentum p (fc) emits (cre- 
ates) at the space-time point {t, r) a photon with momen- 
tum k (an e~^-e~ pair with the electron/positron having 
a momentum p). It is worth pointing out here a con- 
nection between the development of a QED cascade and 
the quantum description of RR. In fact, as has been dis- 
cussed in Sec. VI, from a quantum point of view, RR cor- 
responds to the multiple recoils experienced by the elec- 
tron in the incoherent emission of many photons. Thus, 
in the kinetic approach RR is described by those terms in 
Eqs. (48) and (49), which do not involve e~^-e~ pair pro- 
duction. In fact, in Elkina et at, 2011 it has been shown 
in the ultrarelativistic case that the equation of motion 
for the average momentum of the electron distribution, 
as derived from Eq. (48), coincides with the LL equation 
in the classical regime xo ^ 1 (note that pair production 
is exponentially suppressed at xo 1). 

In Elkina et ai, 2011 the background field is approx- 
imated as a uniform, rotating electric field E{t) and 
f±{r,p,t) ^ f±ip,t) {f^{r,k,t) ^ f-r{k,t)). Analo- 
gously to Duclous et ai, 2011, it is assumed that the 
momentum of the photon emitted in multiphoton Comp- 
ton scattering is parallel to that of the emitting electron 
(positron); in the same way, the momenta of the elec- 
tron and positron created in multiphoton Breit- Wheeler 
pair production are assumed to be parallel to that of the 
creating photon. The evolution of the electron, positron 
and photon distributions has been investigated by nu- 
merically integrating the resulting cascade equations via 



a Monte Carlo method. Whereas, the instants of radia- 
tion and pair production have been randomly generated. 
In particular, the exponential increase of the number of 
e"'"-e~ pairs and the qualitative estimate, for example, of 
the typical energy of the electron at the moment of the 
photon emission carried out in Fedotov et ai, 2010 have 
been confirmed (apart from discrepancies within one or- 
der of magnitude). 

In both papers Elkina et ai, 2011 and Fedotov et ai, 
2010 only the case of a rotating electric field was con- 
sidered. In Bulanov et ai, 2010b it is pointed out that 
the limitation on the maximal laser intensity reachable 
before the cascade is triggered, strongly depends on the 
polarization of the laser beams which create the stand- 
ing wave. The paradigmatic cases of a rotating electric 
field and of an oscillating electric field are compared. 
The estimates presented in Bulanov et ai, 2010b for the 
case of a rotating electric field essentially confirm that 
the avalanche starts at laser intensities of the order of 
lO^s W/cm^. The main physical reason why the cascade 
process in a circularly polarized standing wave starts at 
such an intensity is that in a rotating electric field the 
electron emits photons with typical energies of the order 
of 0.29 a;o7'^; i-c, proportional to the cube of the Lorentz 
factor of the emitting electron 7 (Bulanov et ai, 2010b). 
Whereas, in an oscillating electric field the typical emit- 
ted energy scales as 7^, such that in order to radiate a 
hard photon with a given energy, a much more energetic 
electron is needed. Hence, the authors of Bulanov et ai, 
2010b conclude that in an oscillating electric field RR and 
quantum effects do not play a fundamental role at laser 
intensities smaller than Icr and that avalanche processes 
do not constitute a limitation. It is crucial however, for 
the conclusion in Bulanov et ai, 2010b that the collision 
of the laser beams occurs in vacuum, i.e., the seed elec- 
trons and positrons which would trigger the cascade are 
supposed to be created in the collision itself. 

The question of the occurrence of the avalanche for two 
colliding linearly polarized pulses has also been addressed 
in Nerush et ai, 2011b, where a detailed description of 
the system under investigation has been provided. In 
fact, previous models had assumed the background elec- 
tromagnetic field as given, neglecting in this way the field 
generated by the electrons and positrons. The approach 
followed in Nerush et al, 2011b exploits the existence of 
two energy scales for the photons: one is that of the ex- 
ternal laser field and of the plasma fields which is much 
smaller than m. and the other is that of the photons 
produced by the high-energy electrons which is, by con- 
trast, much larger than m. The evolution of the low- 
energy photons is described by means of Maxwell's equa- 
tions which are solved with a PIC code, i.e., the photons 
are treated as a classical electromagnetic field. Whereas, 
the production of hard photons as well as the creation 
of e~^-e~ pairs is described as a stochastic process em- 
ploying a Monte Carlo method. Unless low-energy ones. 
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FIG. 21 (Color online) Snapshot of the normalized electron 
density pe (part a)), of the normalized photon density p~, (part 
b)), and of the normalized laser intensity pi (part c)) 25.5 laser 
periods after the two laser counterpropagating beams col- 
lided. The normalized density of positrons is approximately 
the same as that of the electrons. Also, the intensity of each 
colliding beam is Jo = 3 x lO'^* W/cm'^ and the common wave- 
length is Ao = 0.8 ^m. Adapted from Nerush et al., 2011b. 



hard photons are treated as particles and their evolution 
is described via a distribution function. It is surpris- 
ing that by considering a single seed electron initially 
at rest at a node of the magnetic field of a linearly po- 
larized standing wave, an avalanche process is observed 
in the numerical simulation already for a laser intensity 
/o = 3 X 10^^ W/cm^ and a laser wavelength Aq = 0.8 ^m 
(see Fig. 21). The figure clearly shows the formation of 
an overdense e^-e~ plasma in the central region |a;| < Aq. 
In the same numerical simulations it is found that after 
about 20 laser periods almost half of the initial energy of 
the laser field has been transferred to the plasma. Dis- 
agreement with the predictions in Bulanov et al., 2010b 
is stated to be due to the formation of the avalanche in 
regions between the nodes of the electric and magnetic 
field, where the simplified analysis of the electron motion 
carried out in Bulanov et al., 2010b is not valid. Fur- 
ther analytical insight into the formation of the cascade 
has been reported in Nerush et al., 2011a by analyzing 
approximate solutions of the cascade equations in the 
presence of a rotating electric field. 



X. MUON-ANTIMUON AND PION-ANTIPION PAIR 
PRODUCTION 

The production of e"'"-e~ pairs in strong laser fields 
has been discussed in Sec. VIII. In view of the ongo- 
ing technical progress the question arises as to whether 
also heavier particles such as muon-antimuon (^~^-fi~) 
or pion-antipion (tt^-tt") pairs can be produced with the 
emerging near-future laser sources. The production of 
fJ-~^-fJ.~ pairs from vacuum in the tunneling regime ap- 
pears rather hopeless, though, since the required field 
needs to be close to Fcr.p. — gfiFd- = 5.6 x 10^" V/cm, 
with the ratio = m^/m ~ 207 between the muon 
mass and the electron mass to. Even by boosting 
the effective laser fields with the Lorentz factors (^ 10^) 
of the most energetic electron beams available (Bamber 
et al, 1999), the value of Fc,,^ seems out of reach. The 
tunneling production of 7r"*'-7r~ pairs is even more diffi- 
cult as = m^/m « 273. However, and tt^-tt^ 
production can occur in microscopic collision processes 
in laser-generated or laser-driven plasmas, as well as by 
few-photon absorption from a high-frequency laser wave. 



A. Muon-antimuon and plon-antlpion pair production in 
laser-driven collisions in plasmas 

Energetic particle collisions in a plasma can in prin- 
ciple drive /i"*"-^" and tt^-tt^ production. The plasma 
may consist either of electrons and ions or of electrons 
and positrons. Both kinds of plasmas can be produced by 
intense laser beams interacting with a solid target. With 
respect to e~^-e~ plasmas, this has been predicted by 
Liang et al., 1998. As has been mentioned in Sec. VIII. 
abundant amounts of e'^-e~ pairs have been recently pro- 
duced in this manner at LLNL with pair densities of 
the order of 10^*^ cm'^ (Chen et al., 2009, 2010) and 
much higher densities of the order of 10^^ cm^^ have been 
also predicted (Shen and Meyer-ter-Vehn, 2001). Theo- 
reticians have therefore started to investigate the prop- 
erties and time evolution of relativistic e^-e~ plasmas 
(Aksenov et al, 2007; Hu and Miiller, 2011; Kuznetsova 
et al., 2008; Kuznetsova and Rafelski, 2012; Mustafa and 
Kampfer, 2009; Thoma, 2009a, b). In particular, it has 
been shown (Kuznetsova et al, 2008; Thoma, 2009a, b) 
that in an e"'"-e~ plasma of 10 MeV temperature, 
pairs, 7r"'"-7r~ pairs as well as neutral 7r° can be created 
in e^-e^ collisions. The required energy stems from the 
high-energy tails of the thermal distributions. 

Also cold e+-e^ plasmas of high density can be gener- 
ated nowadays due to dedicated positron accumulation 
and trapping techniques (Cassidy et al, 2005). When 
such a nonrelativistic low-energy plasma interacts with a 
superintense laser field, pair production can oc- 

cur as weU (Miiller et al, 2006, 2008b). In this case, 
the plasma particles acquire the necessary energy by 
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strong coupling to the external field which drives the elec- 
trons and positrons into violent collisions. The minimum 
laser peak intensity to ignite the reaction e+e^ — > jX^ 
amounts to about 7 x 10^^ W/cm^ at a typical optical 
laser photon energy of = 1 eV, corresponding to 
^o,min = Qfi ~ 207. The rate Re+e-^fj.+ fj,- of the pro- 
cess in the presence of a linearly polarized field reads 

with the number N± of electrons/positrons and the in- 
teraction volume which is determined by the laser fo- 
cal spot size. Equation (50) may be made intuitively 
meaningful by introducing the invariant cross section 
'^e+e-^n+i.i- of production in an e+-e~ collision 

in vacuum (Peskin and Schroeder, 1995): 

where the upper index * indicates quantities in the 
center-of-mass system of the colliding electron and 
positron, as, e.g., the common electron and positron en- 
ergy s* . Now, by exploiting the fact that the quan- 
tity i?e+-e- /y is a Lorentz invariant and that in the 
present physical scenario e* can be estimated as m^Oj 
Eq. (50) implies the usual relation /l^* ~ 

(Te+e-->/i+;i-'^+^- between the number of events per unit 
volume and per unit time R*+^^_^^+^- /V* and the cross 
section crg+g-_j.^+^- ; n^. ~ N±/V* denote here the par- 
ticle densities. The process e"*"e~ — >■ fJ.'^iJ.~ in the pres- 
ence of an intense laser wave has also been considered in 
Nedoreshta et al, 2009. 

In laser-produced electron-ion plasmas resulting from 
intense laser-solid interactions, and tt^-tt" pairs 

can be generated by the cascade mechanism via energetic 
bremsstrahlung, like in the case of e^-e~ pair produc- 
tion mentioned in Sec. VIII. Assuming a laser-generated 
few-GeV electron beam, several hundreds to thousands 
of fJ-'^-n~ pairs arise from bremsstrahlung conversion in 
a high-Z target material (Titov et ai, 2009). The pro- 
duction of 7r+-7r^ pairs by laser-accelerated protons was 
considered in Bychenkov et ai, 2001, where a threshold 
laser intensity of lO'^^ W/cm^ for the process to occur 
was determined. 

B. Muon-antimuon and pion-antipion pair production in 
high-energy XFEL-nucleus collisions 

In this subsection another mechanism of and 
TT~^-TT~ pair creation by laser fields will be pursued, which 
is based on the collision of an x-ray laser beam with an 
ultrarelativistic nuclear beam. This setup is similar to 
the one of Sec. VIII.B. 



In the case of pair creation, by considering an 

x-ray photon energy of wq = 12 keV and a nuclear rela- 
tivistic Lorentz factor of 7„ = 7000, the photon energy in 
the rest-frame of the nucleus amounts to ujq « 27„a;o = 
168 MeV. The energy gap of 2m^ for fi^-fi^ pair pro- 
duction can thus be overcome by two-photon absorption 
(Miiller et al., 2008a, 2009a). Note that because of pro- 
nounced recoil effects, the Lorentz factor which would be 
required for two-photon fi^-^J.^ production by a projec- 
tile electron is much larger: 7 > 10^, corresponding to a 
currently unavailable electron energy in the TeV range. 

At first sight, e+-e^ and A'^'M^ pa-ir production in 
combined laser and Coulomb fields seem to be very sim- 
ilar processes since the electron and muon only differ by 
their mass (and lifetime). In this picture, the corre- 
sponding production probabilities would coincide when 
the laser field strength and frequency are scaled in accor- 
dance with the mass ratio g^, i.e., P^+_^-{Eq,jj_,ujq^^) = 
Pe+-e- (^^o,wo) for £^0,^ = P^^o and wq,^ = gi^WQ. This 
simple scaling argument does not apply, however, as the 
large muon mass is connected with a correspondingly 
small Compton wavelength A^^^ = ^c/Ofi ~ l.SG fm 
(1 fm = lO"^'^ cm), which is smaller than the radius of 
most nuclei. As a result, while the nucleus can be ap- 
proximately taken as pointlike in e'^-e" pair production 
(Ac ~ 386 fm), its finite extension must be taken into 
account in fi^-fJ.^ pair production. Pronounced nuclear 
size effects have also been found for i-i^-j-i^ production in 
relativistic heavy-ion collisions (Baur et at, 2007). 

Muon pair creation in XFEL-nucleus collisions can be 
calculated via the amplitude in Eq. (44), with the nu- 
clear potential Vn{r) arising from an extended nucleus. 
It leads to the appearance of a nuclear form factor F{q^) 
which depends on the recoil momentum q. For example, 
F(q^) ~ exp(— q'^a^/6) for a Gaussian nuclear charge dis- 
tribution of root-mcan-square radius a. Since the typi- 
cal recoil momentum is q ^ m^, the form factor leads to 
substantial suppression of the process. The fully differen- 
tial production rate dR^+_^- = dR^^^_^_ + dR^^I'^^^^ may 
be split into an elastic and an inelastic part, depending 
on whether the nucleus remains in its ground state or 
gets excited during the process. They read dR^'^^ _ = 
dB^_^_Z^F-{^) ^nAdR^^l « ^[l-F^q^)], 

respectively, with dR^^l_^_ being the production rate for 
a pointlike proton. 

Figure 22 shows total production rates R^+_^^ 

in the rest-frame of the nucleus for several nuclei col- 
liding with an intense XFEL beam. For an extended 
nucleus, the elastic rate increases with its charge but de- 
creases with its size. This interplay leads to the emer- 
gence of maximum clastic rates for medium-heavy ions. 
Figure 22 also implies that the total rate = 
r'"'^!^ _ -I- in the laboratory frame saturates at 

high Z values since increases with nuclear charge. 
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FIG. 22 Total rates for fi^-jJ." pair creation by two-photon 
absorption from an intense XFEL beam {ujq — 12 keV, lo ~ 
2.5 X 10^^ W/cm^) colliding with various ultrarelativistic nu- 
clei (7„ = 7000). The triangles show elastic rates, whereas 
the squares indicate total ("elastic -I- inelastic") rates. The 
numerical data are connected by fit curves. The dotted line 
holds for a pointlike nucleus. The production rates are calcu- 
lated in the rest-frame of the nucleus. Adapted from Miiller 
et al, 2008a. 



For highly-charged nuclei the main contribution stems 
from the inelastic channel vifhere the protons inside the 
nucleus act incoherently Z). This implies 

that despite the high charges, high-order Coulomb cor- 
rections in Za are of minor importance. In the collision, 
also e+-e~ pairs are produced by single-photon absorp- 
tion in the nuclear field. However, this rather strong 
background process does not deplete the x-ray beam. 

In XFEL-proton collisions, tt^-tt^ pairs can be gener- 
ated as well. A corresponding calculation has been re- 
ported in Dadi and Miiller, 2011, which includes both 
the electromagnetic and hadronic pion-proton interac- 
tions. The latter was described approximately by a phe- 
nomenological Yukawa potential. It was shown that, 
despite the larger pion mass, tt^-tt^ pair production 
by two-photon absorption from the XFEL field largely 
dominates over the corresponding process of pair 
production in the Doppler-boosted frequency range of 
uJq « 150-210 MeV. This dominance is due to the much 
larger strength of the strong (hadronic) force compared 
with the electromagnetic force. As a consequence, in this 
energy range fJ.~^-fJ.^ pairs are predominantly produced 
indirectly via two-photon tt^-tt^ production and subse- 
quent pion decay, 7r+ — > /i+ + i^,, and tt^ — > /i^ + D^. 

In relativistic laser-nucleus collisions, ^i^-^i" or 7r"*'-7r~ 
pairs can also be produced indirectly within a two-step 
process (Kuchiev, 2007). First, upon the collision, an 
e^-e" pair is created via tunneling pair production. Af- 
terwards, the pair, being still subject to the electromag- 
netic forces exerted by the laser field, is driven by the 
field into an energetic e^-e" collision. If the collision en- 
ergy is large enough, the reaction e+e^ — )■ may 



be triggered. This two-step mechanism thus represents a 
combination of the processes considered in Sec. VIII. B 
and X.A. Besides, it may be considered as a generaliza- 
tion of the well-established analogy between strong-field 
ionization and pair production (see Sec. VIII) to include 
also the recollision step. 

XI. NUCLEAR PHYSICS 

Influencing atomic nuclei with optical laser radiation 
is, in general, a difficult task because of the large nuclear 
level spacing of the order of 1 keV-1 MeV, which 
exceeds typical laser photon energies ~ 1 eV by orders 
of magnitude (Matinyan, 1998). Also the laser's electric 
work performed over the tiny nuclear extension r„ ^ 1-5 
fm is usually too small to cause any sizable effect. In fact, 
the requirement leli^oJ'n ^ can only be satisfied for at 
laser- field amplitudes at least close to F^^. Direct laser- 
nucleus interactions have therefore mostly been dismissed 
in the past. 

On the other hand, laser-induced secondary reactions 
in nuclei have been explored especially in the late 1990s. 
Via laser-heated clusters and laser-produced plasmas, 
various nuclear reactions have been ignited, such as fis- 
sion, fusion and neutron production. In all these cases, 
the interaction of the laser field with the target first 
produces secondary particles such as photo-electrons or 
bremsstrahlung photons which, in a subsequent step, 
trigger the nuclear reaction. For a recent review on this 
subject we refer to Ledingham and Galster, 2010 and 
Ledingham et al., 2003. 

In recent years, however, the interest in direct laser- 
nucleus coupling has been revived by the ongoing techno- 
logical progress towards laser sources of increasingly high 
intensities as well as frequencies. Indeed, when suitable 
nuclear isotopes are considered, intense high-frequency 
fields or superstrong near-optical fields may be capable 
of affecting the nuclear structure and dynamics. 

A. Direct laser-nucleus interaction 

1. Resonant laser-nucleus coupling 

There are several low-lying nuclear transitions in the 
keV range, and even a few in the eV range. Examples of 
the latter are 229Th (Af « 7.6 eV) and ^ssu (Af « 76 
eV) (Beck et al., 2007). These isotopes can be excited 
by the 5th harmonic of a Ti:Sa laser (wo = 1-55 eV) 
and by pulses envisaged at the ELI attosecond source 
(ELI, 2011), respectively. Even higher frequencies can 
be attained by laser pulse reflection from relativistic fly- 
ing mirrors of electrons extracted from an underdense 
plasma (Bulanov et al., 1994) or possibly also from an 
overdense plasma (Habs et al., 2008). Otherwise keV- 
encrgy photons are generated by XFELs, which, as we 
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have seen in Sec. TLB, are presently emerging as large- 
scale facilities, e.g., at SLAC (LCLS, 2011) and DESY 
(European XFEL, 2011), and which could be employed 
with focusing (Mimura et al., 2010) and reflection devices 
(Shvyd'ko et al, 2011). In addition, also XFEL facili- 
ties of table-top size (Griiner et al., 2007; Kneip et al., 
2010) and even fully coherent XFEL sources are envis- 
aged such as the future XFEL Oscihator (XFELO) (Kim 
et al, 2008) or the Seeded XFEL (SXFEL) (Feldhaus 
et al, 1997; LCLS II, 2011). Brilliant gamma-ray beams 
with spectra peaked between 20 KeV and 150 KeV have 
been recently produced from resonant betatron motion of 
electrons in a plasma wake (Cipiccia et al, 2011). A new 
material research center, the Matter-Radiation Interac- 
tions in Extremes (MaRIE) (MaRIE, 2011) is planned, 
allowing for both fully coherent XFEL light with photon 
energy up to 100 keV and accelerated ion beams. The 
photonuclear pillar of ELI to be set up near Bucharest 
(Romania) is planned to provide a compact XFEL along 
with an ion accelerator aiming for energies of 4-5 GeV 
(ELI, 2011). In addition, coherent gamma-rays reach- 
ing few MeV energies via electron laser interaction arc 
envisaged at this facility (ELI, 2011; Habs et al, 2009). 

With these sources of coherent high-frequency pulses, 
driving electric dipole (El) transitions in nuclei is becom- 
ing feasible (Biirvenich et al, 2006a). Table II displays 
a list of nuclei with suitable El transitions. Along with 
an appropriate moderate nucleus acceleration, resonance 
may be induced due to the Doppler shift via the factor 
(1 +Vn)ln in the counterpropagating setup, with and 
7„ being the nucleus velocity and its Lorentz factor, re- 
spectively (Biirvenich et al, 2006a). For example, with 
^^■^Ra and an XFEL frequency of 12.4 keV a factor of 
(1 +Vn)ln — 4 would be sufficient. In general such mod- 
erate pre-accelerations of the nuclei would be of great 
assistance since they increase the number of possible nu- 
clear transitions for the limited number of available light 
frequencies. Note that the electric field strength of the 
laser pulse transforms analogously in the rest-frame of 
the nucleus, such that the applied laser field in the labo- 
ratory frame may correspondingly be weaker for a coun- 
terpropagating setup. 

For optimal coherence properties of the envisaged high- 
frequency facilities, subsequent pulse applications were 
shown to yield notable excitation of nuclei (Biirvenich 
et al, 2006a). In addition, many low-energy electric 
quadrupole (E2) and magnetic dipole (Ml) transitions 
are available. Here it is interesting to note that certain E2 
or Ml transitions can indeed be competitive in strength 
with El transitions (Palffy, 2008; Palffy et al, 2008). 
While the majority of transitions is available for high 
frequencies in the MeV domain (requiring substantial nu- 
cleus accelerations). Fig. 23 displays also numerous suit- 
able nuclear transitions below 12.4 keV along with the 
excitation efhciencics from realistic laser pulses. While 
indeed experimental challenges are high, resonant direct 



TABLE II The transition energy A£^, the dipole moment /i, 
the life-time Tg (r^) of the ground (excited) state of few rel- 
evant nuclear systems and El transitions (Aas et al., 1999; 
NNDC, 2011). Adapted from Biirvenich et al, 2006a. 



Nucleus 


Transition 


[keV] 


/i[e fm] 




Te[ps] 


'^^Sm 


3/2^ 


->3/2+ 


35.8 


> 0.75=^ 


47 h 


< 100 


^**iTa 


9/2- 


^7/2+ 


6.2 


0.04=^ 


stable 


6 X 10'' 


Ac 


3/2+ 


^3/2" 


40.1 


0.24=^ 


10.0 d 


720 


223Ra 


3/2" 


^3/2+ 


50.1 


0.12 


11.435 d 


730 




3/2" 


^1/2+ 


37.9 


b 


18.68 d 


b 




5/2" 


^5/2+ 


186 


0.017 


25.52 h 


1030 



^ Estimated via the Einstein A coefficient from Te and AS 
^ Not listed in the National Nuclear Data Center (NNDC) 
(NNDC, 2011) 
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FIG. 23 (Color online) Number A'^. of signal photons per nu- 
cleus per laser pulse for several isotopes with first excited 
states below 12.4 keV (green squares) and above 12.4 keV 
(red crosses). The results are plotted versus the atomic 
number Z. The considered European XFEL has a pulse 
duration of 100 fs and an average brilliance of 1. 6x10^5 
photons/(s mrad^ mm^ 0.1% bandwidth) (European XFEL, 
2011). Adapted from Pallfy et al, 2008. 



interactions of laser radiation with nuclei is expected to 
pave the way for nuclear quantum optics. Especially con- 
trol in exciting and deexciting certain long-living nuclear 
states would have dramatic implications for nuclear iso- 
mer research (Aprahamian and Sun, 2005; Palffy et al, 
2007; Walker and Dracoulis, 1999). As an obvious ap- 
plication this would be of relevance for nuclear batter- 
ies (Aprahamian and Sun, 2005; Walker and Dracoulis, 
1999), i.e., for controlled pumping and release of energy 
stored in long-lived nuclear states. In atomic physics, the 
STimulated Raman Adiabatic Passage (STIRAP) tech- 
nique has proven to be highly efficient in controlling pop- 
ulations robustly with high precision (Bergmann et al, 
1998). On the basis of currently envisaged accelerators 
and coherent high-frequency laser facilities, it has been 
recently shown that such an efficient coherent population 
transfer will also be feasible in nuclei (Liao et al, 2011). 
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Most recently, a nuclear control scheme with optimized 
pulse shapes and sequences has been developed in Wong 
et ai, 2011. 

Serious challenges are certainly imposed by the nuclear 
linewidths that may either be too narrow to allow for suf- 
ficient interaction with the applied laser pulses or inhibit 
excitations and coherences due to large spontaneous de- 
cay. Decades of research in atomic physics allowing now 
for shaping atomic spectra via quantum interference (Ev- 
ers and Keitel, 2002; Kiffner et ai, 2010; Postavaru et ai, 
2011) raise hopes that such obstacles may be overcome 
in the near future as well. 

Direct photoexcitation of giant dipole resonances 
with few-MeV photons via laser-electron interaction was 
shown to be feasible (Weidenmiiller, 2011) based on en- 
visaged experimental facilities such as ELI. Finally, care 
has to be taken to compare the laser-induced nuclear 
channels with competing nuclear processes via, for ex- 
ample, bound electron transitions or electron captures in 
the atomic shells (Palffy, 2010; Palffy et ai, 2007). 



2. Nonresonant laser-nucleus interactions 

Already decades ago a lively debate was started on 
whether nuclear /3-decay may be significantly affected 
by the presence of a strong laser pulse or not (Akhme- 
dov, 1983; Becker et ai, 1984; Nikishov and Ritus, 1964b; 
Reiss, 1983) and a conclusive experimental answer to this 
issue is still to come. Most recently the notion of affect- 
ing nuclear a-decay with strong laser pulses has been 
discussed showing that moderate changes of such nuclear 
reactions with the strongest envisaged laser pulses are 
indeed feasible (Castaheda Cortes et al, 2012, 2011). 

When the laser intensity is high enough {Iq > 10^^ 
W/cm^) low- frequency laser fields are able to influence 
the nuclear structure without necessarily inducing nu- 
clear reactions. In such ultrastrong fields, low-lying nu- 
clear levels get modified by the dynamic (AC-) Stark shift 
(Biirvenich et ai, 2006b). These AC-Stark shifts are of 
the same order as in typical atomic quantum optical sys- 
tems relative to the respective transition frequencies. At 
even higher, supercritical intensities {Iq > 10^^ W/cm^) 
the laser field induces modifications to the proton root- 
mean-square radius and to the proton density distribu- 
tion (Biirvenich et ai, 2006b). 



B. Nuclear signatures in laser-driven atomic and molecular 
dynamics 

Muonic atoms represent traditional tools for nuclear 
spectroscopy by employing atomic physics techniques. 
Due to the large muon mass compared to that of the 
electron, w 207 m, and because of its correspond- 
ingly small Bohr radius as.fi = Ac./j/a ~ 255 fm, the 



muonic wave-function has a large overlap with the nu- 
cleus. Precise measurements of x-ray transitions be- 
tween stationary muonic states are therefore sensitive to 
nuclear-structure features such as finite size, deforma- 
tion, surface thickness, or polarization. 

When a muonic atom is subjected to a strong laser 
field, the muon becomes a dynamic nuclear probe which 
is periodically driven across the nucleus by the field. This 
can be inferred, for example, from the high-harmonic ra- 
diation emitted by such systems (Shahbaz et al., 2010, 
2007). Figure 24 compares the HHG spectra from muonic 
hydrogen versus muonic deuterium subject to a very 
strong XUV laser field. Such fields are envisaged at the 
ELI attosecond source (see Fig. 1). Due to the differ- 
ent masses A/„ of the respective nuclei (M„ ~ mp for a 
hydrogen nucleus and M„ « mp+nin for a deuterium nu- 
cleus by neglecting the binding energy, with irip/n being 
the proton/neutron mass), muonic hydrogen gives rise to 
a significantly larger harmonic cut-off energy. The reason 
can be understood by inspection of the ponderomotive 
energy 

which depends in the present case on the reduced mass 
Air = TTifiMn/{m^ + Mn) of the muon-nucleus system. 
The reduced mass of muonic hydrogen (« 93 MeV) is 
smaller than that of muonic deuterium (« 98 MeV) and 
this implies a larger ponderomotive energy and an en- 
larged plateau extension. The influence of the nuclear 
mass can also be explained by the separated motions of 
the atomic binding partners. The muon and the nucleus 
arc driven by the laser field into opposite directions along 
the laser's polarization axis. Upon recombination their 
kinetic energies sum up as indicated on the right-hand 
side of Eq. (52). Within this picture, the larger cut-off 
energy for muonic hydrogen results from the fact that, 
due to its smaller mass, the proton is more strongly ac- 
celerated by the laser field than the deuteron. 

Due to the large muon mass, very high harmonic cut- 
off energies can be achieved via muonic atoms with charge 
number Z in the nonrelativistic regime of interaction. 
Since the harmonic-conversion efficiency as well as the 
density of muonic atom samples are rather low, it is im- 
portant to maximize the radiative signal strength. A siz- 
able HHG signal requires efficient ionization on the one 
hand, and efficient recombination on the other. The for- 
mer is guaranteed if the laser's electric field amplitude 
Eq lies just below the border of over-barrier ionization, 

QeS 16 

Here, QcS ~ \e\{Z/Mn + l/m^)Mr represents an effec- 
tive charge (Reiss, 1979; Shahbaz et al., 2010). Efficient 
rccoUision is guaranteed if the magnetic drift along the 
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FIG. 24 (Color) HHG spectra emitted from muonic hydrogen 
(black) and muonic deuterium (red) in a laser field of intensity 
Iq w 10^^ W/cm^ and pfioton energy uio ~ 60 eV. "Arb. u." 
stands for "arbitrary units". From Shahbaz et ai, 2007. 



are quite small, however, because of the large difference 
between the laser photon energy and the nuclear tran- 
sition energy. Nuclear excitation can also be triggered 
by intense laser-induced recollisions of ficld-ionizcd high- 
energy electrons (Kornev and Zon, 2007; Milosevic et ai, 
2004; Mocken and Keitel, 2004). 

Finally, muonic molecules are of particular interest for 
nuclear fusion studies. Modifications of muon-catalyzed 
fusion in strongly laser-driven muonic D J molecules have 
been investigated (Chelkowski et ai, 2004; Paramonov, 
2007). It was found that applied field intensities of the 
order of 10^'^ W/cm^ can control the molecular recol- 
lision dynamics by triggering the nuclear reaction on 
a femtosecond time scale. Similar theoretical studies 
have recently been carried out on aligned (electronic) 
HT molecules, i.e., involving a Tritium atom (Zhi and 
Sokolov, 2009). 



laser propagation direction can be neglected. Equation 
(15) indicates that this is the case here, provided 

/Qcff^\% 16^ (54) 

V MrUJo J ^ y/2MrIp ■ 

The two above inequalities define a maximum laser in- 
tensity and a minimum laser angular frequency which are 
still in accordance with the conditions imposed. At these 
laser parameters, the maximum harmonic cut-off ener- 
gies are attained and an efficient ionization-recollision 
process is guaranteed. For muonic hydrogen the corre- 
sponding lowest frequency lies in the Vacuum Ultraviolet 
(VUV) range {uiq « 27 eV) and the maximum field in- 
tensity amounts to 1.6 x 10^"^ W/cm^. At these values, 
the harmonic spectrum extends to a maximum energy of 
approximately 0.55 MeV. For light muonic atoms with 
nuclear charge number Z > 1, the achievable cut-off 
energies are even higher, reaching several Me Vs. This 
holds prospects for the production of coherent ultrashort 
gamma-ray pulses (see also Xiang et ai, 2010). 

In principle, also nuclear size effects arise in the HHG 
signal from muonic atoms. This has been shown qualita- 
tively in ID numerical simulations, where a 50% enhance- 
ment of the harmonic plateau height has been obtained 
for muonic hydrogen compared with muonic deuterium 
(Shahbaz et ai, 2010, 2007). This has been attributed to 
the enhanced final muon acceleration towards the hydro- 
genic core. For more precise predictions, 3D calculations 
are desirable. 

Existing results also indicate that muonic atoms in 
high-intensity, high-frequency laser fields can be utilized 
to dynamically gain structure information on nuclear 
ground states via their high-harmonic response. Besides, 
the laser-driven muonic charge cloud, causing a time- 
dependent Coulomb field, may lead to nuclear excita- 
tion (Shahbaz et ai, 2009). The excitation probabilities 



XII. LASER COLLIDERS 

The fast advancement in laser technology is opening up 
the possibility of employing intense laser beams to effi- 
ciently accelerate charged particles and to make them col- 
lide for eventually even initiating high-energy reactions. 



A. Laser acceleration 

Strong laser fields provide new mechanisms for par- 
ticle acceleration alternative to conventional accelerator 
technology (Tajima and Dawson, 1979). With presently- 
available laser systems an enormous electron acceleration 
gradient ~ 1 GeV/cm can be achieved, which exceeds by 
three orders of magnitude that of conventional accelera- 
tors and which raises prospects for compact accelerators 
(Clayton et ai, 2010; Faure et ai, 2004; Geddes et al, 
2004; Hafz et al, 2008; Mangles et ai, 2004). Differ- 
ent schemes of laser-electron acceleration have been pro- 
posed. These include the laser wakefield accelerator, the 
plasma beat wave accelerator, the self-modulated laser 
wakefield accelerator and plasma waves driven by mul- 
tiple laser pulses (see the recent reviews Esarey et at, 
2009 and Malka, 2011). High-gradient plasma wakefields 
can also be generated with an ultrashort bunch of pro- 
tons (Caldwell et al., 2009), allowing electron acceler- 
ation to TeV energies in a single stage. The achiev- 
able current and emittance of presently-available laser- 
accelerated electron beams is sufficient to build syn- 
chrotron radiation sources or even to aim at compact 
XFEL lasers (Schlenvoigt et ai, 2008). 

Laser acceleration of ions provides quasimonoenergetic 
beams with energy of several MeVs per nucleon (Fuchs 
et ai, 2007; Haberberger et ai, 2012; Hegelich et ai, 
2006; Schwoerer et al, 2006; Toncian et al., 2006). It 
mostly employs the interaction of high-intensity lasers 
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with solid targets. One of the main goals of laser-ion 
acceleration is to create low-cost devices for medical ap- 
plications, such as for hadron cancer therapy (Combs 
et at, 2009). Several regimes have been identified for 
laser-ion acceleration (see also the forthcoming review 
Macchi et ai, 2012). For laser intensities in the range 
10^^-10^^ W/cm^ and for solid targets with a thickness 
ranging from a few to tens of micrometers, the so-called 
target-normal-sheath acceleration is the main mecha- 
nism (Fuchs et ai, 2006). A further laser- ion interac- 
tion process is the skin-layer ponderomotive acceleration 
(Badziak, 2007). By contrast, the radiation-pressure ac- 
celeration regime operates when the target thickness is 
decreased (see Esirkepov et al., 2004 and Macchi et al., 
2009 for the so-called "laser piston" and "light sail" 
regimes, respectively). In Galow et al., 2011 a chirped 
ultrastrong laser pulse is applied to proton acceleration 
in a plasma. Chirping of the laser pulse ensures optimal 
phase synchronization of the protons with the laser field 
and leads to efficient proton energy gain from the field. 
In this way, a dense proton beam (with about lO'' protons 
per bunch) of high energy (250 MeV) and good quality 
(energy spread ~ 1%) can be generated. 

An alternative promising way for particle acceleration 
is direct laser acceleration in a tightly focused laser beam 
(Salamin and Keitel, 2002) or in crossed laser beams 
(Salamin et al., 2003). Especially efficient accelerations 
(Bochkarev et al., 2011; Gupta et al., 2007; Salamin, 
2006) can be achieved in a radially polarized axicon laser 
beam (Dorn et al., 2003). For example, the generation 
of mono-encrgctic GcV electrons from ionization in a ra- 
dially polarized laser beam is theoretically demonstrated 
in Salamin, 2007, 2010. A setup for direct laser acceler- 
ation of protons and bare carbon nuclei is considered in 
Salamin et al., 2008. It has been shown that laser pulses 
of 0. 1-10 PW can accelerate the nuclei directly to energies 
in the range required for hadron therapy. Simulations 
in Galow et al., 2010 and further optimization studies in 
Harman et al., 2011 indicate that protons stemming from 
laser-plasma processes can be efficiently post-accelerated 
employing single and crossed pulsed laser beams, focused 
to spot radii of the order of the laser wavelength. The 
protons in the resulting beam have kinetic energies ex- 
ceeding 200 MeV and small energy spreads of about 1%. 
The direct-acceleration method has proved to be efficient 
also for other applications. In Salamin, 2011 it is shown 
that 10 keV helium and carbon ions, injected into 1 TW- 
power crossed laser beams of radial polarization, can be 
accelerated in vacuum to energies of hundreds of kcV 
necessary for ion lithography. 



B. Laser-plasma linear collider 

Laser-electron accelerators have already entered the 
GeV energy domain where the realm of particle physics 



starts (Leemans et al., 2006). In fact, the strong interac- 
tion comes into play at distances d of the order of d ~ 1 
fm, which for relativistic processes corresponds to ener- 
gies e ^ 1/d ~ 1 GeV. Thus, a laser-based collider is in 
principle suitable for performing particle-physics exper- 
iments. However, in order to initiate high-energy reac- 
tions with sizable yield, not only GeV energies are re- 
quired but also collision luminosities C at least as high 
as 10^^-10^^ cm^^s^^. Meanwhile, for the ultimate goal 
of being competitive with the next International Linear 
Collider (ILC, 2011), energies on the order of 1 TeV and 
luminosities of the order of 10'^'^ cm^^s^^ are required 
(EUis and Wilson, 2001). 

The potential of the Laser-Plasma Accelerator (LPA) 
scheme to develop a laser-plasma linear collider is dis- 
cussed in Schroeder et al., 2010. Two LPA regimes are 
analyzed which are distinguished by the relationship be- 
tween the laser beam waist size wq and the plasma fre- 
quency ujp = Aimpe^ /m, with Up being the plasma 
density: 1) the quasilinear regime at large radius of the 
laser beam ujIwI > 2^l/^/T+£U2 (^o ^ 1); 2) the bub- 
ble regime, at ujpWo < 2v^ (^o > !)• The latter is 
less suitable for the collider application because the bub- 
ble cavity leads to dcfocusing for positrons; besides, the 
focusing forces and accelerating forces are not indepen- 
dently controllable here as both depend on the plasma 
density. Whereas, in the quasilinear regime this is pos- 
sible due to the existence of a second control parameter 
given by the laser beam waist size. In the following, 
we will discuss qualitatively the scaling properties of the 
quasilinear regime. For more accurate expressions the 
reader is referred to Schroeder et ai, 2010. 

In the standard LPA scheme, the electron plasma wave 
is driven by an intense laser pulse with duration tq of 
the order of the plasma wavelength Xp = 27r/wp, which 
accelerates the electrons injected in the plasma wave by 
wave breaking (Esarey et ai, 2009; Leemans and Esarey, 
2009; Malka et al., 2008). The accelerating field Ep of 
the plasma wave can be estimated from 

muj„ I/O , , 

Ep^-^^ ny^. (55) 

In fact, in the plasma wave, the charge separation occurs 
on a length scale of the order of \p, producing a sur- 
face charge density Op ^ \e\np\p and a field Ep ^ 47rcrp, 
which corresponds to Eq. (55). The number of elec- 
trons A'e that can be accelerated in a plasma wave is 
given approximately by the number of charged particles 
required to compensate for the laser-excited wakeficld, 
having a longitudinal component E\\ . From the relation 
^ 47rAfe|e|/7rwg, it follows that 

iVe - ^ (X (56) 

because WQiOp '-^ 1 in the quasilinear regime. The inter- 
action length of a single LPA stage is limited by laser- 
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diffraction effects, dephasing of the electrons with re- 
spect to the accelerating field, and laser-energy deple- 
tion. Laser-diffraction effects can be reduced by employ- 
ing a plasma channel, and plasma-density tapering can 
be utilized to prevent dephasing. Therefore, the LPA 
interaction length will be determined by the energy de- 
pletion length Ld- We can estimate the latter by equat- 
ing the energy spent for accelerating the Ng electrons 
along Ld Ne\e\EpLd) to the energy of the laser pulse 
(~ EqItwItq / 8tt) . Recalling that woOJp ^ 1, this yields 

Ld-4^pOcn;y\ (57) 

A staging of LPA is required to achieve high current den- 
sities along with high energies. The electron energy gain 
Aes in a single-stage LPA is 

Ae, - \e\EpLd - m-| cx n'^. (58) 

Therefore, the number of stages Ng to achieve a 
total acceleration energy Eq is A^^ = Eo/Aes ^ 
{Eo/m){u!p/ujQ)'^ cx np. It corresponds to a total collider 
length Lc of 

L,^ NsLd^-Xp(xnZ^/\ (59) 

When two identical beams each with N particles and 
with horizontal (vertical) transverse beam size i'^y) 
collide with a frequency /, the luminosity C is defined as 
C = N^f/AnaxCTy. In LPA N = Ne and / is the laser 
repetition rate, then 



647r Lj^axCTy 

The laser energy Ws required in a single stage in 
the LPA collider is approximately given by Wg ^ 
m{Xp/rQ){u!o/u!p)'^ at ^ 1 and the total required power 
Pt amounts to Pt NsWsf ~ eof{Xp/ro). 

The above estimates show that, although the number 
Ae of electrons in the bunch as well as the single-stage 
energy gain Ae^ increase at low plasma densities, the 
accelerating gradient Asg/Ld nevertheless decreases be- 
cause the laser depletion length Ld and the overall col- 
lider length Lc increase as well. Limiting the total length 
of each LPA in a collider to about 100 m will require 
a plasma density Up ~ 10^^ cm^"^ to provide a center- 
of mass energy of ~ 1 TeV for electrons, with ~ 10 
GeV energy gain per stage (Schroeder et ai, 2010). For 
a number of electrons per bunch of Ag ~ 10^, a laser 
repetition rate of 15 kHz and a transverse beam size of 
about 10 nm would be required to reach the goal-value 
of 10'^^ cm^^s^^ for the accelerator luminosity. At the 
usual condition wq ^ 1/ijJp ^ 10 /xm, instead, the lu- 
minosity amounts to £ ^ 10^^ cm~^s~^. In the above 



conditions each acceleration stage would be powered by 
a laser pulse with an energy of 30 J corresponding to 
an average power of about 0.5 MW at the required rep- 
etition rate of 15 kHz. Such high-average powers are 
beyond the performance of present-day lasers. The fu- 
ture hopes for high-average-power lasers are connected 
with diode-pumped lasers and new amplifier materials. 
Further challenges of LPA colliders are their complexity, 
as they involve plasma channels and density tapering, 
and, most importantly, the problem of how to accelerate 
positrons. 



C. Laser micro-collider 

Wc turn now to another scheme for a laser collider, 
which is based on principles quite different than those of 
the LPA scheme discussed above. In the LPA scheme, 
the electron is accelerated due to its synchronous motion 
with the propagating field. In this way the symmetry 
in the energy exchange process between the electron and 
the oscillating field is broken, as required by the Lawson- 
Woodward theorem (Lawson, 1979; Woodward, 1947). 
Another way to exploit the energy gain of the electron 
in the oscillating laser field is to initiate high-energy pro- 
cesses in situ, i.e., inside the laser beam (McDonald and 
Shmakov, 1999). In this case the temporary energy gain 
of the electron during interaction with a half cycle of the 
laser wave is used to trigger some processes during which 
the electron state may change (in particular, the electron 
may annihilate with a positron) and the desired asymme- 
try in the energy exchange can be achieved. In fact, this 
approach is widely employed in the nonrelativistic regime 
via laser-driven recollisions of an ionized electron with its 
parent ion (see Sec. IV. B). 

The question arises also as to whether the temporary 
energy gain of the electron in the laser beam can also be 
employed in the relativistic regime at ultrahigh energies. 
As pointed out in Sec. IV, an extension of the established 
recoUision scheme with normal atoms into the relativistic 
regime is hindered by the relativistic drift. However, the 
drift will not cause any problem when positronium (Ps) 
atoms arc used because its constituent particles, electron 
and positron, have the same absolute value of the charge- 
to-mass ratio (see Sec. III. A). 

A corresponding realization of high-energy e"^-e~ rec- 
ollisions in the GeV domain aiming at particle reactions 
has been proposed in Hatsagortsyan et at, 2006. It re- 
lies on (initially nonrelativistic) Ps atoms exposed to 
super-intense laser pulses. After almost instantaneous 
ionization of Ps in the strong laser field, the free electron 
and positron oscillate in opposite directions along the 
laser electric field and experience the same ponderomo- 
tive drift motion along the laser propagation direction. In 
this way, the particles acquire energy from the field and 
arc driven into periodic e'^-e~ collisions (Henrich et at, 
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FIG. 25 (a) In conventional e+-e^ colliders bunches of ac- 
celerated electrons and positrons are focused to collide head- 
on-head incoherently, i.e., the bunches collide head-on-head 
but electrons and positrons in the bunch do not. (b) In the 
recoUision-based collider, the electron and positron originat- 
ing from the same Ps atom may collide head-on-head coher- 
ently (Henrich et al., 2004). From Hatsagortsyan et al., 2006. 



2004). Provided that the applied laser intensity is large 
enough, elementary particle reactions like heavy lepton- 
pair production can be induced in these recoUisions. The 
common center-of-mass energy s* of the electron and the 
positron at the recollision time arises mainly from the 
transversal momentum of the particles and it scales as 
e* ~ ni^o- A basic particle reaction which could be trig- 
gered in a laser-driven collider is e+-e^ annihilation with 
production of a pair, i.e., e~^e~ — > The 

energy threshold for this process in the center-of-mass 
system is 2m^ w 210 MeV. It can be reached with a laser 
field such that ^ 200, corresponding to laser intensities 
of the order of 10^^ W/cm^, currently within reach. 

In addition, the proposed recollision-based laser col- 
lider can yield high luminosities compared to conven- 
tional laser accelerators. In the latter, bunches of elec- 
trons and positrons are accelerated and brought into 
head-on-head collision. However, the particles in the 
bunch are distributed randomly such that each micro- 
scopic e'^-e~ collision is not head-on-head but has a mean 
impact parameter bi ^ ab determined by the beam ra- 
dius flb, characterizing the collision as incoherent (see 
Fig. 25). Instead, in the recollision-based collider the 
electron and the positron stem from the same Ps atom 
with initial coordinates being confined within the range 
of one Bohr radius as « 5.3 x 10~^ cm. Since they are 
driven coherently by the laser field, they can recollide 
with a mean impact parameter be a^p of the order of 
the electron wave packet size a^p (see Fig. 25). Con- 
sequently, the luminosity contains a coherent component 
(Hatsagortsyan et al., 2006): 



Np{Np-l) Np 
5? ^ 6? 



/, 



(61) 



where Np is the number of particles in the bunch, and / 
is the bunch repetition frequency. The coherent com- 
ponent {Np/b1)f can lead to a substantial luminosity 
enhancement in the case when the particle number is 
low and the particle's wave packet spreading is small, 
NpO^p < a^. Note that the reaction Ps — fJ,~^fJ,~ arising 



in a strongly laser-driven e^-e" plasma may be consid- 
ered as the coherent counterpart of the incoherent process 
e+e^ — > /i^/i^ discussed in Sec. X.A. 

Rigorous quantum-electrodynamical calculations for 
pair production in a laser field have been per- 
formed in Miiller et al., 2006, 2008b, c. In agreement 
with Eq. (61), they enabled the development of a simple- 
man's model in which the rate of the laser-driven pro- 
cess can be expressed via a convolution of the rescat- 
tering electron wave packet with the field-free cross sec- 
tion CTe+e^-^u+u- (see Eq. 51). The latter attains the 



maximal value a 



(max) 
' e+e^ 



at 



£* sa 260 MeV (Peskin and Schroeder, 1995). However, 
when the field driving the Ps atoms is a single laser wave, 
the e^-e" recollision times are long and the pro- 
duction process is substantially suppressed by extensive 
wave packet spreading. This obstacle can be overcome 
when two counterpropagating laser beams are employed. 

The role of the spreading of the electron wave packet in 
counterpropagating focused laser beams of circular and 
linear polarization has been investigated in detail in Liu 
et al., 2009. The advantage of the circular-polarization 
setup is the focusing of the recoUiding electron wave 
packet. However, this advantage is reduced by a spa- 
tial offset in the e'^-e~ collision when the initial coordi- 
nate of the Ps atom deviates from the symmetric position 
between the laser pulses. The latter imposes a severe re- 
striction on the Ps gas size along the laser propagation 
direction. Thus, the linear-polarization setup is prefer- 
able when the offset at the recollision is very small and 
the wave packet size at the recollision is within accept- 
able limits. Results from a Monte Carlo simulation of 
the e"'"-e~ wave-packet dynamics in counterpropagating 
linearly polarized laser pulses are shown in Fig. 26. 

The luminosity C and the number of reaction 
events M for the recollision-based collider with counter- 
propagating laser pulses can be estimated as: 

C^Np,-^Trf, (62) 



wp 
(max) 



-TrNp^Nr 



(63) 



respectively, where Nps is the number of Ps atoms, 
the recollision time of the order of the lasers period, Nl 
the number of laser pulses, and / the laser repetition 
rate. Taking Np^ « 10* (Cassidy and Mills, 2005), / = 1 
Hz and the spatial extension of the e+-e~ pair from Fig. 
26, one estimates a luminosity of £ '--^ 10^^ cm~^s~^ and 
about one fJ,~^-n~ pair production event every 10'^ laser 
shots at a laser intensity of 4.7 x 10^^ W/cm^. 

In conclusion, the scheme of the recollision-based laser 
collider allows to realize high-energy and high-luminosity 
collisions in a microscopic setup. However, it is not easily 
scalable to the parameters of the ILC, namely, to TeV 
energies and luminosities of the order of 10"^ 
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FIG. 26 (Color) The coordinate-space distributions of the 
electron and the positron wave-packets at the recoUision time 
in focused counterpropagating pulses along the z direction 
with Wo = 10 ^m, Ao = 0.8 and with 7o = 4.7 x 
lO^^W/cm^ (parts (a) and (c)) and Jo = 1.4 x lO^^W/cm^ 
(parts (b) and (d)). The Ps atom is initially located at the 
origin. Spatial coordinates are given in "atomic units", with 
1 a.u. = 0.05 nm. Adapted from Liu et al, 2009. 



XIII. PARTICLE PHYSICS WITHIN AND BEYOND THE 
STANDARD MODEL 

The sustained progress in laser technology towards 
higher and higher field intensities raises the question as 
to what extent ultrastrong laser fields may develop into 
a useful tool for particle physics beyond QED. Below we 
review theoretical predictions regarding the influence of 
super-intense laser waves on electroweak processes and 
their potential for probing new physics beyond the Stan- 
dard Model. 



A. Electroweak sector of the Standard Model 

The energy scale of weak interactions is set by the 
masses of the and exchange bosons, mw ~ mz ^ 
100 GeV. Therefore, the influence of external laser fields, 
even if strong on the scale of QED, is generally rather 
small. An overview of weak interaction processes in the 
presence of intense electromagnetic fields has been given 
in Kurflin, 1999. 

Various weak decay processes in the presence of in- 
tense laser fields have been considered. They can be di- 
vided into two classes: 1) laser-assisted processes which 
also exist in the absence of the field but may be modi- 
fied due to its presence; 2) field-induced processes which 
can only proceed when a background field is present, 
providing an additional energy reservoir. With respect 
to processes from the first category, tt — > /i -t- and 
fi~ ^ e~ + Ve + have already been examined (Ritus, 
1985). Laser-assisted muon decay has also been revis- 
ited recently (Dicus et ai, 2009; Farzinnia et ai, 2009; 
Narozhny and Fedotov, 2008). and boson decay 



into a fermion-antifermion pair was calculated in Kurilin, 
2004, 2009. In all cases, the effect of the laser field was 
found to be small. As a general result, the presence of the 
laser field modifies the field- free decay rate Rm,o of a par- 
ticle with mass M to Rm = -Ro,m(1 + A), with the correc- 
tion A being of the order of Xo m range of param- 
eters ^o.M = iom/M > 1 and xo,m = Xo{m/M)^ < 1. 
The elastic scattering of a muon neutrino and an electron 
in the presence of a strong laser field has been considered 
in Bai et ai, 2012 and multiphoton effects in the cross 
section are predicted. 

External fields can also induce decay processes, which 
are energetically forbidden otherwise. In Kurilin, 1999, 
the field-induced lepton decay l~ — >■ W~ + vi was 
considered. Since the mass of the initial-state parti- 
cle is smaller than the mass of the decay products, 
the process is clearly impossible in vacuum. The pres- 
ence of the field does allow for such an exotic decay, 
but the probability Pi- +ui remains exponentially 
suppressed, i.e., Pi-^w-+ui ^ cxp(— l/xo,iy), where 
Xo,w = Xoi'm/mwf- 

Finally, the production of an e+-e~ pair by high-energy 
neutrino impact on a strong laser pulse has been cal- 
culated in Tinsley, 2005. The setup is similar to the 
e^-e~ pair production processes in QED discussed in 
Sees. VIII. A and VIII. B. However, as it was shown in 
Tinsley, 2005, the laser-induced process v ^ u + e^ + 
is extremely unlikely. At a field intensity of about 
/o = 3 X 10^® W/cm^, the production length is on the 
order of a light year, even for a neutrino energy of 1 PeV. 

B. Particle physics beyond the Standard Model 

Recently, a lot of attention has been devoted to the 
possibility of employing intense laser sources to test as- 
pects of physical theories which go even beyond the Stan- 
dard Model. In Heinzl et al, 2010b, for example, it is 
envisaged that effects of the noncommutativity of space- 
time modify the kinematics of multiphoton Compton 
scattering by inducing a nonzero photon mass. We recall 
that in noncommutative quantum field theories operators 
X^^ are associated to spacetime coordinates x^, which do 
not commute, but rather satisfy the commutation rela- 
tions [AT^jX*^] = iQ^'^ , with B'^" being an antisymmetric 
constant tensor (Douglas and Nekrasov, 2001). 

On a different side, one of the still open problems of 
the Standard Model is the so-called strong CP problem 
(Kim and Carosi, 2010). The nontrivial structure of 
the vacuum, as predicted by Quantum Chromodynam- 
ics (QCD), allows for the violation within QCD of the 
combined symmetry of charge conjugation (C) and par- 
ity (P). This implies a value for the neutron's electric 
dipole moment which, however, is already many orders 
of magnitude larger than experimental upper limits. One 
way of solving this problem was suggested in Peccei and 
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Quinn, 1977 which required the existence of a massive 
pseudoscalar boson, called axion. The axion has never 
been observed experimentally although some of its prop- 
erties can be predicted on theoretical grounds: it should 
be electrically neutral and its mass should not exceed 
1 eV in order of magnitude. Although being electri- 
cally neutral, the axion is predicted to couple to the elec- 
tromagnetic field T^'^{x) through a Lagrangian-density 
term 

^a^{x) = ^aix)T'"'{x)f^,{x), (64) 

with g being the photon- axion coupling constant and a{x) 
the axion field. 

The photon-axion Lagrangian density in Eq. (64) has 
mainly two implications: 1) the existence of axions in- 
duces a change in the polarization of a light beam pass- 
ing through a background electromagnetic field; 2) a 
photon can transform into an axion (and vice versa) 
in the presence of a background electromagnetic field. 
The first prediction has been tested in experiments 
like the Brookhaven-Fcrmilab- Rochester- Trieste (BFRT) 
(Cameron et ai, 1993) and the Polarizzazione del Vuoto 
con LAScr (PVLAS) (PVLAS, 2011), where a linearly 
polarized probe laser field crossed a region in which a 
strong magnetic field was present of 3.25 T and 5 T 
at BFRT and at PVLAS, respectively. Testing the sec- 
ond prediction is the aim of the so-called "light shining 
through a wall" experiments like the Any Light Particle 
Search (ALPS) (ALPS, 2010), the CERN Axion Solar 
Telescope (CAST) (CAST, 2008) and Gamma to miUi- 
eV particle search (GammeV) (GammeV, 2011) (see also 
the detailed theoretical analysis in Adler et ai, 2008). 
In the GammeV experiment, for example, the light of a 
Nd:YAG laser passes through a region in which a 5 T 
magnetic field is present. A mirror is positioned behind 
that region in order to reflect the laser light. The axions 
which would eventually be created in the magnetic-field 
region pass through the mirror undisturbed and can be 
reconverted to photons by means of a second magnetic 
field, activated after the mirror itself. So far these exper- 
iments have given negative results. An interesting exper- 
imental proposal has been put forward in Rabadan et al. , 
2006, where the high-energy photon beam delivered by 
an XFEL facility has been suggested as a probe beam to 
test regions of parameters (like the axion mass nia or the 
photon-axion coupling constant g) which are inaccessible 
via optical laser light. 

The perspective for reaching ultrahigh intensities at 
future laser facilities has stimulated new proposals for 
employing such fields to elicit the photon-axion interac- 
tion (Gies, 2009). In fact, an advantage of using strong 
uniform magnetic fields is that they can be kept strong 
for a macroscopically long time (of the order of hours) 
and on a macroscopic spatial region (of the order of 1 m) 
(for QED processes occurring in a strong magnetic field. 



see the standard review Erber. 1966 and the very current 
overview paper Dunne, 2012 for recent progresses in the 
field). On the other hand, laser beams deliver fields much 
stronger than those employed in the mentioned experi- 
ments (the magnetic field strength of a laser beam with 
the available intensity of 10^^ W/cm^ amounts to about 
6.5 X 10^ T) but in a microscopic space-time region. How- 
ever, it has been first realized in Mendonga, 2007 that en- 
visaged ultrahigh intensities at future laser facilities may 
compensate for the tiny space-time extension of the laser 
spot region. In Mendonga, 2007 the coupled equations of 
the electromagnetic field J'^''{x) and the axion field a{x) 

d^d^a^mla=\gT^^f,,^ ^^^^ 

are solved approximately. It is shown that if a probe 
laser field propagates through a strong plane- wave field, 
the axion field "grows" at the expense mainly of the 
probe field itself, whose intensity should be observed to 
decrease. Laser powers of the order of 1 PW have al- 
ready been shown to provide stronger hints for the pres- 
ence of axions than magnetic-field-based experiments like 
the PVLAS. More realistic Gaussian laser beams are 
considered in Dobrich and Gies, 2010 where the start- 
ing point is also represented by Eq. (65). The sug- 
gested experimental setup assumes a probe electromag- 
netic beam with angular frequency Wp^in passing through 
a strong counterpropagating Gaussian beam with angu- 
lar frequency Wq,!! ^-iid another strong Gaussian beam 
propagating perpendicularly and with angular frequency 
a;o,±. By choosing wo,_l = 2^0, y, it is found that after 
a photon-axion-photon double conversion, photons are 
generated with angular frequencies Wp,out = ^p,i-a ± "^o.ll- 
The amplitudes of these processes are shown to be 
peaked at specific values of the axion mass ma^± = 

2 y/wp^nWoJI + '^o ll^-*- -'-)/2- Since the optical photon 
energies are of the order of 1 eV, this setup allows for the 
investigation of values of the axion mass in this regime. 
This is very important because such a region of the ax- 
ion mass is inaccessible to experiments based on strong 
magnetic fields, which can probe regions at most in the 
mcV range. 

In addition to electrically neutral new particles such 
as axions, yet unobserved particles with nonzero charge 
may also exist. The fact that they have so far escaped 
detection implies that they are either very heavy (ren- 
dering them a target for large-scale accelerator experi- 
ments), or that they are light but very weakly charged. 
In the latter case, these so-called minicharged particles, 
i.e., particles with absolute value of the electric charge 
much smaller than |e|, are suitable candidates for laser- 
based searches (Gies, 2009). Let and Qc = ee, with 
< e <C 1; denote the minicharged particle mass and 
charge, respectively. Then, the corresponding critical 
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field scale fcr,e = 'm1/\Qf\ can be much lower than 
-Fcr- As a consequence, vacuum nonlinearities associ- 
ated with minicharged particles may be very pronounced 
in an external laser field with intensity much less than 
lev ~ 10-^^ W/cm^. Moreover, even at optical photon 
energies ~ 1 eV the efi'ective Lagrangian approach might 
become inappropriate to describe the relevant physics if 
m, < 1 eV (see Sec. VILA). In Gies et al, 2006, vac- 
uum dichroism and birefringence effects due to the ex- 
istence of minicharged particles were analyzed when a 
probe laser beam with ujp > 2me traverses a magnetic 
field. It was shown that polarization measurements in 
this setup would provide much stronger constraints on 
minicharged particles in the mass range below 0.1 eV 
than in previous laboratory searches. 

XIV. CONCLUSION AND OUTLOOK 

The fast development of laser technology has been 
paving the way to employ laser sources for investigating 
relativistic, quantum electrodynamical, nuclear and high- 
energy processes. Starting with the lowest required in- 
tensities, relativistic atomic processes are already within 
the reach of available laser systems, while the proposed 
methods to compensate the deteriorating effects of the 
relativistic drift still have to be tested experimentally. 
Moreover, a fully consistent theoretical interpretation of 
recent experimental results on correlation effects in rela- 
tivistic multielectron tunneling is still missing. 

Concerning the interaction of free electrons with in- 
tense laser beams, we have seen that experiments have 
been performed to explore the classical regime. Only 
the E-144 experiment at SLAC has so far been realized 
on multiphoton Compton scattering, although presently 
available lasers and electron beams would allow for prob- 
ing this regime in full detail. We have also pointed out 
that at laser intensities of the order of 10^^-10^'^ W/cm^, 
RR effects come into play at electron energies of the order 
of a few GeV. It is envisaged that the quantum radiation- 
dominated regime, where quantum and RR effects sub- 
stantially alter the electron dynamics, could be one of 
the first extreme regimes of light-matter interaction to 
be probed with upcoming petawatt laser facilities. On 
the theoretical side, most of the calculations have been 
performed by approximating the laser field as a plane 
wave, as the Dirac equation in the presence of a focused 
background field cannot be solved analytically. Certainly, 
new methods have to be developed to calculate photon 
spectra including quantum effects and spatio-temporal 
focusing of the laser field in order to be able to quanti- 
tatively interpret upcoming experimental results. 

Nonlinear quantum electrodynamical effects have been 
shown to become observable at future multipetawatt 
laser facilities, as well as at ELI and HiPER. Here the 
main challenges concern the mcasurability of tiny effects 



on the polarization of probe beams and on the detection 
of a typically very low number of signal photons out of 
large backgrounds. Similar challenges are envisaged to 
detect the presence of light and weakly-interacting hypo- 
thetical particles like axions and minicharged particles. 
The physical properties (mass, coupling constants, etc.) 
of such hypothetical particles are, of course, unknown. 
In this respect intense laser fields may be employed here 
to set bounds on physical quantities like the axion mass 
and, in particular, to scan regions of physical parameters, 
which are inaccessible to conventional methods based, for 
example, on astrophysical observations. 

Different schemes have been proposed to observe e"''-e~ 
pair production at intensities below the Schwinger limit, 
which seems now to be feasible in the near future, at least 
from a theoretical point of view. Corresponding studies 
would complement the results of the pioneering E-144 
experiment and deepen our understanding of the QED 
vacuum in the presence of extreme electromagnetic fields. 
This is also connected with the recent investigations on 
the development of QED cascades in laser-laser collisions. 
In addition to being intrinsically interesting, the devel- 
opment of QED cascades is expected to set a limit on the 
maximal attainable laser intensity. However, the study of 
quantum cascades in intense laser fields has started rela- 
tively recently and is still under vivid development. More 
advanced analytical and numerical methods are required 
in order to describe realistically and quantitatively such 
a complex system as an electron-positron-photon plasma 
in the presence of a strong driving laser field. 

Nuclear quantum optics is also a new exciting and 
promising field. Since the energy difference between nu- 
clear levels is typically in the multi-keV and MeV range, 
high-frequency laser pulses, especially in combination 
with accelerators, are preferable in controlling nuclear 
dynamics. As pointed out, especially table-top highly co- 
herent x-ray light beams, envisaged for the future, open 
up perspectives for exciting applications including nu- 
clear state preparation and nuclear batteries. 

Finally, we want to point out that most of the consid- 
ered processes have not yet been observed or tested ex- 
perimentally. This is in our opinion one of the most chal- 
lenging aspects of upcoming laser physics, not only from 
an experimental point of view, but also from the point of 
view of theoretical methods. Experimentally, the main 
reason is that in order to test, for example, nonlinear 
quantum electrodynamics or to investigate nuclear quan- 
tum optics, high-energy particle beams (including photon 
beams) are required to be available in the same labora- 
tory as the strong laser. On the one hand, the combined 
expertise from different experimental physical communi- 
ties is required to perform such complex but fundamental 
experiments. On the other hand, the fast technological 
development of laser-plasma accelerators is very promis- 
ing and exciting, as this seems the most feasible way 
towards the realization of stable, table-top high-energy 
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particle accelerators. Combining such high-energy probe 
beams with an ultrarelativistic laser beam in a single all- 
optical setup will certainly result in a unique tool for 
advancing our understanding of intense laser-matter in- 
teractions. 



LIST OF FREQUENTLY-USED SYMBOLS 



Eo 

F,,^my\e\ 



I,, = F^,/4n 

Z 

e 

fc^" = (w, k) = unt" 

k'l" = {uj',k') =Lo'n''' 

fcg = ujoriQ = a;o(l, no) 
m 

Po = (£o,Po) = "^7o(l,/3o) 



Xo = {po,-/m){EQ/Fcr) 
= {k^/m){Eo/F„) 

Ao = To = 27r/wo 
Ac = l/m 

£.0 = \e\EQ/muja 



laser field amplitude 

critical electromagnetic field 

(« 1.3 X IQi^ V/cm) 

(« 4.4 X 10" G) 

laser peak intensity 

critical laser intensity 

(« 4.6 X 10^9 W/cm2) 

nuclear charge number 

electron charge 

(w -1/VT37« -0.085) 

initial or incoming photon 

four- momentum 

final or outgoing photon 

four- momentum 

laser photon four-momentum 

electron mass 

(« 0.511 MeV) 

initial or incoming electron 

four- momentum 

fine-structure constant 

(« 1/137 « 7.3 X 10-3) 

nonlinear electron quantum 

parameter 

nonlinear photon quantum 
parameter 

laser wavelength and period 
Compton wavelength 
(« 3.9 X 10-" cm) 
classical relativistic 
parameter 

laser angular frequency 
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